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OBSERVATIONS  OF  OH  RADIO  EMISSION  SOURCES 


ABSTRACT 

The  1 8 -cm  OH  radio  emission  sources  were  studied  to  determine  their  character 
and  their  relation  to  other  objects  in  our  galaxy.  The  Haystack  and  Millstone  an¬ 
tennas  were  used  with  a  100-channel  digital  autocorrelator  to  obtain  Stokes  param¬ 
eter  spectra  of  many  of  the  brighter  sources.  Time  variations  in  the  OH  emission 
from  some  of  these  sources  were  detected  and  studied  and  maps  of  the  regions 
around  OH  emission  sources  were  obtained  in  the  continuum  near  2  and  4  cm.  The 
Haystack -Millstone  spectral -line  interferometer  was  used  to  check  the  positions 
of  a  pair  of  OH  emission  sources  and  to  search  for  angular  variation  in  the  OH 
absorption  across  Cas  A.  Evidence  for  phenomenological  associations  between 
OH  emission  and  H  II  regions,  supernovae  remnants,  and  infrared  objects  is  dis¬ 
cussed  in  the  report  and  a  kinematic  model  of  the  OH  emission  sources  presented. 
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OBSERVATIONS  OF  OH  RADIO  EMISSION  SOURCES 


L  INTRODUCTION 
A.  OH  Molecule 

The  OH  molecule  is  an  electrically  neutral  combination  of  an  oxygen  and  a  hydrogen  atom. 

The  internal  energy  of  such  a  molecule  may  be  divided  into  three  types,  rotational  (i.e.f  energy 
of  rotation  of  the  constituent  atoms  about  their  barycenter),  vibrational  (i.e.,  energy  of  periodic 
variation  of  the  distance  between  the  two  atoms),  and  electronic  (i.e.,  energy  of  the  electrons 
bound  to  the  molecule).  In  molecules  such  as  OH,  transitions  between  rotational  states  give  rise 
to  microwave  and  far  infrared  lines,  transitions  between  vibrational  states  give  rise  to  near  in¬ 
frared  and  visible  lines,  and  transitions  between  electronic  states  give  rise  to  visible  and  ultra¬ 
violet  lines. 

The  ground  state  of  OH  is  ^  =  2  ^eve^  *s  higher  in  energy).  In  this  no¬ 

tation,  the  superscript  is  2S  +  1  where  S  is  the  unpaired  electron  spin  (there  are  nine  electrons, 
one  unpaired,  giving  S  =  l/2),  and  the  n  indicates  that  A,  the  projection  of  the  electronic  or¬ 
bital  angular  momentum  onto  the  internuclear  axis,  is  one.  The  subscript  following  the  n  is 
the  vector  sum  of  S  and  A,  i.e.,  the  total  electronic  angular  momentum  component  on  the  inter¬ 
nuclear  axis,  and  J  is  the  total  angular  momentum  (except  for  nuclear  spin),  i.e.,  the  sum  of 
the  electronic  and  rotational  angular  momenta. 

The  transition  from  the  ground  state  to  the  first  rotationally  excited  state  (  n_  J  =  5/2) 
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and  the  "fine  structure"  transitions  between  ^  states  and  II ^  ^  states  are  in  the  infrared. 

However,  the  ground  state  of  OH  is  A  doubled  with  a  separation  of  about  1666  MHz  in  the  radio 
microwave  region.  This  A  doubling  occurs  because  of  an  interaction  between  the  electronic  or¬ 
bital  angular  momentum  and  the  molecular  rotation,  and  may  be  thought  of  as  a  partial  uncoupling 
of  the  electronic  angular  momentum  from  the  internuclear  axis  with  a  tendency  to  couple  it  in¬ 
stead  to  the  rotation  axis.  The  ground  state  of  OH  is  intermediate  between  Hund's  coupling  case 
(a)  for  which  the  spin  and  orbital  electronic  angular  momenta  are  coupled  to  the  internuclear 
axis,  and  Hund's  case  (b)  for  which  the  orbital  electronic  angular  momentum  is  still  coupled  to 
the  internuclear  axis,  but  the  spin  is  coupled  instead  to  the  resultant  of  the  orbital  electronic 
and  rotational  angular  momenta.  A  doubling  represents  a  tendency  toward  Hund's  case  (d)  for 
which  the  orbital  angular  momentum  is  coupled  to  the  rotation  axis.  The  A  doublet  states  may 
properly  be  labeled  with  their  parity  (+  or  — )  but  not  by  A,  because  each  state  is  a  combination 
of  A  =  +1  and  —  1 . 

These  A  doublet  states  are  further  split  with  a  separation  of  about  53  MHz  by  an  interaction 
with  the  magnetic  moment  of  the  nuclear  spin.  These  substates  may  be  labeled  F,  the  total 
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angular  momentum  including  nuclear  spin.  The  resulting  energy  level  diagram  for  the  ground 
state  of  OH  is  shown  schematically  in  Fig.  1.  The  four  transitions  of  interest  are  at  1612231, 
1665401,  1667358,  and  1720527  kHz.  The  relative  line  strengths  (1 : 5: 9:1 )  are  also  indicated  in 
Fig.  1.  The  transitions  in  the  figure  are  electric  dipole  and  the  transition  probability  (the 
Einstein  A)  is  7.7  x  10  ^/sec  for  the  1667  transition. 

F'2 
FM 

Fig.  1.  Energy  level  diagram  of  OH. 
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Further  details  about  the  OH  molecule  may  be  found  in  Townes  and  Schawlow  (1955)  and 
Herzberg  (1  950). 

B.  History  of  OH  Radio  Astronomy 

The  original  proposal  that  the  radio  frequency  transitions  of  OH  should  be  looked  for  from 
extraterrestrial  sources  is  apparently  due  to  Shklovsky  (1953).  Townes  (1957)  also  made  such 
a  suggestion  and  in  1956  Barrett  and  Lilley  (1957)  tried  and  failed  to  detect  OH.  At  that  time, 
the  frequencies  of  the  lines  were  known  only  with  low  accuracy.  But  as  increasingly  accurate 
laboratory  measurements  became  available  [Dousmanis,  et  al.  (1955),  Ehrenstein,  et  al. .  (1959), 
Radford  (1961)  and  (1964)],  two  attempts  were  made  in  1963  [Penzias  (1964),  Weinreb,  et  al. 
(1963)]  Weinreb,  Barrett,  Meeks  and  Henry,  using  Lincoln  Laboratory's  84-foot  Millstone  an¬ 
tenna,  detected  the  1665-  and  1667-MHz  absorption  against  Cas  A  late  in  1963  [Weinreb,  et  al. 
(1963)]  This  was  confirmed  by  Dieter  and  Ewen  (1964)  and  Bolton,  et  al .  (1964),  and  Weaver 
and  Williams  (1964)  found  OH  absorption  against  Sag  A.  Barrett,  et  al .  (1964)  studied  Cas  A 
further  and  found  that  the  OH  absorption  near  zero  velocity  is  a  double  feature  without  a  counter¬ 
part  in  H  I  absorption  Robinson,  et  al.  (1964)  studied  the  Sag  A  OH  absorption  in  more  detail 
and  found  additional  strong  OH  absorption  with  almost  no  corresponding  feature  in  the  H  I  ab¬ 
sorption.  They  gave  evidence  for  believing  that  some  of  the  OH  is  not  distributed  uniformly  over 
the  source.  Gardner,  et  al.  (1964)  detected  OH  absorption  at  1612-  and  1720-MHz  against  Sag  A. 
Goldstein,  et  al.  (1964)  discovered  additional  OH  absorption  against  Sag  A,  again  without  any 
strong  counterpart  in  the  H  I  absorption. 

At  this  point  the  observations  seemed  to  fit  fairly  well  into  the  pattern  of  accepted  ideas 
about  the  galaxy.  The  ratios  of  the  absorptions  in  the  various  OH  lines  agreed  fairly  well  with 
the  values  predicted  by  the  transition  probabilities  with  the  possibility  of  a  high  optical  depth 
and  a  mottled  distribution  of  OH  [Gunderman  (1965)].  The  width  of  the  absorption  lines  ranges 
from  ten  to  several  hundred  kHz. 

In  1965  Gunderman  (1965)  and  Weaver,  Williams,  Dieter  and  Lum  [(1965)  and  (1966)]  dis¬ 
covered  remarkably  narrow  and  intense  emission  features  while  looking  for  absorption  against 
certain  Galactic  H  II  regions.  These  results  were  confirmed  [Zuckerman,  et_al.  (1965),  McGee, 
et  al.  (1965)]  and  it  was  found  that  some  of  the  features  are  strongly  linearly  [Weinreb,  et  al. 
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(1965)]  and  circularly  [Davies,  e_t  al.  (1966),  Barrett  and  Rogers  (1966)]  polarized.  The  polar¬ 
ization  properties  of  some  of  the  stronger  emission  sources  were  measured  in  detail  [Palmer 
and  Zuckerman  (1967),  Williams,  et  al.  (1966),  Meeks,  e_t  al .  (1966)]  and  it  was  found  that  strong 
circular  polarization  is  very  common,  linear  polarization  is  less  common,  and  completely  un¬ 
polarized  features  are  rare.  Very  little  correlation  of  features  between  1665  and  1667  MHz  was 
found  and  there  is  perhaps  evidence  of  anti-correlation.  No  simple  pattern  of  Zeeman  split  fea¬ 
tures  will  fit  the  data.  Although  the  kinetic  temperatures  implied  by  the  widths  of  the  emission 
features  are  typically  10°tol00°K,  the  brightness  temperatures  were  already  known  to  be  over 
1000°K  and  this  was  only  a  lower  limit  because  the  angular  sizes  of  the  sources  were  unresolved. 

To  add  to  the  list  of  anomalies.  Dieter,  Weaver,  and  Williams  (1966)  reported  large  time 
variations  within  a  span  of  weeks  in  some  of  the  OH  emission  features  coming  from  NGC  6334 
and  Orion,  and  these  results  have  been  at  least  partially  confirmed  [Palmer  and  Zuckerman 
(1967),  Ball  and  Meeks  (1968)],  although  this  is  clearly  an  area  which  needs  much  more  obser¬ 
vational  work. 

Because  the  angular  sizes  and  hence  the  brightness  temperatures  could  not  be  found  from 
single  antenna  measurements,  a  series  of  interferometer  measurements  were  undertaken 
[Cudaback,  e_t  al, .  (1966),  Rogers,  eR  al.  (1966),  Rogers,  et  al .  (1967),  Davies,  et  al .  (1967), 
Radhakrishnan  and  Whiteoak  (1967),  Moran,  et  al.  (1967  a  and  b)].  Most  of  the  OH  emission 
sources  do  not  exactly  coincide  with  strong  continuum  sources.  These  interferometer  measure¬ 
ments  have  succeeded  in  separating  some  of  the  sources  into  spatial  components,  but  some  of 
the  individual  features  have  not  yet  been  resolved  in  angle. 

The  latest  results  [Moran,  et,  al.  (1968)],  show  that  the  W3  OH  emission  is  coming  from  a 
region  a  few  thousand  au's  in  diameter.  Since  this  size  is  characteristic  of  solar-system  di¬ 
mensions  rather  than  interstellar  distances,  we  are  led  to  think  that  we  are  observing  OH  emis¬ 
sion  from  a  circumstellar  envelope  around  a  single  star,  or  perhaps  a  protostar.  Furthermore, 
within  this  region  the  individual  OH  emission  features  originate  in  much  smaller  spots  typically 
only  a  few  au  in  size.  The  apparent  brightness  temperature  of  one  of  the  features  in  W3  is  now 

known  to  be  at  least  10*2°K  [Moran,  et  ad.  (1968)]. 
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Meanwhile  other  OH  lines  were  sought.  O  H  was  seen  in  absorption  [Barrett  and  Rogers 
(1964),  Robinson  (1967)]  but  not  in  emission  [Wilson  and  Barrett  (1968)].  The  2n  ^  J  =  l/2 
transition  [Zuckerman,  et  al.  (1968)],  the  ^3/ 2 »  ^  =  5/2  transition  [Yen,  et_  al .  (1969)],  and  the 
2n  4  J  =  5/2  transition  (Schwarz,  P. ,  private  communication)  were  found  in  OH  emission 
The  recent  discovery  of  water  vapor  in  emission  from  some  of  the  same  sources  seen  in  OH 

emission  [Cheung,  et  al.  (1969),  Knowles,  et  al .  (1969)]  promises  to  help  us  understand  the  na- 

2  2  0 

ture  of  these  sources.  The  n  —  2  bands  of  OH  occur  around  3000  A  and  although  this  wave¬ 
length  is  rather  strongly  absorbed  by  the  earth's  atmosphere,  these  bands  have  been  looked  for 
in  stellar  spectra  [Goss  and  Spinrad  (1966)],  but  without  success. 

From  the  myriad  of  problems  concerning  the  OH  emission  sources,  the  following  questions 
were  selected  for  this  study,  primarily  because  the  technology  for  making  significant  measure¬ 
ments  was  available  or  in  prospect.  (1)  What  additional  observational  constraints  can  be  placed 
on  the  models  of  the  OH  emission  sources  and  the  kinematic  and  spatial  relationships  within  and 
among  the  sources9  To  approach  this  question,  the  details  of  the  Stokes  parameters  of  the  OH 
emission  features  from  some  of  these  sources  were  studied  and  time  variations  in  these  param¬ 
eters  were  monitored.  (2)  What  associations  and  interactions  exist  between  the  OH  emission 
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sources  and  nearby  objects  (H  II  regions,  SNR's,  etc.)?  To  approach  this  problem,  the  regions 
around  some  of  the  OH  emission  sources  were  mapped  in  the  radio  continuum  near  2  and  4cm. 

The  material  presented  here  can  only  be  considered  a  progress  report;  some  progress  has 
been  made  on  these  problems,  but  some  new  problems  have  been  uncovered.  Sections  II  through 
IV  present  the  data  acquisition  and  data  reduction  techniques  used  in  this  study.  Some  of  these 
techniques  are  new  and  not  previously  discussed  in  print;  however,  the  reader  who  is  primarily 
interested  in  the  astronomical  results  may  wish  to  skip  to  Sec.  V. 

n.  EQUIPMENT  -  INSTRUMENTAL  EVALUATION  AND  MEASUREMENTS 
A.  Haystack  Antenna 

All  of  the  continuum  measurements  and  many  of  the  spectral  line  measurements  described 
in  this  report  were  made  with  the  120-foot-diameter  antenna  of  the  Haystack  Research  Facility 
of  Lincoln  Laboratory.  This  is  a  paraboloidal  antenna  enclosed  in  a  radome  and  located  in 
Tyngsboro,  Massachusetts.  It  is  azimuth-elevation  mounted  and  both  the  surface  accuracy  and 
the  pointing  accuracy  permit  operation  down  to  about  8-millimeter  wavelengths.  This  antenna 
has  been  described  in  some  detail  by  Weiss  (1965)  and  additional  information  can  be  found  in 
Lincoln  Laboratory  Report  "Radio  Physics  and  Astronomy"  1965:1,  Allen  (1967),  and  Stuart 
(1968).  The  antenna  is  pointed  in  real  time  by  computer  programs  written  for  the  Univac  490 
(U490).  'Hiese  programs  are  described  in  a  series  of  Lincoln  Laboratory  Technical  Notes  [see 
Heart,  et  aL  (1965)  and  other  references  therein]. 

Allen  (1967)  measured  the  gain  and  beamwidth  of  the  Haystack  antenna  as  a  function  of  ele¬ 
vation  at  8  and  15  GHz  and  analyzed  the  behavior  of  these  parameters  as  a  function  of  the  tem¬ 
perature  of  the  antenna  structure.  These  measurements,  however,  were  made  before  the  re¬ 
cent  rerigging  [Stuart  (1968)].  In  the  present  study  a  number  of  Discrete  Source  Scan  (DSS)  runs 
(see  Sec.  Ill- A)  were  made  to  determine  the  gain,  beamwidth  and  pointing  accuracy  of  this  an¬ 
tenna  at  various  observing  frequencies  and  with  different  feeds.  The  results  of  this  antenna  cal¬ 
ibration  are  summarized  in  Table  I. 

The  aperture  efficiency  at  1.6  and  1.7  GHz  is  based  on  continuum  measurements  of  Cas  A. 
The  Cas  A  flux  was  calculated  from  a  formula  suggested  by  Allen  (1967),  namely 

S  =  6500(^/0.4)'0-785  fu 

for  v  in  GHz  for  the  epoch  1964.4,  and  a  secular  decrease  of  1.1  percent  per  year.  This  gives, 

—  7  A  7 

for  example,  2055  fu  for  1.66  GHz  for  1967.4.  One  fu  is  10  watts/M  /Hz.  The  aperture  effi¬ 
ciency  at  8  and  15  GHz  is  based  on  continuum  measurements  of  W75-DR21,  which  was  taken  to 
be  a  point  source  and  for  which  a  flux  of  20 fu  at  both  these  frequencies  was  assumed  [see  Ryle 
and  Downs  (1967)  and  Mezger,  et  al .  (1967)]. 

£  is  the  ratio  of  the  source  flux  in  a  given  sense  of  polarization  to  the  observed  antenna 
temperature,  and  is  equal  to  k/A^  where  k  is  Boltzmann's  constant  and  is  the  collecting  area 
of  the  antenna.  However,  since  the  flux  published  for  a  continuum  source  is  usually  the  total 
flux  (Sq),  i.e.,  the  sum  of  the  fluxes  in  two  orthogonal  polarizations,  a  value  for  the  conversion 
factor  from  antenna  temperature  to  flux  equal  to  2|  is  often  used  [e.g.,  by  Allen  (1966)]  to  con¬ 
vert  a  measured  antenna  temperature  in  one  sense  of  polarization  into  a  total  flux  for  a  source 
presumed  to  be  unpolarized. 
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Allen  (1967)  found  some  dependence  of  the  Haystack  antenna  gain  on  elevation  angle  and  on 
differential  temperature  across  the  antenna  structure.  The  gain  variation  with  elevation  was 
relatively  small  (perhaps  10  percent  in  going  from  the  zenith  to  an  elevation  angle  of  20°)  and 
also  the  antenna  system  has  been  changed  considerably  since  Allen's  measurements  were  made. 
The  temperature  differential  across  the  antenna  surface  has  been  substantially  reduced  [Stuart 
(1968)].  For  these  reasons,  these  effects  have  been  neglected. 

The  pointing  of  the  Haystack  antenna  was  calibrated  in  November  1965  and  again  in  November 
1967  (after  rerigging)  by  doing  a  large  number  of  DSS  runs  (see  Sec.  III-A)  on  a  number  of  small- 
diameter  continuum  radio  sources  whose  positions  are  accurately  known.  Each  OSS  produces 
an  azimuth  and  an  elevation  offset  which  equal  the  differences  between  the  command  positions 
and  the  position  at  which  the  program  actually  finds  the  source. 

In  order  to  understand  these  errors  and  predict  them  for  subsequent  observations,  they  were 
fit  in  a  least-square  sense  with  equations  that  take  into  account  azimuth  and  elevation  axis  tilts 
and  skews  as  well  as  the  dominant  sag  or  flexure  term  in  elevation.  These  equations  are 

A  A  =  [ (p  cos  (A  —  K  ^)  +  e]  tan  E  +  D  sin  E  + 


and 


where 


AE  =  —  <p  sin  (A  —  K^)  +  3  cos  E  + 


A  =  azimuth 


E  =  elevation 


<t> 

90°  -  Ka 


e 


D 

d  +  ca 


tilt  of  the  antenna-azimuth  axis  with  respect  to  the  zenith 

azimuth  of  the  antenna-azimuth  axis  (on  the  celestial 
sphere) 

offset  of  the  antenna  elevation  axis  from  the  plane  per¬ 
pendicular  to  the  antenna  azimuth  axis 

offset  of  the  electrical  boresight  axis  from  the  plane 
perpendicular  to  the  antenna  elevation  axis 

offset  of  the  zero  position  of  the  azimuth  encoders  from 
true  North 


3  -  sag  or  flexure  coefficient 


3  +  Cp  =  offset  of  the  zero  position  of  the  elevation  encoders  from 
true  horizontal  (at  A  =  K^). 


These  equations  are  derived  and  the  program  discussed  in  Meeks,  et  al .  (1968).  The  results 
of  the  two  pointing  correction  measurements  are  shown  in  Table  II.  The  1967  program  was  per¬ 
formed  at  15  GHz  while  the  1965  program  was  carried  out  at  8  GHz  and  so  had  a  larger  beamwidth 
and  used  some  radio  sources  whose  angular  structure  is  not  negligible  compared  with  the  beam- 
width.  For  these  reasons  the  1967  program  should  be  more  accurate  in  addition  to  showing 
changes  in  the  antenna  system  since  1965.  The  parameters  determined  by  this  process  were 
used  to  prepare  a  table  which  was  incorporated  into  the  computer  pointing  program  for  subse¬ 
quent  observations. 


B.  Millstone  Antenna 

Some  of  the  spectral  line  measurements  described  in  this  report  were  made  with  the  84-foot- 
diameter  antenna  of  the  Lincoln  Laboratory's  Millstone  Facility  in  Westford,  Massachusetts,  a 
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TABLE  II 

HAYSTACK  POINTING  CALIBRATION 

Parameter 

Nov.  1965  Value 

Nov.  1967  Value 

Units 

4> 

from  azimuth  errors 

4.51 

5.30 

mdeg 

ka 

25 

41 

deg 

e 

-4.  77 

-3.  98 

mdeg 

D 

9.02 

8.55 

mdeg 

CA 

-11.  13 

-25.  58 

mdeg 

9 

from  elevation  errors 

3.01 

4.50 

mdeg 

ka 

27 

27 

deg 

P 

-51. 70 

-53.  38 

mdeg 

CE 

61.62 

57.08 

mdeg 

<(6A  cosE)2)’/2 

3.  1 

2.9 

mdeg 

<(6E)2>1/2 

3.6 

3.6 

mdeg 

half  mile  from  Haystack.  The  Millstone  antenna  is  also  azimuth-elevation  mounted  and  is  usable 
down  to  about  6-centimeter  wavelengths.  The  results  of  the  gain  and  beamwidth  calibration  of 
this  antenna,  using  the  horn  feed  at  the  prime  focus,  are  listed  in  Table  I. 

A  pointing  calibration  program  similar  to  that  done  on  Haystack  was  carried  out  on  the 
Millstone  antenna  in  November  1967.  The  digital  encoders  which  monitor  the  pointing  position 
of  the  Millstone  antenna  are  quantized  in  40  arcsecond  (11  millidegree)  steps.  Therefore,  as  a 
rule  of  thumb,  we  cannot  expect  the  pointing  errors  to  be  reduced  much  below  this  level.  From 
19  DSS  runs  the  elevation  errors  were  found  to  have  an  average  value  of  41  millidegrees  with  a 
scatter  around  this  value  from  +17  to  —15  millidegrees.  There  is  no  obvious  pattern  to  these 
errors  when  plotted  against  either  azimuth  or  elevation,  so  a  40  millidegree  constant  elevation 
bias  was  used  in  the  pointing  program  for  subsequent  observations.  The  azimuth  errors  scaled 
by  cosine  of  elevation  showed  an  average  value  of  7  millidegrees  and  a  scatter  similar  to  the 
elevation  errors.  These  azimuth  errors  were  neglected  in  subsequent  observations. 

C.  Front  Ends 

This  report  contains  measurements  made  at  three  frequencies  —  8  and  15  GHz  continuum  and 
1.6  GHz  spectral  line  —  using  at  least  two  different  front  ends  at  each  frequency.  Measurements 
were  made  both  before  and  after  rerigging  [Stuart  (1968)].  Table  III  contains  a  list  of  the  front 
ends  used  and  some  of  their  important  characteristics. 

The  figure  of  merit  is  the  calculated  rms  noise  for  1 -second  integration.  In  calculating  the 
figure  of  merit,  the  bandwidth  for  the  front  ends  used  for  continuum  measurements  was  taken 
from  the  table;  however,  for  the  front  ends  used  for  spectral  line  work  (1.6  GHz),  the  bandwidth 
was  taken  to  be  1  kHz  to  correspond  to  the  usual  resolution  used.  To  convert  each  figure  of 
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merit  to  brightness  temperatures  and  flux  units,  the  numbers  from  Table  1  for  the  corresponding 
feed  and  antenna  were  used.  For  those  front  ends  used  with  more  than  one  feed,  the  numbers 
from  Table  1  were  averaged. 

The  numbers  in  Table  III  are  only  estimates;  in  particular  the  bandwidth  and  system  tem¬ 
perature  of  the  maser  and  the  paramps  depend  on  the  tuning,  and  the  system  temperature  rises 
at  low  elevations. 

D.  The  Spectral  Line  Receiver 

Figure  2  is  a  simplified  block  diagram  of  the  spectral  line  receiver  in  a  typical  configuration. 
The  dashed  lines  enclose  optional  items  which  were  used  for  some,  but  not  all,  of  the  observa¬ 
tions  described  herein. 

The  coaxial  switch  following  the  antenna  selects  right  or  left  circular  polarization,  or,  with 
the  system  reconfigured,  horizontal  or  vertical  linear  polarization,  when  used  with  the  Andrew 
feed  on  Haystack  or  the  horn  feed  on  Millstone.  With  the  Clavin  feed  on  Haystack,  this  coaxial 
switch  is  not  used  and  the  feed  can  be  rotated  by  hand  to  any  linear  polarization  position  angle. 

The  calibration  noise  tube  inserts  45 °K  additional  noise  into  the  antenna  side  of  the  ferrite 
switch.  This  45°K  was  measured  (by  J.  C.  Carter)  using  a  room  temperature  and  a  liquid- 
nitrogen-cooled  load  substituted  for  the  antenna.  The  parametric  amplifier  itself  provides  image 
rejection;  however,  a  preceding  bandpass  filter  was  sometimes  used  to  prevent  the  Millstone 
radar  transmitter  (at  about  1295  MHz)  from  interfering. 

The  local  oscillator  was  phase  locked  to  a  signal  derived  from  the  station  frequency  standard 
using  a  general  purpose  frequency  synthesizer.  Frequency  switching  is  possible  with  the  Hay¬ 
stack  system,  but  was  not  used  for  any  of  the  observations  reported  here.  Frequency  errors 
from  this  source  are  less  than  1  Hz.  Sometimes  a  28-MHz  crystal  oscillator  was  used  as  well, 
but  its  frequency  was  known  to  about  50  Hz.  The  30- MHz  to  video  frequency  converter  contains 
additional  crystal  oscillators  whose  frequencies  could  be  monitored.  After  some  initial  difficul¬ 
ties  with  this  unit,  we  now  believe  that  the  total  frequency  error  from  hardware  causes  is  less 
than  about  20  Hz.  Some  of  the  early  measurements  were  made  before  this  was  achieved,  and  so 
may  be  somewhat  worse.  Section  1V-B  contains  a  discussion  of  frequency  errors  that  may  arise 
in  the  program  that  calculates  the  observer's  velocity. 

The  video  frequency  converter  selects  a  frequency  window  (typically  40  kHz)  centered  at 
30  MHz  and  converts  it  to  video,  e.g.,  to  0  to  40  kHz.  This  signal  is  clipped,  sampled,  and  fed 
to  the  autocorrelator  at  a  clock  rate  two  and  a  half  times  the  window  width,  e.g.,  100  kHz  for  a 
40-kHz  window.  The  correlator  accumulates  the  autocorrelation  functions  for  100  ms  and  then 
dumps  them  into  the  U490  computer.  Additional  details  about  this  system  can  be  found  in  Weinreb 
(1963),  and  Rogers  (1967a). 

When  the  receiver  is  Dicke  switching,  the  synchronous  and  gated  detectors  provide  informa¬ 
tion  to  the  U490  computer  to  calculate  and  monitor  the  effective  system  temperature  which  is 
needed  to  scale  the  final  spectra.  When  the  receiver  is  not  Dicke  switching  a  total  power  detec¬ 
tor  provides  this  information. 

in.  DATA  ACQUISITION  -  TECHNIQUES  AND  PROCEDURES 
A.  Continuum 

The  Discrete  Source  Scan  (DSS)  program  was  written  for  the  U490  computer  by  Mrs.  P.  P. 
Crowther  to  specifications  by  Dr.  M.  L.  Meeks.  It  has  been  described  elsewhere  [Allen  and 
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Simplified  block  diagram  of  spectral  line  receiver  in  a  typical  configuration. 


Barrett  (1967)  and  Meeks,  et  al.  (1968)]  so  will  only  be  outlined  here.  The  DSS  program  involves 
an  interaction  between  real-time  pointing  of  the  antenna  and  real-time  continuum  data  reduction 
which  take  place  in  the  same  computer.  The  DSS  program  first  commands  the  antenna  to  make 
a  declination  scan  through  the  pointing  position  of  the  source.  The  declination  position  of  the 
peak  temperature  and  the  half-power  width  in  declination  are  determined  by  fitting  a  pair  of 
straight  lines  to  the  sides  of  the  response  curve.  The  program  then  commands  the  antenna  to 
make  a  right  ascension  scan  through  this  peak  declination.  The  right  ascension  scan  is  proc¬ 
essed  in  the  same  way  as  the  declination  scan.  Then,  using  the  peak  position  just  found,  the 
program  commands  an  11  on-off”  in  right  ascension  such  that  the  same  track  in  azimuth-elevation 
coordinates  is  followed  for  the  "off”  as  for  the  "on."  The  results  of  the  DSS,  which  are  typed 
out  immediately  after  the  scan  is  finished,  consist  of  the  offsets  from  the  command  position  to 
the  position  where  the  program  actually  found  the  peak  temperature,  and  the  value  of  this  peak 
temperature  from  the  "on-off."  These  offsets  are  given  in  both  the  azimuth-elevation  and  right 
ascension-declination  coordinate  systems.  The  program  also  types  the  response  widths  in  right 
ascension  and  declination. 

All  the  continuum  maps  in  this  report  were  made  using  the  drift-curve  technique  with  the 
Haystack  antenna.  The  antenna  remains  stationary  while  the  earth's  rotation  carries  the  source 
through  the  antenna  beam.  This  technique  has  several  outstanding  advantages  and  perhaps  only 
one  minor  disadvantage  which  are  summarized  below. 

(1)  The  pointing  errors  are  minimized.  It  is  easier  to  make  the  antenna 
remain  stationary  at  a  specified  position  than  to  track  at  a  specified 
speed. 

(2)  The  contribution  to  the  antenna  temperature  due  to  side  lobes  and  back 
lobes  of  the  antenna  pattern  striking  the  ground,  and  also  the  contri¬ 
bution  due  to  the  earth's  atmosphere,  are  almost  exactly  constant  during 
a  drift  curve  and  so  they  can  be  subtracted  out.  Allen  (1967)  reported 
that  a  component  of  the  background  temperature  contribution  varied 
rapidly  with  azimuth  or  elevation  and  he  attributed  this  effect  to  the 
radome.  Such  a  contribution  has  no  effect  on  a  drift  curve. 

(3)  The  baseline  temperature,  which  corresponds  to  the  "off"  temperature 
in  an  "on-off"  technique,  is  determined  from  the  average  of  many  data 
points.  Therefore,  it  has  a  lower  noise  component  than  the  "off"  of  an 
"on-off"  that  contains  the  same  integration  time  as  the  "on."  In  See.  III-B, 
an  almost  exactly  analogous  technique  for  spectral  line  observing  is  an¬ 
alyzed  and  shown  to  be  superior. 

The  only  real  disadvantage  of  the  drift  curve  technique  has  to  do  with  the  effective  integration 
time  which  is  equal  to  the  time  required  for  the  earth's  rotation  to  carry  a  point  through  half  of 
the  antenna's  half-power  beamwidth.  The  observer  may  get  any  effective  integration  time  for  a 
map  which  is  an  integral  multiple  of  this  time  for  a  single  drift  curve,  but  he  cannot  spend  less 
time  in  making  a  map.  This  restriction  was  no  disadvantage  for  the  maps  in  this  report. 

One  other  point  deserves  mention.  In  making' a  drift  curve  map  one  tacitly  assumes  that 
the  baseline  region  of  the  sky  outside  the  region  to  be  mapped  has  negligible  emission  With 
possibly  one  exception  mentioned  in  Sec.  V-H,  this  condition  was  made  to  hold  for  the  maps  in 
this  report  by  taking  sufficiently  long  drift  curves.  However,  this  point  may  be  a  serious  prob¬ 
lem  with  drift  curve  maps  made  at  lower  frequencies  or  with  more  sensitive  systems. 

B.  Spectral  Line 

The  spectral  line  data  for  this  report  were  taken  with  the  100-channel  one-bit  digital  auto¬ 
correlator  described  in  Sec.  II-D.  Several  modes  of  operation  were  possible  including  Dicke 
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switching  and  frequency  switching  (with  Haystack).  During  the  course  of  the  observations,  the 
nonswitched  mode  of  operation  also  became  possible.  The  autocorrelator  accumulates  autocor¬ 
relation  functions  for  only  100  ms  before  dumping  them  into  the  computer.  Originally,  this  in¬ 
terval  corresponded  to  half  a  cycle  of  Dicke  or  frequency  switching.  However,  it  was  found  that 
the  system  noise  temperature  as  a  function  of  frequency  across  the  window  width  and  referred 
to  the  antenna  (which  is  the  relevant  parameter  with  a  clipped  system),  did  not  change  appreci¬ 
ably  over  a  period  of  many  minutes  or  hours.  It  is  therefore  possible  to  accumulate  "off"  data 
for  many  minutes  with  the  antenna  pointing  off  but  near  the  source,  and  then  move  the  antenna 
onto  the  source  for  "on"  data.  This  procedure,  referred  to  as  the  nonswitched  mode  of  opera¬ 
tion,  should  probably  be  called  "antenna  switching."  It  has  the  following  advantages: 

(1)  The  frequency  dependent  VSWR  of  the  antenna  as  compared  with  the  VSWR 
of  the  dummy  load,  which  may  produce  a  sloping  or  curved  baseline  with 
Dicke  or  frequency  switching,  has  almost  no  effect  on  the  data  taken  in  the 
nonswitched  mode.  It  has  no  effect  at  all  if  the  continuum  temperature 
seen  by  the  antenna  during  the  "off"  is  the  same  as  during  the  "on." 

(2)  The  ferrite  (Dicke)  switch  may  be  completely  removed  from  the  circuit  to 
reduce  the  system  temperature.  (This,  of  course,  is  also  true  of  fre¬ 
quency  switching.) 

(3)  The  "offs"  and  "ons"  need  not  contain  the  same  integration  time  and  one 
"off"  may  be  used  with  several  "ons"  under  certain  circumstances.  The 
rule  for  using  a  single  "off"  with  several  "ons"  is  that  the  resulting  spec¬ 
tra  should  not  be  averaged  or  added  unless  first  shifted  by  one  or  more 
equivalent  filter  widths.  Between  two  spectra  using  the  same  "off,"  the 
receiver  noise  is  partly  correlated  in  the  same  equivalent  filters,  but 
uncorrelated  in  adjacent  filters.  Since  more  than  one  frequency  window 
was  wanted  in  almost  all  the  spectra  taken  for  this  report,  this  technique 
was  quite  valuable  in  reducing  the  integration  time  necessary  for  a  given 
signal-to-noise  ratio. 

Given  that  one  "off"  will  be  used  with  several  "ons,"  the  question  arises  as  to  how  one  should 
divide  up  a  given  block  of  observing  time  among  the  "ons"  and  the  "off"  to  produce  minimum 
noise  on  the  resulting  spectra.  This  problem  is  analyzed  below. 

Let  N  be  the  number  of  "ons"  wanted  (to  be  used  with  a  single  "off"),  t  be  the  time  for  each 

"on,"  ck  be  the  ratio  of  "off"  time  to  t,  and  T  be  the  total  observing  time  which  is  (a  +  n)  t.  Then 

the  noise  for  each  "on"  is  «t  *^2  and  for  the  "off"  is  a(Qft)  1  2  The  total  noise  for  the  "on-off" 

is  ^[(o  +  1  )/of t] ^  2  which  is  (a  +  N  +  1  +  N/a)*  2.  We  differentiate  this  with  respect  to  o,  set 

1  /2 

the  result  to  zero,  and  get  a  =  N  '  .  This  means  that  we  minimize  the  noise  on  each  of  the  re- 

1  /2 

suiting  spectra  by  spending  N  times  as  much  time  on  the  "off"  as  on  each  "on."  For  example, 
we  may  divide  a  5-hour  observing  session  into  a  one-hour  "off"  and  sixteen  15-minute  "ons." 

For  this  example,  the  noise  on  the  resulting  spectra  will  be  reduced  by  a  factor  of  0.584  com¬ 
pared  with  Dicke  or  frequency  switching  for  the  same  total  observing  session. 

A  similar  analysis  may  be  applied  to  the  drift  curve  technique  in  continuum  observations. 

The  baseline  interval  corresponds  to  the  "off"  and  each  X/(2D)  section  of  the  record  corresponds 
to  an  "on." 


C.  Interferometer 

The  Haystack-Millstone  interferometer  has  been  described  in  some  detail  by  Rogers  (1967) 
and  Rogers  (1968)  and  so  will  only  be  outlined  here.  The  antennas  and  front  ends  for  this  inter¬ 
ferometer  are  described  in  Sec.  11.  A  reference  frequency,  from  which  the  local  oscillator  fre¬ 
quency  was  synthesized,  was  sent  on  a  coaxial  cable  from  Haystack  to  Millstone  and  the  Millstone 
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IF  signals  were  sent  back  to  Haystack  on  the  same  cable.  The  electrical  length  of  this  cable  was 
held  constant  by  a  motor-driven  mechanical  line  stretcher  (trombone  section).  This  technique 
virtually  eliminates  instrumental  phase  drift. 

The  most  desirable  way  of  obtaining  the  fringe  amplitude  and  phase  for  a  spectral  line  inter¬ 
ferometer  involves  cross-correlating  the  received  signals.  No  cross-correlator  was  available, 
so  Rogers  developed  an  add-subtract-autocorrelation  technique  which  yields  the  same  informa¬ 
tion  with  a  slightly  poorer  signal-to-noise  ratio.  The  Haystack  and  Millstone  IF  signals  are 
first  added  for  100  ms  and  the  autocorrelator  processes  the  resulting  bit  string  and  dumps  the 
results  into  the  U490  computer.  The  IF  signals  are  then  subtracted  for  100  ms,  etc.  A  program 
in  the  U490  computer  (a  modification  by  Mrs.  P.  P.  Crowther  of  AUTOCOR)  writes  these  autocor- 
relation  numbers  onto  magnetic  tape,  along  with  timing  and  other  information,  for  subsequent 
processing  by  a  post-real-time  program  in  the  CDC  3300  computer.  This  post-real-time  pro¬ 
gram  (MILSTAK  or  MILSTAKS,  written  by  J.  M.  Moran,  Jr.  and  A.  E.  E.  Rogers)  produces  the 
fringe  amplitude  and  phase  as  a  function  of  frequency  across  the  spectral  window  for  each  inte¬ 
gration  period  which  was  typically  10  to  30  minutes. 

The  original  observations  with  this  interferometer  by  Rogers,  e_t  al .  (1966)  were  of  unre¬ 
solved  OH  emission  sources.  Positions  and  upper  limits  for  the  angular  sizes  of  such  sources 
were  obtained.  Observations  of  continuum  sources  (as  calibrators)  were  also  made.  Rogers 
shows  that  for  observations  of  continuum  sources  for  which  the  total  delay  is  sufficient  to  pro¬ 
duce  many  interferometer  fringes  across  the  spectral  window,  the  correct  continuum  fringe  am¬ 
plitude  and  phase  may  be  obtained  from  a  weighted  vector  average  over  the  spectral  window.  If, 
however,  there  is  insufficient  total  delay  relative  to  the  spectral  window,  then  a  clipped  inter¬ 
ferometer  cannot  give  the  correct  answer.  This  is  easy  to  see  when  there  is  no  delay.  For  this 
case,  only  the  power  level  input  to  the  clipper  changes  with  phase  switching  and  with  time.  Since 
the  clipper  is  insensitive  to  the  power  level  at  its  input,  no  information  about  this  power  level  is 
available  to  the  computer  unless  it  is  recorded  separately. 

This  problem  can  be  looked  at  as  one  of  recovering  the  true  autocorrelation  numbers  from 
the  clipped  ones.  There  are  two  parts  to  this  clipping  correction:  (1)  taking  the  trigonometric 
sine  of  the  clipped  autocorrelation  functions  and  (2)  Munnormalizing"  by  multiplying  by  the  equiv¬ 
alent  interferometer  system  temperature.  When  there  are  not  many  interferometer  fringes 
across  the  frequency  window,  this  normalization  factor  varies  on  a  second-by-second  basis 
(i.e.,  it  has  time  fringes),  rather  than  remaining  nearly  constant  through  a  run. 

These  variations  of  the  equivalent  system  temperature  will  be  referred  to  as  continuum 
fringes,  and  the  process  of  recording  the  continuum  fringes  and  using  the  correct  variable  nor¬ 
malization  factor  as  continuum  fringe  compensation.  It  is  necessary  to  do  this  in  order  to  ob¬ 
serve  a  continuum  source  with  a  narrow  spectral  window,  even  if  the  continuum  information  it¬ 
self  is  not  wanted  The  only  other  alternative  would  be  to  introduce  additional  delay  into  the  long 
leg  of  the  interferometer,  to  produce  many  fringes  across  the  spectral  window. 

To  implement  the  continuum  fringe  compensation  scheme  on  the  Haystack-Millstone  inter¬ 
ferometer,  two  gated  detectors  were  used  to  record  the  quantities: 


and 


D1  =  (eH  +  eM> 


D2  (eH-eM> 


(1) 


(2) 
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where  e^  and  e^  are  the  IF  voltages  from  the  Haystack  and  Millstone  receivers,  respectively. 
The  bar  indicates  a  time  average  long  compared  with  a  cycle  of  the  IF  but  short  compared  with 
a  cycle  of  the  interferometer  fringes;  three  seconds  were  used  in  practice.  The  quantities  D ^ 
and  are  proportional  to  the  factor  which  should  be  used  to  scale  those  clipped  autocorrelation 
functions  that  are  taken  during  the  "add”  and  the  "subtract1’  portions  of  the  switch  cycle,  respec¬ 
tively.  The  quantity  is  proportional  to  the  constant  component  and  is  propor¬ 

tional  to  the  variable  component  of  the  equivalent  system  temperature.  In  practice  the  variable 
component  was  fit  with  a  sinusoid  in  order  to  extrapolate  ahead  in  time,  since  the  continuum 
data  were  delayed  on  the  tape.  This  process  also  reduces  the  noise  contribution. 

The  "sine"  portion  of  the  clipping  correction  i&  somewhat  more  of  a  problem,  because  to 
actually  calculate  the  trigonometric  sine  in  the  computer  for  each  of  the  autocorrelation  numbers 
would  require  an  excessive  amount  of  computer  time  (one  thousand  sines  would  be  needed  for 
each  second  of  data).  However,  the  normalized  autocorrelation  numbers  are  always  small  and 
it  is  a  good  approximation  to  take  the  argument  for  the  sine  (i.e.,  just  multiply  by  tt/2).  In  the 
absence  of  continuum  fringes,  only  the  difference  between  the  two  autocorrelation  numbers  is 
needed  [see  Rogers  (1967b),  p.  14];  in  the  presence  of  continuum  fringes  both  the  sum  and  dif¬ 
ference  are  needed. 

It  was  found  in  obtaining  the  data  in  this  way  that  a  residual  slope  (i.e,,  a  variation  of  fringe 
phase  with  frequency  across  the  spectral  window)  was  sometimes  seen  (see  See.  \  -J)  and  this 
instrumental  effect  may  possibly  be  due  to  the  omission  of  the  sine  clipping  correction. 

IV.  DATA  REDUCTION  -  TECHNIQUES  AND  PROCEDURES 

A.  Continuum 

1.  Continuum  Maps 

The  real-time  program  [Stylos  (1965)]  in  the  U490  provides  a  preliminary  reduction  of  the 
incoming  raw  data  by  converting  the  voltages  to  temperatures  which  are  printed  and  recorded 
on  magnetic  tape.  A  map  reduction  program  written  for  the  CDC  3300  computer  by  N.M.  Brenner 
to  specifications  by  Dr.  M.  L.  Meeks  reads  these  tape  records  and  processes  the  data  to  eventu¬ 
ally  obtain  a  map.  The  details  of  this  program  are  given  by  Brenner  [N.  M.  Brenner  (1969),  in 
preparation]  and  arc  outlined  here. 

The  program  is  divided  into  three  parts,  lettered  A,  B,  and  C.  Program  part  A  reads  the 
data  tapes  made  by  the  U490  and  locates  the  drift  curves  on  them.  The  tape  record  contains  the 
time,  right  ascension,  declination,  azimuth,  and  elevation  as  well  as  the  temperatures,  so  the 
program  checks  that  the  azimuth,  elevation  and  declination  are  constant,  and  the  right  ascension 
is  drifting  at  the  earth’s  rotation  rate  (within  prescribed  tolerances).  These  criteria  define  a 
drift  curve  When  the  complete  drift  curve  is  located,  program  part  A  subtracts  off  a  baseline 
temperature  and  filters  the  data  by  convolving  the  function  sin  (Ax)/ (Ax),  whose  Fourier  trans¬ 
form  is  unity  at  all  frequencies  up  to  the  cut-off  frequency  of  the  antenna  and  zero  for  all  higher 
frequencies  [see  Bracewell  and  Roberts  (1954)  and  Bracewcll  (1965)).  This  convolution  filters 
out  receiver  noise  components  that  cannot  represent  signals  because  they  arc  at  a  frequency  that 
cannot  be  accepted  by  the  antenna.  And  this  convolution  also  provides  an  exact  method  of  inter¬ 
polation  to  get  temperatures  at  intermediate  positions  between  the  original  data  points.  We  neg¬ 
lect  the  effect  of  the  true  integrator  and  sampler  which  filter  the  data  ahead  of  the  computer. 

This  is  a  good  approximation,  provided  that  the  integration  time  is  small  compared  with  a  half 
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cycle  of  the  antenna  cut-off  frequency.  For  a  drift  curve  with  Haystack  the  cut-off  frequency  is 
0.1  3  Hz  at  15  GHz  and  0.  07  Hz  at  8  GHz  In  practice,  1  second  of  integration  for  15  GHz  obser¬ 
vations  and  2  seconds  of  integration  for  8  GHz  observations  were  used. 

Program  part  A  plots  the  original  data  points  and  the  smoothed  curve  fit  through  them  on 
the  cathode  ray  tube  display  (CDC-280)  and  the  operator  may  accept  or  reject  a  scan,  or  modify 
it  before  accepting  it.  Accepted  scans  are  recorded  on  magnetic  tape  for  further  processing  by 
program  part  B. 

Program  part  B  reads  these  processed  scans  output  from  A  and  aligns  and  averages  scans 
taken  at  the  same  declination.  Again,  they  are  plotted  on  the  cathode  ray  tube  display  and  the 
operator  may  accept  or  reject  scans.  In  program  part  B  a  correction  for  atmospheric  atten¬ 
uation  and  precession  corrections  back  to  1950  may  be  applied  as  the  operator  directs.  The 
correction  for  atmospheric  attenuation  was  based  on  ground-level  values  for  temperature,  pres¬ 
sure,  and  humidity,  and  was  calculated  from  formulas  suggested  by  Allen  (1967).  The  correction 
for  precession  is  computed  for  the  center  of  each  scan.  For  very  long  scans  in  certain  parts  of 
the  sky  there  may  be  some  appreciable  error  in  the  precession  correction  for  points  on  the  ends 
of  the  scans.  This  error  is  a  small  part  of  an  arcminute  on  the  maps  in  this  report.  The  aver¬ 
aged  scans  are  recorded  on  magnetic  tape  for  further  processing  by  program  part  C. 

Program  part  C  reads  the  averaged  scans  output  from  B  and  further  processes  them  to  ob¬ 
tain  a  two-dimensional  map.  The  tape  from  B  contains  only  one  averaged  scan  for  each  step  in 
declination  and  these  are  placed  in  a  two-dimensional  matrix  of  numbers  which  is  filtered  by 
being  convolved  by  J^(Ax)/(Ax).  This  function  is  the  two-dimensional  analogue  of  the  filter  func¬ 
tion  used  in  part  A.  Again  the  effect  is  to  filter  out  only  receiver  noise  and  to  provide  inter¬ 
polated  temperatures. 

The  filtered  temperatures  represent  the  best  estimate  of  the  antenna  temperature  over  the 
mapped  region  that  is  possible  with  the  given  data.  Restoration  (see  Bracewell  (1958)]  could  be 
performed  at  this  time  by  choosing  a  different  filter  function  instead  of  J^(Ax)/(Ax);  however, 
this  was  not  done  for  any  of  the  maps  presented  in  this  report.  If  the  data  are  spaced  in  decli¬ 
nation  in  steps  larger  than  \/(2D),  where  \  is  the  wavelength  of  observation  and  D  is  the  diam¬ 
eter  of  the  antenna,  then  one  must  increase  A  to  correspond  to  a  smaller  antenna  and  larger 
beamwidth.  If  the  data  are  spaced  in  declination  in  nonuniform  steps,  then  there  is  no  way  to 
process  the  data  with  this  program.  The  critical  step  interval  X/(2D)  for  Haystack  is  0°016  at 
15  GHz,  and  0*029  at  8  GHz. 

It  is  also  possible  to  increase  A  in  the  filter  formula  to  smooth  a  map  to  a  larger  beamwidth 
The  equivalent  antenna  pattern  for  this  modified  map  will  not  be  the  same  as  the  antenna  pattern 
of  the  same  antenna  used  at  a  longer  wavelength;  to  produce  such  an  equivalent  antenna  pattern 
would  require  a  more  complicated  filter  function.  However,  such  a  modified  map  is  useful  in 
showing  weak  extended  features  which  otherwise  might  be  hidden  in  noise.  A  rough  comparison 
with  longer  wavelength  maps  is  also  possible. 

The  filtered  temperatures  may  be  presented  in  many  forms  as  a  matrix  of  numbers;  as  a 
shaded  map,  as  a  ruled  surface  map;  or  as  a  contour  diagram,  which  is  the  form  used  for  all 
the  maps  in  this  report.  These  maps  are  presented  on  the  cathode  ray  tube  display  and  simul¬ 
taneously  on  another  cathode  ray  tube  to  be  photographed. 

2.  Continuum  Fluxes  and  Other  Source  Parameters 

There  are  at  least  two  possible  ways  to  obtain  the  flux  of  an  individual  source  from  such  a  map. 
One  may  fit  the  observed  distribution  with  a  source  model  function  such  as  a  two-dimensional 
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Gaussian  and  take  the  volume  under  the  model  to  represent  tne  integrated  antenna  temperature. 
For  a  two-dimensional  Gaussian  model  the  integrated  antenna  temperature  is  just 

J  Ta(0)  do  =  (1.133)  T  e6eo  (3) 

source  r 

where  T  is  the  peak  temperature  of  the  Gaussian  and  ©^  and  ©^  are  its  two  widths  to  half¬ 
power  in  declination  and  right  ascension. 

Or,  one  may  add  up  the  actual  temperature  numbers  in  the  two-dimensional  matrix  that  cor¬ 
responds  to  the  map  of  the  source.  This  summation  may  be  performed  by  program  part  C. 
Bracewell  and  Roberts  (1954)  have  shown  that  the  integral  under  the  surface  defined  by  the  an¬ 
tenna  temperature  map  is  given  exactly  by  the  sum  of  the  numbers  in  such  a  matrix;  no  more 
complicated  integration  formula  is  needed. 

The  sources  observed  for  this  study  fall  into  two  classes:  (1)  Those  with  an  angular  struc¬ 
ture  that  is  small  or  comparable  to  the  antenna  beam.  For  these  sources  the  Gaussian  model 
is  most  appropriate  and  this  method  yields  lower  error  estimates  for  the  flux.  (2)  Those  sources 
with  an  angular  structure  that  is  large  compared  with  the  antenna  beam.  None  of  these  sources 
is  adequately  approximated  by  a  two-dimensional  Gaussian  model  or  any  other  simple  model,  so 
the  summation  method  was  used  to  calculate  the  flux.  The  two  methods  were  compared  on  the 
source  W75-DR21  (see  Table  I).  For  a  complex  source,  the  integrated  antenna  temperature  as 
determined  from  a  Gaussian  model  is  almost  always  smaller  than  the  corresponding  value  deter¬ 
mined  by  the  summation  method. 

When  the  integral  of  T^  is  known,  one  may  use  the  well-known  formula 


_2k  i  r 

■V  2  Tj  n  J 


source 


Ta(Q)  dQ 


(4) 


to  get  the  flux  density  (S^)  of  a  source  if  the  beam  efficiency  (i?^)  is  known,  or  one  may  use  this 
formula  to  get  the  beam  efficiency  as  measured  on  a  source  whose  flux  is  known.  In  this  for¬ 
mula  k  is  Boltzmann's  constant,  X  is  the  wavelength  of  observation  and  d H  is  an  element  of  solid 
angle  on  the  source.  The  constant  2k/(X2r^)  is  the  (mapping  efficiency)”1  listed  in  Table  I.  In 
this  formula,  the  l/7?^  factor  converts  the  integrated  antenna  temperature  to  what  would  be  ob¬ 
served  using  a  lossless  antenna  that  accepted  radiation  only  in  its  main  beam,  the  2k/X^  factor 
converts  a  brightness  temperature  to  a  specific  intensity  according  to  the  Rayleigh-Jeans  approx¬ 
imation  to  the  Planck  radiation  law,  and  the  integral  of  the  specific  intensity  over  the  solid  angle 
of  the  source  gives  the  flux  density.  The  beam  efficiency  r\ ^  is  defined  by 


,f  .  .  G(e.*)g 

main  beam 


(5) 


where  G  is  the  antenna  power  gain.  The  region  over  which  this  integral  is  to  be  taken  —  the 
"main  beam"  —  is  somewhat  arbitrary.  For  most  antennas,  this  integral  will  change  only  very 
slowly  as  additional  solid  angle  is  included  in  the  integral  after  the  main  lobe  and  perhaps  the 
first  sidelobes  of  the  power  pattern  are  included.  Equation  (4)  is  only  strictly  correct  if  this 
variation  is  negligible. 

Since  Eq.  (4)  is  identical  with  the  corresponding  formula  for  source  flux  in  terms  of  bright¬ 
ness  temperature  except  for  the  1/17^  factor,  we  may  define  an  apparent  brightness  temperature 
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by  T^(f2)/rj^.  This  apparent  brightness  temperature  will  be  equal  to  the  true  brightness  tem¬ 
perature  for  any  source  whose  angular  structure  is  not  smeared  by  the  antenna  beam,  i.e.,  for 
a  source  that  is  completely  resolved.  It  is  thus  analogous  to  the  apparent  flux  density  in  the 
spectral  plots  given  in  this  report.  The  apparent  flux  density  will  be  equal  to  the  true  flux  den¬ 
sity  for  a  source  whose  spectral  features  are  completely  resolved  by  the  spectral  resolution  used. 

Equation  (4)  is  applicable  to  either  a  point  or  an  extended  source;  however,  for  a  source 
whose  angular  structure  is  small  compared  with  the  antenna  beam,  an  alternative  formula  is 
usually  more  useful,  namely: 


S 


(6) 


18-3-11893 


where  is  the  effective  collecting  area  of  the  antenna  and  T^  is  the  antenna  temperature  seen 

in  one  sense  of  polarization  on  the  source.  The  constant  k/A^  is  £  listed  in  Table  I.  The  factor  2 

in  both  these  formulas  is  necessary  whenever  T^ 

refers  to  the  antenna  temperature  measured  in  one 

sense  of  polarization  and  S  refers  to  the  total  flux  of 

r  V  a  ■  O,  S  -  CONSTANTv 

a  source  assumed  to  be  unpolarized  For  all  the  spec¬ 
tral  plots  in  this  report,  the  form  of  Eq.  (6)  was  used 
without  the  2  to  convert  measured  antenna  tempera¬ 
ture  to  the  corresponding  Stokes  parameters. 

H  II  regions  emit  thermal  continuum  by  free-free 
transitions  A  typical  spectrum  from  such  a  region 
is  shown  in  Fig.  3  In  the  low  frequency  part  of  the 
spectrum  the  opacity  is  large  and  the  radiation  cor¬ 
responds  to  an  opaque  black  body.  The  brightness 

temperature  is  constant  over  this  region  of  the  spectrum  and  just  equal  to  T  ,  the  electron  ki¬ 
netic  temperature.  As  the  frequency  increases,  the  opacity  drops  according  to  the  relationship 
[Oster  (1961 )] 


va  *  2,  T  ■  T  *  CONSTANT 

B  • 


Fig.  3.  Typical  radio  continuum  spectrum 
of  an  H  II  region. 


Ne' 


log 


,3/2 


const. 


(7) 


We  have  assumed  a  uniform  medium  inside  the  H  II  region,  i.e.,  constant  T^  and  N^,  the  elec¬ 
tron  density.  For  frequencies  high  enough  so  that  «  1  we  have 


S  »  T  r  S7  —  constant  (8) 

v  2  e  v  s 
c 

if  we  neglect  the  v  dependence  in  the  logarithm.  is  the  solid  angle  subtended  by  the  source 

and  i  is  its  thickness  along  the  line  of  sight.  The  —  log  ^  term  means  that  the  flux  will  droop  a 
bit  from  the  constant  line.  The  flux  in  this  high  frequency  portion  of  the  spectrum  is  approxi¬ 
mately  proportional  to  figN^//T^  Therefore,  if  the  source  is  resolved  and  if  is  known 
from  low  frequency  measurements  or  from  excited  hydrogen  emission  line  measurement,  it  is 
possible  to  use  the  high  frequency  measurements  to  calculate  the  emission  measure  (N^/).  If, 
furthermore,  the  distance  to  the  source  is  known  and  it  is  assumed  that  l  is  the  same  as  the 
transverse  dimensions,  we  may  estimate  N  . 
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2 

The  knee  of  the  curve  in  Fig.  3  corresponds  to  t  al,  and  this  occurs  at  EM/i'J'  **  2  5  [EM 

3  2  " 

in  (el/cm  )  pc  (parsecs);  v  in  MHz]  for  an  assumed  T  of  10,000°K  [see  Kraus  (1966)].  A  typ- 

6  3  2  ^ 

ical  value  for  EM  might  be  10  (el/cm  )  pc  and  this  corresponds  to  =  630  MHz.  Turnover 
frequencies  higher  than  this  (corresponding  to  higher  EM)  have  sometimes  been  observed  [Ryle 
and  Downes  (1967)]. 

In  principle  then,  two  widely  spaced  points  on  the  spectrum  of  such  a  source  together  with 
resolution  in  angle  suffice  for  a  calculation  of  the  electron  temperature,  the  emission  measure, 
and  perhaps  the  electron  density.  In  practice  the  uncertainties  will  be  large  unless  points  below 
the  knee  of  the  curve  can  be  measured. 

Mezger  and  Henderson  (1967)  have  shown  how  to  carry  out  this  analysis  in  detail  for  several 
assumed  source  models.  This  has  not  been  done  for  this  study  for  the  following  reasons:  (1)  The 
8-GHz  measurements  are  almost  never  low  enough  in  frequency  to  define  the  knee  of  the  curve  or 
to  get  Tg;  (2)  As  previously  noted,  those  sources  that  are  resolved  in  angle  are  complex  and 
cannot  be  adequately  described  by  any  of  the  models  suggested  by  Mezger  and  Henderson  (1967); 
(3)  For  those  sources  to  which  this  analysis  might  be  applicable  (such  as  W3  and  W49)  because 
of  incipient  but  not  complete  angular  resolution,  the  error  estimates  for  the  electron  densities 
and  masses  are  so  large  as  to  make  the  measurements  meaningless.  For  a  discussion  of  such 
problems,  see  Aikman  (1968). 

Table  VI  (page  33)  contains  fluxes  determined  as  described  above,  angular  sizes  uncorrected 
for  antenna  beam  smearing,  and  apparent  spectral  indices,  i.e.,  straight  line  slopes  between  the 
two  measured  points. 

B.  Spectral  Line 

1  Real  Time  (AUTOCOR) 

The  real-time  spectral  line  data  reduction  program  (AUTOCOR)  was  written  by  G.  H. 
Conant,  Jr.,  to  specifications  by  Dr.  Sander  Weinreb.  It  has  been  described  in  some  detail  by 
Conant  and  Meeks  (1968)  and  so  will  only  be  outlined  here.  The  description  below  refers  to  the 
latest  version  of  this  program  which  has  been  modified  by  Mrs.  P.  P.  Crowther 

The  AUTOCOR  program  operates  in  the  U490  computer  in  the  control  environment  of  the 
antenna  pointing  system  [Heart,  et  al.  (1965)].  AUTOCOR  accepts  the  100  autocorrelation  num¬ 
bers  from  the  digital  autocorrelator  every  100  ms  and  sums  them  into  either  of  two  arrays  which 
represent  "on"  and  "off"  source  data.  It  may  switch  back  and  forth  between  these  two  arrays 
every  100  ms  for  Dicke  or  frequency  switching  or  it  may  sum  exclusively  into  one  array  for  an 
entire  run  in  the  unswitched  mode.  AUTOCOR  also  accepts  low  speed  data  from  the  synchronous 
and  signal  power  detectors.  During  a  calibration  sequence  this  information  allows  AUTOCOR 
to  compute  the  system  noise  temperature,  and  while  running,  to  monitor  changes  in  antenna 
temperature.  The  length  of  a  run  may  be  selected  in  half-minute  steps  from  \  to  500  minutes. 

At  the  end  of  this  time  AUTOCOR  performs  the  Fourier  transform  necessary  to  get  from  the 
accumulated  autocorrelation  numbers  to  the  powrer  spectrum  (antenna  temperature  vs  frequency). 
The  total  system  temperature  is  used  as  a  scale  factor  for  this  transform. 

Four  temperature  points  are  calculated  per  frequency  resolution  interval  (i.e.,  one  point 
per  250  Hz  for  the  40-kHz  window  used  for  most  of  the  spectra  given  in  this  report).  Although 
two  points  per  resolution  interval  suffice  to  specify  the  function  exactly,  experience  has  shown 
that  it  is  very  difficult  to  draw  in  the  proper  curve  by  hand  if  only  two  points  per  resolution 
interval  are  plotted.  And  straight  lines  connecting  the  points  are  even  less  accurate. 
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AUTOCOR  prints  the  autocorrelation  and  temperature  numbers  on  the  high-speed  printer 
and  makes  a  printer  plot  of  the  temperatures.  The  speetrum  is  thus  available  to  the  observer 
within  a  few  seeonds  after  the  end  of  each  run.  A  magnetic  tape  record  is  also  written  by 
AUTOCOR  at  the  end  of  eaeh  run.  This  record  contains  the  accumulated  autocorrelation  num¬ 
bers,  the  temperature  numbers,  the  system  temperature,  the  date  and  the  times  of  the  start  and 
end  of  the  run,  and  numerous  other  parameters  pertaining  to  the  run. 

2.  Post  Real  Time  (SLRKDUCE) 

The  post-real-time  spectral  line  processing  program  (SLREDUCE),  written  for  the  CDC  3300, 
reads  the  magnetic  tape  records  made  by  AUTOCOR,  and  further  processes  them  to  produce  av¬ 
eraged  spectra,  Stokes  parameters,  and  other  polarization  parameters.  SLREDUCE  has  been 
described  by  Ball  (1968)  and  will  only  be  outlined  here. 

The  instructions  to  SLREDUCE  are  punched  on  cards  and  tell  the  program  which  spectra  to 
search  for  on  the  data  tape  and  what  to  do  with  them.  For  circular  polarization  runs  or  runs  on 
an  unpolarized  source,  SL REDUCE  may  be  directed  to  just  average  the  runs  with  a  frequency 
shift  if  necessary.  If  desired,  this  frequency  shift  may  be  recalculated  by  the  program  using  the 
local  oscillator  synthesizer  setting  and  a  velocity  calculation  subroutine.  This  procedure  was 
followed  for  all  the  data  presented  in  this  report.  Errors  in  the  pre-observation  velocity  calcu¬ 
lation  program  are  therefore  unimportant.  The  largest  error  in  this  velocity  calculation  sub¬ 
routine  is  due  to  the  motion  of  the  earth  with  respeet  to  the  earth-moon  barycenter,  which  is 
neglected,  and  which  amounts  to  about  0.015  km/scc  at  most.  In  order  to  align  runs  exactly,  it 
may  be  necessary  to  shift  some  of  them  by  a  fraction  of  the  separation  of  points  on  the  frequency 
axis,  i.e.,  to  calculate  interpolated  temperatures.  This  may  be  done  either  by  convolving  the 
temperatures  by  sin  (Ay  )/(Ac)  or  by  redoing  the  Fourier  transforms  with  a  shift.  The  latter 
method  was  used  for  aligning  individual  runs,  but  the  convolution  method  was  used  to  align  av¬ 
eraged  runs  for  restoration  (see  See.  IV-R-4). 

3  Stokes  Parameters 

The  Stokes  parameter  definitions  used  in  this  report  are  those  of  Born  and  Wolf  (1964)  ex¬ 
cept  that  the  electrical-engineering  definition  of  the  sense  of  circular  polarization  is  used  rather 
than  the  optical  definition.  This  section  contains  abbreviated  definitions  and  some  relationships 
used  in  the  data  reduction. 

The  four  Stokes  parameters  are  defined  as  follows: 

So  =  TR  +  T[  =  T(0)  +  T(G  +  9°0) 

Sj  -  T(O)  -  T ( 9 0 ° ) 

S2  =  T(45°)  -  T(1  35 0  ) 


where  Tj  and  Tj  are  the  antenna  temperatures  measured  at  right  and  left  circular  polarization, 
respectively,  and  where  T(0)  is  the  antenna  temperature  measured  in  linear  polarization  with  a 
polarization  position  angle  0.  For  a  point  source,  these  Stokes  parameters  defined  in  terms 
of  antenna  temperature  differ  from  the  Stokes  parameters  defined  in  terms  of  flux  by  a  constant 
factor  (see  Sec.  II).  The  polarization  position  angle  for  linear  polarization  is  the  angle  from  the 
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line  of  constant  right  ascension  through  the  source  around  to  the  line  corresponding  to  the  E 
vector  of  the  antenna  projected  onto  the  sky.  Counterclockwise  rotation  is  taken  as  positive 
Values  of  9  differing  by  180°  are  equivalent  except  possibly  for  instrumental  effects.  Right  cir¬ 
cular  polarization  denotes  clockwise  rotation  of  the  E  vector  as  seen  along  the  direction  of  prop¬ 
agation,  i.e.,  the  E  vector  traces  a  right-handed  screw  thread. 

The  above  equations  are  used  to  process  circular  polarization  data  to  determine  and  for 
one  of  two  statistically  independent  determinations  of  SQ.  To  process  the  linear  polarization 
data,  the  relation 

T(9 )  =  Ta  +  Tr  cos  2(0  -  0) 

was  used  The  linear  data  consist  of  a  number  of  values  of  antenna  temperature  (T^)  each  with 
an  associated  polarization  position  angle  (9.).  These  were  fit  to  the  above  equation  in  a  least 
square  sense,  i.e.,  the  quantity 


£[Ti“TA_TB  c°s2(9.  -  0)]2 

was  minimized  for  variations  of  T^,  T^  and  0,  to  give  a  best  estimate  of  these  three  quantities 
at  each  frequency  point.  Then  the  Stokes  parameters  were  calculated  from 

Sd  =  2T^  sin  20 
S2  =  2T g  cos  20 


This  value  of  Sq  from  the  linear  polarization  data  was  compared  with  the  value  of  from  the 
circular  polarization  data  as  a  consistency  check  and  as  a  check  on  the  relative  gains  of  the  feeds 
used  for  the  different  polarizations  (see  Table  I) 

Other  useful  quantities  are  the  fraction  linear  polarization 


Pi  = 


J' 


s2  +  s| 


2T 


B 


the  fraction  circular  polarization 


pc  =  S' 


T  —  T 
R  L 

T  +  T 
R  L 


and  the  total  fraction  polarization 


i 


2  x  2 

Pr-  +  Pi 


R 


2  2  2 
+  S  2  +  S  2 


Note  that  with  this  definition  a  positive  pc  corresponds  to  dominantly  right  circular  polarization 
and  a  negative  pc  to  left.  The  ellipticity  angle  x  is  obtained  from 


sin  2x  = 


R 


2  *  s25  .  S33 
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Fig.  4.  Three-dimensional  Poincare  sphere:  (a)  Tap  view;  (b)  side  view  (section  A-A  above). 

Figure  4,  which  is  a  sketch  of  the  three-dimensional  Poincare  sphere,  is  useful  in  visualizing 
these  relations.  The  fraction  polarization  is  represented  by  the  fraction  of  the  radius  of  the  cir¬ 
cle  inside  P,  the  fraction  linear  polarization  is  represented  by  the  fraction  of  the  radius  inside 
the  projection  of  P  onto  the  equatorial  plane,  and  the  fraction  circular  polarization  is  repre¬ 
sented  by  the  fraction  of  the  radius  inside  the  projection  of  P  along  the  vertical  axis. 

The  angle  20  is  the  longitude  and  the  angle  2\  is  the  latitude  on  the  Poincare  sphere.  Fur¬ 
ther  details  on  the  Stokes  parameters  may  be  found  in  Born  and  Wolf  (1965)  and  Ball  (1  968)  gives 
details  on  how  the  Stokes  parameters  and  other  polarization  parameters  were  obtained  from  the 
data. 


4.  Restoration 

Two  of  the  OH  emission  sources  studied  for  this  report  (NGC  6334  and  W75)  are  actually 
pairs  of  small  diameter  OH  emission  sources  far  enough  apart  to  be  distinguished  with  the  22  arc- 
minute  beam  of  Haystack.  When  one  of  the  two  OH  emission  sources  in  such  a  region  is  being 
observed,  the  other  emission  source  is  on  the  side  of  the  beam  so  that  features  from  the  off- 
beam  position  are  present  on  the  spectrum  but  with  reduced  amplitude.  Using  spectra  taken  with 
the  antenna  pointed  at  each  of  the  emission  points,  and  using  a  knowledge  of  the  antenna  beam 
pattern  and  the  positions  of  the  sources,  it  is  possible  to  calculate  the  spectra  which  would  have 
been  seen  if  the  sources  could  have  been  observed  separately.  This  process  will  be  referred  to 
as  "restoration,”  and  the  resulting  spectra  as  "restored"  spectra.  The  equations  used  in  resto¬ 
ration  are  discussed  below. 

R  R 

Let  T^  (y)  and  T^  (y)  be  the  restored  antenna  temperatures,  i.e.,  the  antenna  temperatures 
(in  a  given  sense  of  polarization)  that  would  be  measured  with  the  other  nearby  source  turned  off. 
And  let  T^(^)  and  T^(^)  be  the  actual  antenna  temperatures  measured.  Then  we  may  write 

T1  (i^)  =  T^V)  +  aT2V) 

T2(v)  =  T*V)  +  oT(V) 

where  a  is  the  fractional  reduction  in  the  antenna  power  pattern  at  the  position  of  the  other 
nearby  source,  a  will  not  be  a  function  of  v  for  the  measurements  at  a  given  OH  transition. 
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We  have  assumed  that  the  antenna  beam  is  symmetrical  and  that  the  antenna  was  pointing  exactly 
at  each  source  in  turn;  if  either  of  these  conditions  fails  to  hold,  then  somewhat  more  compli¬ 
cated  equations  result.  We  also  assume  that  the  polarization  properties  of  the  antenna  beam  are 
the  same  at  the  off-axis  and  on-axis  positions.  We  invert  the  above  equations  to  obtain 


Tfoo  = 


Tj  (y)  -  «T2W 


1  —  a 


and 


T?(v)-aT  Av) 

T2  W  =  - 2^ 

1  —  a 


However,  these  equations  do  not  tell  the  whole  story.  Between  tw'o  sources  even  this  close 
together  there  is  a  noticeable  differential  Doppler  shift,  i.e.,  a  difference  in  the  components  of 
the  observer's  motion  projected  onto  the  lines  of  sight  to  the  sources.  Since  the  differential 
Doppler  varies  as  a  function  of  the  time  of  year  (and  even  varies  through  a  day,  although  the 
daily  variation  is  negligible  for  the  measurements  described  herein)  it  is  necessary  to  restore 
spectra  taken  months  apart  before  they  can  strictly  be  intercompared.  Failure  to  do  this  may 
lead  to  apparent  time  variations  wrhen  in  fact  the  source  itself  is  unchanged.  This  effect  also 
makes  it  necessary  to  calculate  for  each  source  separately  the  conversion  from  frequency  to 
velocity  for  the  abscissa  axis. 

Successful  restoration  depends  on  knowing  the  antenna  pattern  and  the  source  positions  ac¬ 
curately  and  controlling  the  pointing  to  within  a  small  fraction  of  a  beamwidth.  Wc  have  also 
tacitly  assumed  that  the  sources  remain  constant  during  the  time  interval  between  the  measure¬ 
ments  on  the  two  sources.  No  noticeable  difficulties  were  experienced  with  any  of  these  require¬ 
ments  except  with  respect  to  the  pointing  of  the  Millstone  antenna.  This  problem  is  discussed 
under  the  individual  sources  in  Sec.  V. 


C.  Interferometer  Data  Reduction 

1.  Introduction 

This  section  describes  a  procedure  which  may  be  used  to  process  interferometer  data  taken 
on  an  extended  source.  The  data  obtained  from  a  two-antenna  fixed-baseline  interferometer  may 
be  presented  in  the  form  of  a  map,  such  as  w'ould  be  obtained  from  a  single  antenna  with  a  par¬ 
ticular  pattern.  The  pattern  of  this  equivalent  single  antenna  depends  on  the  source  declination 
as  well  as  the  baseline  parameters  and  for  noncircumpolar  sources  the  sidelobes  arc  so  bad  that 
this  technique  is  not  recommended. 

The  procedure  for  calculating  this  equivalent  antenna  pattern  and  processing  the  data  to  pro¬ 
duce  such  a  map  is  given  in  this  section.  This  procedure  may  be  thought  of  as  aperture  synthesis 
with  only  two  fixed  antennas  involved.  An  indication  of  how  to  generalize  these  results  for  any 
number  of  antennas  and  spacings  (true  aperture  synthesis)  is  presented  at  the  conclusion  of  this 
section. 

The  results  of  this  section  were  used  to  write  a  computer  program  which  was  used  to  cal¬ 
culate  several  equivalent  antenna  patterns.  The  program,  however,  was  not  actually  used  to  re¬ 
duce  any  of  the  interferometer  data  for  the  reasons  indicated  in  Sees.  V-A  and  V-J. 
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2.  Fringe-Amplitude  Plane  Coverage 
From  Rogers  (1967b),  Eqs.  (7.2-51)  and  (7.2-52), 

k  =  N  cos  6  g  sin  (Lg  —  Lp)  (9) 

and 

t  =  N  sin<5^  cos<5g  —  N  cos6R  sin6g  cos  (Lg  —  LR)  (10) 

give  the  coverage  of  the  fringe-amplitude  plane  for  a  fixed-baseline  interferometer,  where  Lg 
and  are  the  local  hour  angles  of  the  source  and  baseline,  respectively,  and  6 g  and  6^  are  the 
declinations  of  the  source  and  baseline,  respectively.  N  is  the  baseline  length  in  wavelengths. 
These  equations  may  be  written  as 


k  =  A  sin  t 

(11) 

l  =  B  —  C  cos  t 

(12) 

A  =  N  cos  <5  £ 

(13) 

B  =  N  sin  <5  ^  cos  <5 g 

(14) 

C  =  N  cos  <5  p  sin  6 g 

(15) 

t  =  LS-LB 

(16) 

These  equations  represent  an  ellipse  in  the  kf  plane  as  shown  in  Fig.  5.  The  whole  ellipse,  which 
corresponds  to  t  running  through  a  complete  cycle,  can  be  observed  only  for  a  circumpolar 
source. 

Several  special  cases  are  worth  noting  in  passing.  If  <5^  =  0,  which  corresponds  to  an  east- 
west  baseline,  B  =  0,  and  the  ellipse  is  centered  on  the  origin.  In  this  case,  as  we  shall  see 
below,  it  is  sufficient  to  take  only  12  hours  of  data  instead  of  the  full  24  hours.  If  6g  =  0,  C  =  0, 
and  the  ellipse  degenerates  into  a  horizontal  line  corresponding  to  an  unending  series  of  lobe 


Fig.  5.  Typical  interferometer  fringe 
plane  coverage. 
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ambiguities  in  declination,  or,  in  our  terminology,  an  equivalent  antenna  pattern  which  has  a 
series  of  zero  decibel  sidelobes  in  the  declination  direction.  If  6g  =  6g,  we  have  B  =  C  and  the 
ellipse  goes  through  the  origin.  This  case  has  special  significance  for  sources  with  assumed 
symmetry.  And  finally,  if  =  6g  =  0,  B  =  C  =  0  and  the  ellipse  becomes  a  horizontal  line 
through  the  origin. 

3.  Equivalent  Antenna  Pattern 

The  equivalent  antenna  pattern  of  any  instrument  is  the  response  of  that  instrument  to  a  point 
source  at  the  origin  of  coordinates.  For  an  interferometer,  the  data  are  usually  reduced  in  such 
a  way  that  a  point  source  at  the  reference  point  on  the  sky  produces  a  fringe  amplitude  of  unity 
and  a  fringe  phase  of  zero  for  all  local  hour  angles. 

If  f(x,  y)  represents  the  sky  brightness  temperature  inferred  from  the  data,  and  F(k,f)  rep¬ 
resents  the  fringe  amplitude  and  phase,  then  the  two  are  connected  by  a  two-dimensional  Fourier 
transform 


«*.y>  ■  4-.  f  r«'11"'1  F(k,i)  dkdi 

Z7r  J-OO  QO 


(17) 


The  integral  in  practice  can  be  carried  out  only  over  the  range  of  k  and  i  to  which  the  instru¬ 
ment  responds.  If  we  take  F  to  be  zero  for  those  values  of  k  and  /  which  cannot  be  measured 
and  a  constant  weighting  function,  then  our  result  will  correspond  to  the  "principal  solution"  of 
Bracewell  [Bracewell  and  Roberts  (1954),  Bracewell  (1958)]. 

We  therefore  integrate  along  the  ellipse  of  Fig.  5 


where 


f(x,y)=27  \  e‘l(kx+/y)  K(t)  2*R(t)  o(t)  dt 


R2 (t)  =  tZ  +  k2  =  A2  sin2  t  +  B2  -  2BC  cos  t  +  C2  cos2  t 


(18) 


(19) 


and  the  other  quantities  are  indicated  in  Fig.  6.  We  think  of  o  tending  to  zero  and  K  to  00  while 
their  product  remains  finite.  The  product  KcrR  represents  the  data  (which  we  have  taken  to  be 


Fig.  6.  Area  differential  in  the  k-1  plane. 


unity  for  all  t)  times  an  arbitrary  weighting  function.  This  arbitrary  weighting  function  is  anal¬ 
ogous  to  the  feed  illumination  for  a  single  antenna.  We  have  rather  limited  control  over  the 
equivalent  antenna  pattern  by  varying  this  weighting  function. 
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If  we  define 


W(t)  =  K(t)  c( t)  R(t)  (20) 

then  we  may  write 

p27T 

f(x,  y)  =  2  ft*  \  exp  [—  i(xA  sin  t  +  yB  —  yC  cos  t)J  W(t)  dt  .  (21 ) 

J0 

Note  that  the  integral  in  Eq.  (21)  will  be  complex,  but  f,  since  it  represents  a  brightness  or 
brightness  temperature,  must  be  real.  The  problem  is  that  we  have  taken  the  integral  only  over 
a  region  of  the  ki  plane  that  is  not  symmetrical  about  the  origin.  If  we  demand  that  f  be  real, 
we  have  from  Eq.  (17)  that 

F*  (—  k,  — i)  =  F(k,f)  (22) 

where  *  represents  complex  conjugation.  In  effect,  we  have  measured  the  fringe  phase  and 
amplitude  not  only  along  the  original  ellipse,  shown  in  Fig.  5  by  a  solid  line,  but  also  along  an¬ 
other  ellipse,  obtained  from  the  first  ellipse  by  reflection  through  the  origin.  This  second  el¬ 
lipse  is  shown  in  Fig.  5  by  a  dashed  line.  If  we  add  the  integrals  along  the  two  ellipses,  we  ob¬ 
tain  the  expression  in  Eq  (21). 

Now  if  we  substitute,  according  to  Fig.  7 


s2  =  (XA)2  +  (yC)2 

(23) 

and 

xA  =  S  cos  0 

(24) 

and 

yC  =  S  sin  0 

(25) 

then 

xA  sint  —  yC  cos  t  = 

S  cos  0  sin  t  —  S  sin  0  cos  t  -  — S  sin  (0  —  t) 

(26) 

and  we  get 

f(x,  y)  =  2  e"iyB  ( 

'  2  7T 

exp[+iS  sin  (0  —  t)]  W(t)  dt 

n 

(27) 

For  W(t)  =  1,  this  integral  can  be  done: 

f.  (x*  y)  =  2  ft*  \  exp[iS  sin  (0  —  t)]  dt  =  4n  cos  (yB)  J  (S)  .  (28) 

l  JQ  o 


Fig.  7.  Typical  point  on  the  source 
in  the  xA-yC  plane. 
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So  for  this  case  we  have  an  equivalent  antenna  pattern  with  a  width  to  half-power  of  S  =  1.52 
(along  the  x  axis)  and  a  first  (negative)  sidelobe  of  —0.403  at  S  =  3.83,  y  =  0.  We  shall  see  later 
what  this  means  for  certain  specific  values  of  A  and  C. 

We  return  now  to  the  general  case  in  Eq.  (27)  and  assume  that  W(t)  is  expanded  in  a  Fourier 
series, 

k 

W (t)  =  X  Mn  elnt  (29) 

n=—  k 


W(t)  and  may  be  complex.  So  we  get 


f(x 


,  y)  =  2  e  £  M  i  exp  fi  S[sin  (<t>  —  t)  +  nt] }  dt 
n 


or 


f(x,  y)  =  2*.  e_iyB  £  M  ei[S  sin9  +  n(<A  ~  0)1  de 

n 


which  is 


f(x.y)  =  47 r  S.  e"lyB  Yj  MnJn(S)  e 


i  n<p 


(30) 


(31) 


(32) 


So  we  see  that  an  expansion  of  W(t)  in  a  Fourier  series  corresponds  to  an  expansion  of  f(x,  y)  in 
Bessel  functions,  and  this  relationship  may  be  used  in  the  data  processing. 

Partial  coverage  of  the  ellipse  of  Fig.  5  corresponding  to  a  noncircumpolar  source,  may  be 
handled  by  setting  W(t)  =  0  for  that  range  of  t  for  which  the  source  is  not  observable,  and  then 
performing  the  Fourier  series  expansion  of  Eq.  (29)  with  sufficient  terms  to  approximate  the  dis¬ 
continuous  function  It  is  not  correct,  however,  to  determine  the  Fourier  coefficients  to  fit  the 
data  over  the  observable  range  and  allow  W(t)  to  be  nonzero  over  the  unobservable  range,  be¬ 
cause  this  corresponds  to  arbitrarily  introducing  k*  plane  coverage  that  was  not  in  fact  observed 


4.  Sampling  the  Data 

We  know  from  the  sampling  theorem  that  a  band-limited  function  need  be  measured  only  at 
uniform  intervals  equal  to  one-half  of  the  wavelength  at  the  maximum  frequency.  If  W(t)  is  a 
band-limited  function,  then  we  may  specify  At,  the  maximum  time  interval  between  observations, 
and  the  appropriate  k  in  Eq.  (29). 

It  is  immediately  clear  that  W(t)  cannot  be  a  band-limited  function  of  t  unless  F(k,f)  is  band- 
limited,  and  this  can  only  be  the  case  if  f(x,  y)  is  spatially  limited,  i.e.,  only  if  f  =  0  for  S  >  S  . 
But  this  will  always  be  approximately  true  in  practice  because  f  will  be  limited  by  the  beam- 
width  of  a  single  antenna,  if  not  by  the  source  itself. 

So  consider  a  point  source  at  S  ,  4>q,  which  we  take  to  be  a  representative  point  on  the  edge 
of  an  extended  but  limited  source.  Then  since 

/MJO  /»  OO 

F(k,  /)  =  Zn  \  \  f(x,y)  el(kX+yi)  dxdy  (33) 

00  -  OO 
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we  get 


i(kx  +1  y  ) 

F(k,  I)  =  2w  e 


or 

W(t)  =  27 r  exp  [  iy^B  —  iSQ  sin  —  t)] 

for  this  case. 

To  find  whether  this  is  a  band-limited  function  of  t,  we  consider  F(c o)  defined  by 
C2n 

F(u>)  =  \  W(t)  elaH  dt 
J0 


which  is 


_  iy  B  p27r 

F(co )  =  27r  e  °  \  exp  {— i[S  sin  (<£  —  t)  +  cot])  dt 

o  o 


(34) 


(35) 


(36) 


(37) 


or 

i(y  B+0  co)  /*27r  —  i[ S  sinu+cou] 

F(co)  -  2tt  e  °  °\e°  du 

J0 


(38) 


where  u  =  <t>Q  —  t.  This  integral  gives  us 
__  2  i(y  B+< p  c o) 

F(co)  -  —  (2 7r )  e  °  °  J^(SQ)  +  other  terms  for  co  noninteger  (39) 

Because  J  (S  )  does  not  become  precisely  zero  for  all  to  >  co  ,  W(t)  is  not  strictly  a  band-limited 
to  o  o 

function.  But  there  is  a  practical  cut-off  frequency  coq  for  which  J . .  (S  )  <  €  for  all  co  >  coq,  where 
f  may  be  chosen  as  small  as  we  please.  For  e  ~  10  ^  we  have  as  an  approximation  that 

3S  +8 

wo  “  -  and  <40) 

o  o 

This  practical  cut-off  frequency  corresponds  approximately  to  the  naive  cut-off  that  we  may 
calculate  from  Fig.  8  which  shows  an  f(x,  y)  specified  at  the  critical  interval  corresponding  to  the 
maximum  values  of  k  and  l.  We  get  A0  »  7r/SQ  from  Fig.  8  and  for  large  Sq  this  is  just  3/2  of 
the  At  above. 


Fig.  8.  Sampling  in  the  xA-yC  plane. 
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TABLE  IV 

HAYSTACK-MILLSTONE  INTERFEROMETER  ON  CAS  A 

Hordwore  pyrometers 

60  =  43?  2667 
B 

Ld  =  -  153P3083 

B 

N  =  3800  (approximate,  at  1665  MHz) 

Source  fxjrameters 

6  »  58? 6 

s 

x  or  y  C*  3' 
max  max 

System  parameters 

A  =  2767 

B  =  1357 

C  =  2362 

S  ~  3.28 
o 

Conversion  factors 

Units  of  S  to  arcminutes  along  x  oxis 
along  y  axis 

Equivalent  antenna  pattern  parameters 

Width  between  half-pawer  paints  in  declination  (y)  =  3!  1 

in  right-ascension  (x)  =  3177 
Location  of  first  sidelobe  (—0.403):  y  =  0,  x  =  4175 

Data  characteristics 

Moximum  sampling  interval  At 
Maximum  Fourier  component  k 


1440 

3S  +  8 
o 


=  81 


=  18 


J_  10800' 

A  TT 

1  10800' 
C  TT 


1124 

1146 
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5.  Specific  Numeric  Example 


Table  IV  lists  the  relevant  parameters  for  the  source  Cas  A  as  observed  with  the  Haystack- 
Millstone  interferometer  system. 


Fig.  9.  Interferometer  equivalent  antenna  pattern 
(Haystack-Mi I Istane  on  Cas  A  at  1665  MHz). 


RIGHT  ASCENSION  OFFSET  (  Arcmmutes ) 


A  contour  diagram  of  one  quadrant  of  the  equivalent  antenna  pattern  for  this  case  is  shown 
in  Fig.  9. 


6.  Conclusion 

Because  it  is  possible  in  principle  to  recover  the  original  fringe  amplitude  and  phase  data 
from  the  map,  it  is  clear  that  in  presenting  the  data  in  the  form  of  a  map  we  have  not  discarded 
any  information  contained  in  the  data.  However,  we  may  not  have  made  full  use  of  other  knowl¬ 
edge  about  the  source.  For  example,  we  assumed  in  Sec.  IV-C-4  that  the  source  was  spatially 
limited,  but  we  have  not  incorporated  this  assumption  fully  into  the  map.  Specifying  that  f(x,  y) 
shall  be  zero  outside  some  region  of  x  and  y  is  equivalent  to  convolving  F(k,f)  by  the  Fourier 
transform  of  a  function  which  is  zero  outside  the  prescribed  region  of  x  and  y.  This  convolution 
will  smooth  our  original  F(k,f),  which  was  essentially  a  line  function  along  the  ellipses,  into  a 
band.  It  will  also  perform  some  smoothing  along  the  ellipses.  However,  there  is  no  necessity 
for  actually  performing  this  convolution;  it  is  preferable  to  manipulate  f(x,  y). 

The  flux  in  the  map 

S  >  OC  OC 

Flux  =  \  \  f(x,  y)  dxdy  =  F(0,  0)  =  0  (41) 

OO  '  _  OC 

is  zero  because  we  have  not  incorporated  the  information  from  single  antenna  measurements  into 
the  map.  However,,  this  may  be  done  by  adding  a  constant  to  f,  unless  the  source  is  so  large  that 
it  is  resolved  by  the  single  antennas  used. 
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RIGHT  ASCENSION 


Fig.  10.  W3  region  map  made  at  7.  795  GHz 
at  Haystack. 


62°  00*  00" 


Fig.  1 1.  W3  region  map  made  at  15.  5  GHz 
at  Haystack. 
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In  more  general  terms,  we  may  incorporate  any  additional  kl  plane  coverage  by  adding  the 
corresponding  functions  in  the  xy  plane,  provided  that  there  is  no  overlap  in  the  corresponding 
kf  plane  coverages.  This  allows  us  to  do  true  aperture  synthesis  by  adding  a  series  of  individual 
maps  made  by  this  technique. 

V.  SOURCES  -  DATA  AND  PRELIMINARY  INTERPRETATION 

Data  of  one  kind  or  another  were  taken  for  this  study  on  nine  regions  of  the  sky  that  either 
contain  OH  emission  sources  or  are  interesting  in  relation  to  the  problem  of  understanding  the 
anomalous  OH  emission.  These  data  are  presented  in  this  section  together  with  a  discussion  of 
some  of  the  other  published  data  on  these  sources. 

A  number  of  continuum  maps  of  some  of  these  sources  were  made  and  some  essential  char¬ 
acteristics  of  the  maps  are  summarized  in  Table  V.  All  of  the  maps  were  made  with  the  Hay¬ 
stack  antenna  using  the  drift  curve  technique.  A  number  of  spectra  of  the  OH  emission  from 
some  of  these  sources  were  made  and  some  essential  characteristics  of  the  spectra  are  sum¬ 
marized  in  Table  VI.  These  maps  and  spectra  are  presented  in  this  section  under  the  individual 
sources. 

A  number  of  individual  sources  were  separated  out  on  the  continuum  maps,  and  their  perti¬ 
nent  parameters  are  summarized  in  Table  VII.  Similarly,  a  number  of  OH  emission  features 
were  separated  out  on  the  spectra,  and  the  pertinent  parameters  for  these  features  are  summa¬ 
rized  in  Table  VIII.  As  an  aid  in  this  process,  we  used  a  computer  program  that  compares  the 
spectral  data  with  a  model  consisting  of  a  series  of  uniformly  polarized,  Gaussian-shaped  fea¬ 
tures  that  add  incoherently  in  any  region  of  overlap  [Ball  and  Meeks  (1968)]. 

A.  W3  (IC1795/IC1805) 

W3  is  one  of  the  most  intensively  observed  of  the  sources  studied  for  this  report.  Radio 
continuum  observations  have  been  made  by  Akabane,  et  al.  (1967),  Mezger  and  Henderson  (1967), 
Mezger,  et  al.  (1967),  Aikman  (1968),  and  others  [see  Aikman  (1968)  for  a  list].  Excited  hydro¬ 
gen  lines  have  been  observed  in  this  region  by  Mezger  and  Hoglund  (1967),  Dieter  (1967),  and 
Reifenstein  (1968).  The  OH  emission  from  this  region  has  been  observed  with  a  single  antenna 
by  Weaver,  et  al.  (1965),  Weinreb,  et  al.  (1965),  Dieter,  et  al .  (1966b),  Barrett  and  Rogers 
(1966),  Davies,  et  al.  (1966),  Meeks,  et  al.  (1966),  Weaver,  et  al.  (1968),  and  Coles,  et  al. 
(1968),  and  with  an  OH  interferometer  by  Rogers,  et  al.  (1966),  Cudaback,  et  ah  (1966),  Davies, 

et  aj..  (1967),  Moran,  et  al.  (1967),  Moran,  et  aj. .  (1968),  and  by  Raimond  and  Eliasson  (1 969). 

2 

The  TTj  /£,  J  =  l/2  line  of  excited  OH  has  been  detected  from  this  source  by  Zuckerman,  et  al. 
(1968),  the  27t^  J  =  5/2  line  by  Yen,  et  al.  (1969),  and  the  27r ^  J  =  5/2  line  by  P.  Schwartz 
(1969,  private  communication). 

This  region  has  been  photographed  and  studied  optically.  The  W3  region  from  the  Palomar 
Sky  Survey  prints  was  reproduced  in  Rogers,  et  al.  (1966)  and  in  Weaver,  et  al.  (1968).  Ishida 
and  Akabane  (1968)  show  a  photograph  taken  at  the  Tokyo  Astronomical  Observatory.  Aikman 
(1968)  studied  in  some  detail  the  association  of  the  radio  continuum  sources  in  this  region  with 
their  optical  counterparts. 

A  contour  map  made  at  Haystack  at  7.795  GHz  is  shown  in  Fig.  10  and  one  made  at  15.  5  GHz 
is.  shown  in  Fig.  11.  Although  this  is  clearly  an  extended  complex  region,  it  is  divided  into  three 
principal  sources  listed  in  Table  VII.  The  position  of  the  OH  emission  source  is  coincident 
within  the  error  bars  with  W3y  (G133.9  +  1.0).  This  continuum  source  has  been  discussed  in 
some  detail  by  Mezger,  et  al.  (1967)  and  by  Aikman  (1968). 
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TABLE  V 

SPECIFICATIONS  FOR  CONTINUUM  MAPS 

Minimum 

Calculated  * 

Peah 

Nominal 

Center 

Number  t 

Filler 

Effective 

rmt  Nolle 

Negative 

Contour*  » 

Center 

Band- 

Wove- 

Number 

Scant  per 

Beom- 

Integration 

Antenna 

Antenna 

of  Antenna 

Source 

Frequency 

width 

length 

Front 

of 

66 

Dec  1  (notion 

width 

Time 

Temperature 

Temperature 

Temperature 

Figure 

Nome 

(GHz) 

(MHz) 

(cm) 

Feed 

End 

Polorizotion 

Date 

Scan* 

(d*0) 

Step 

(deg) 

(**c) 

CK) 

CK) 

CK) 

10 

W3 

7.  795 

10 

3.846 

Radar 

Mater 

Left 

2/  9/67 

23 

0.025 

1 

0.073 

3.7 

0.02 

0.06 

0.05, 

0.5,  0.05; 

horn 

circular 

0.5  , 

1.0,  0.  1  ; 

1.0  , 

2.0,  0.  2  ; 

2  0  , 

4.0,  10 

n 

W3 

15.5 

1000 

1.934 

Rotatable 

New 

Vertical 

9/19/67 

83 

0.015 

6,  3,  2 

0.036 

3.7 

0.03 

0.21 

0  2  , 

2.6,  0  2  ; 

horn 

TDA 

9/20/67 

2  0  , 

3.2,  0  4 

15 

W28,  W29 

7.  795 

10 

3.846 

Radar 

Mater 

Left 

6/  8/68 

125 

0.03 

2 

0.072 

3.8 

0.02 

0  07 

0.1  , 

1.0,  0.  1  » 

horn 

circular 

6/21/68 

10  , 

2.0,  0.2  ; 

V  4/68 

2.0  , 

3.2,  0  4 

16 

W28 

15.5 

1000 

1.934 

Rotatable 

New 

Vertical 

5/  2/68 

54 

0.015 

3,  2 

0.036 

1.9 

0.029 

0.  17 

0.2 

horn 

TDA 

5/  3/68 

17 

W28 

15.5 

1000 

1.934 

Rototobl* 

New 

Vertical 

5/  2/68 

54 

0.015 

3,  2 

0.072 

3.8 

0.026 

0.06 

0  1  . 

0  2 

horn 

TDA 

5/3/68 

19 

W43 

7.  795 

10 

3.846 

Radar 

Mater 

Left 

7/  9/68 

68 

0  03 

2,  1 

0.072 

1.7 

0.03 

0  09 

0.  1  , 

10,  0  1  ; 

horn 

circular 

1.0  , 

20,  0  2  ; 

2.0  , 

5.4,  0.2 

20 

W43 

15.5 

1000 

1.934 

Rotatable 

New 

Vertical 

5/  8/68 

64 

0.015 

2 

0.036 

1.7 

0.04 

0.  19 

horn 

TDA 

5/  9/68 

6/26/68 

22 

W44 

7.795 

10 

3  846 

Radar 

Mater 

Left 

6/  6/68 

58 

0  03 

2 

0.079 

3.8 

0.02 

0  07 

0-  1  , 

10,0  1  ; 

horn 

circulor 

6/  7/68 

1.0 

1  8,  0.2 

23 

W44 

15.5 

1000 

1.934 

Rotatable 

New 

Vertical 

5/  8/68 

45 

0.015 

1 

0.036 

0.9 

0.05 

0.21 

0.2  , 

1.2,  0.2 

horn 

T0A 

25 

NGC  6334 

7.  795 

10 

3.846 

Radar 

Mo**r 

Left 

4/23/68 

45 

0.03 

2 

0.073 

4.3 

0.018 

0.18  * 

0.2  , 

2.0,  0.2  ; 

horn 

circulor 

6/22/68 

2.0  , 

5.0,  0.5 

26 

NGC  6334 

7. 795 

10 

3.846 

Radar 

Mater 

Left 

4/23/68 

45 

0.03 

2 

0.073 

4.3 

0.018 

0.  18  * 

0.2  , 

2.0,  0.2  j 

horn 

circulor 

6/22/68 

2.0  , 

5.0,  0.5 

27 

NGC  6334 

15.5 

1000 

1.934 

Rotatable 

Old 

Verticol 

3/4/66 

53 

0.015 

2,  1 

0.036 

1.1 

0.  13 

0.28 

0.2  , 

1.0,  0  2 

horn 

TO  A 

4/  3/66 

38 

W49 

B 

1000 

3.747 

Clavin 

T0A 

Vertical 

4/  4/66 

20 

0  03 

4,  3 

0.073 

5.3 

0.035 

0.07 

0  1  , 

1.0,  0. 1  ; 

9/17/66 

1.0  , 

2.0,  0.2  S 

2  0  , 

3.5,  0,5 

39 

W49 

15.5 

1000 

1.934 

Rotatable 

Old 

Verticol 

2/27/66 

27 

0.015 

5,  3,  2 

0.036 

1.8 

0.  105 

0.  16 

0.2  , 

2.2,  0.2 

horn 

TDA 

4/3/66 

43 

W51 

8 

1000 

3  747 

Clavin 

T0A 

Verticol 

9/16/66 

18  ’ 

0.03 

2 

0.072 

3.6 

0.042 

0.07 

0  2  , 

1.0,  0.2  ; 

1.0  , 

3.0,  05  j 

3.0  , 

6  0,  1.0 

44 

W51 

15.5 

1000 

1,934 

Rotatable 

Old 

Verticol 

10/  3/66 

22 

0.015 

3,  2 

0  036 

1.8 

0.  105 

0  24 

0.4  , 

2.0,  0  4  ; 

horn 

TDA 

20  , 

4  0,  0.5 

46 

W75 

7.  795 

10 

3.846 

Radar 

Moter 

Left 

7/18/68 

130 

0  03 

7-2 

0.072 

4  7 

0  017 

0.  1 

0.05, 

0.  1; 

horn 

circular 

7/19/68 

(0.05)  * 

0  2  , 

1.0.  0.2  ,* 

1.0  , 

3.0,  0  5 

47 

W75 

15.5 

1000 

1.934 

Rotatable 

New 

Verticol 

8/  4/68 

193 

0.015 

6-2 

0  036 

2.3 

0  033 

0.  16 

0.2  , 

1.0,  0  2  ; 

ham 

TO  A 

(0. 0B)  * 

1.0  , 

1.8,  0.4 

*  The  coordinate*  an  all  the  maps  or*  epoch  1950. 

t  For  the  central  region 

,  lometimei  the  edge*  hove  fewer  icon*. 

t  Approximate  value*,  bated  on  the  number*  ir 

>  Table  III. 

LTh*  notation  o.,o~,o„ 

1  2  xi 

j,  ,  where  o 

ttand*  for  fh*  number*  in  thi*  column,  meant  that  contour*  or*  plotted  from 

a,  taa2 

in  ttepi  af  a^.  The  unitt  are  *K  of  antenna  temperature. 

T  But  four  of  the  team  u*#d  only  0.  5-GHz  bandwidth. 

1  The  not**  i*  larger,  probably  becoute  of  the  low  elevation. 

#Th*  number  in  porenth#***  oppliei 

to  the  top  holf  of  the  mop  olanr 
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TABLE  VI 

SPECIFICATIONS  FOR  STOKES  PARAMETER  PLOTS 

Painting  Potition 

Tran*  it  ion 
(MHz) 

a  (1950) 
hr  min  lec 

6  (1950) 
deg  min 

i11 

b" 

^Syitem 

CK) 

Number  of  Rune 

Figure 

Source  Nome 

lec 

(<Ug) 

(deg) 

Antenno 

Feed 

Polarization 

Mode 

Date 

Total 

At  Each  v 

13  t 

W3 

1665 

2  23  16.3 

61  38 

57 

133.949 

1.064 

Hoyt  tack 

Clovin 

Unecrr 

370 

Dkke 

11/29/65 

11/30/65 

15 

7 

Hoyt  toe  k 

Andrew 

Circular 

370 

Dlcke 

12/  7/65 
3/  7/66 
3/12/66 
10/  6/66 

15 

6 

14  t 

W3 

1665 

2  23  16.3 

61  38 

57 

133.949 

1.064 

Hayttock 

Clovin 

Linear 

260 

Dkke 

3/  3/69 
3/  6/69 

29 

6 

Hayttock 

Andrew 

Circular 

250 

Dkke 

9/21/67 

9/23/67 

18 

7 

18  t 

W28 

1720 

17  58  48 

-23  17 

6.76 

-0.27 

Mlllitone 

Horn 

Uneor 

260 

Unew  itched 

4/  9/68. 

6 

6 

Mlllttone 

Horn 

Circular 

260 

Drew  Itched 

3/30/68 

3/31/68 

6 

6 

21 

W43 

1612 

16  45 

-1  45.5 

30.99 

0.08 

Miiiitone 

Horn 

Circular 

230 

U  new  itched 

7/11/68 

14 

7 

24  * 

W44 

1720 

18  53  34 

1  14 

34.63 

-0.45 

Millitone 

Horn 

Uneor 

260 

Drew  itched 

4/  9/68 

4 

4 

Millstone 

Horn 

Circular 

260 

Untwitched 

3/30/68 

3/31/68 

6 

6 

28  f 

NOC  6334 N 

1665 

17  17  33.5 

-35  45 

35 

351.37 

+0.63 

Hayttock 

Clovin 

Llneor 

420 

Dlcke 

9/  2/66 
9/23/66 
10/  2/66 

15 

14 

Hayttock 

Andrew 

Circular 

420 

Dkke 

10/  4/66 
10/  7/66 

7 

5 

29  1 

NGC  6334S 

1665 

17  16  50.1 

-35  53 

17 

351.  18 

+0  67 

Hayttock 

Clovin 

Linear 

420 

Dlcke 

9/21/66 
9/23/66 
10/  2/66 

14 

6 

Hayttock 

Andrew 

Circular 

420 

Dkke 

10/  4/66 
10/  7/66 

10 

8 

30  * 

NGC  6334N 

1667 

17  17  33.5 

-35  45 

35 

351.37 

+0.63 

Hayttock 

Clovin 

Linear 

420 

Dkke 

9/  2/66 
9/30/66 
10/  1/66 

10 

9 

Hayttock 

Andrew 

Clrcuior 

420 

Dicke 

10/  5/66 
10/  8/66 

12 

6 

31  t 

NGC  A334S 

1667 

17  16  50.1 

-35  53 

17 

351.  18 

+0. 67 

Hayttock 

C  levin 

Uneor 

420 

Dicke 

9/24/66 
9/30/66 
10/  1/66 

10 

5 

Hayttock 

Andrew 

Circulor 

420 

Dicke 

10/  5/66 
10/  8/66 

11 

8 

32 

NGC  6334N 

1665 

17  17  32.1 

-35  44 

15 

351.  16 

+0.64 

Hayttock 

Andrew 

Clrcuior 

260 

Unew  Itched 

5/  1/68 
5/20/68 

10 

6 

33 

NGC  6334 S 

1665 

17  16  36.3 

-35  54 

57 

351.  13 

+0.70 

Hayttock 

Andrew 

Circulor 

260 

Unew  Itched 

5/  1/68 
5/20/68 

14 

6 

34 

NGC  6334N 

1667 

17  17  32.1 

-35  44 

15 

351.38 

+0.64 

Hayttock 

Andrew 

Circular 

260 

Untwitched 

V  '/68 

5/20/68 

6 

4 

35 

NGC  6334 S 

1667 

17  16  36.3 

-35  54 

57 

351.  13 

+0.74 

Hayttock 

Andrew 

Circular 

260 

Untwl  tched 

5/  1/68 
5/20/68 

7 

5 

36 

NGC  6334S 

1667 

17  16  36.3 

-35  54 

57 

351.  13 

+0.70 

Millttone 

Horn 

Circulor 

230 

Unew  Itched 

1/21/68 

1/24/68 

12 

8 

40  * 

W49 

1665 

19  07  54 

9  01 

46 

43.  19 

-0,00 

Hayttock 

Clovin 

Uneor 

360 

Dkke 

7/5/66 
7/  6/66 

24 

6 

Hayttock 

Andrew 

Circulor 

370 

Dicke 

3/12/66 
10/  6/66 
10/  7/66 
10/  0/66 

16 

4 

41  » 

W49 

1667 

19  07  54 

9  01 

46 

43.  19 

-0.00 

Hayttock 

Clovin 

Uneor 

360 

Dicke 

4/12/66 

4/13/66 

3/13/66 

12 

6 

Hayttock 

Andrew 

Circulor 

370 

Dicke 

9 

2 

42 

W49 

1720 

19  07  54 

9  00 

38 

43.  19 

-0.00 

Millttone 

Horn 

Circular 

260 

Unew  itched 

7/16/66 

4 

4 

45  t 

W51 

1665 

19  21  33 

14  23 

49  48 

-0  38 

Hayttock 

Andrew 

Uneor 

380 

Dicke 

9/20/66 

9/22/66 

11 

9 

Hayttock 

Andrew 

Circulor 

380 

Dicke 

10/  5/66 

6 

6 

48 

W75N 

1665 

20  36  51 

42  25 

30 

81-86 

0.76 

Millttone 

Horn 

Circulor 

200 

Unew  itched 

7/22/68 

6 

2 

49 

W75S 

1665 

20  3  7  28 

42  12 

81.75 

0.54 

Millttone 

Horn 

Circular 

200 

Unew  itched 

7/22/68 

3 

3 

•  The  pointing  petition  li  the  pet  1  tl on  octually  used  for  the  observe  tlons,  bu*  ii  not  neeettorily  the  boil  known  petition  for  the  OH  emi 

ittion.  Refer  to  the  text  for  ea 

ch  source. 

f  The  S 
'  a 

curvet  ore  from  the  linear  data,  the  curvet  from  the  circular, 

t  The  plat  ii  the  average  af  circular  and  linear. 
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TABLE  VIII 

OH  EMISSION  FEATURE  SUMMARY  * 

Velocity  * 
(km/sec) 

S  Peok 
o 

Flux 

Feature  Width  ® 

Instrumental 

Resolution 

Pc 

\ 

Center  Polorizotion 
Position  Angle 

Nome 

(MHz) 

No. 

<fu) 

(kHz) 

(km/sec) 

(kHz) 

(km/sec) 

(percent)  ' 

(percent)  * 

(deg) 

Remork* 

W3 

1665 

1 

-49.1 

±0.1 

54  ±5 

1.5  ±0.5 

0.27  ±0.1 

1 

0.18 

-50  ±  15 

14  ±  U 

170  ±  90 

2 

-46.4 

±0.1 

137  ±  15 

3.7  ±0.5 

0.68  ±0.1 

1 

0.18 

-94  to -100 

<10 

- 

3 

-45.4 

±0.1 

170  ±20 

3.1  ±0.5 

0.55  ±0.1 

1 

0.18 

-95  to  —100 

17*  5 

65  ±40  1 

4 

-45.1 

±0.  1 

245  ±  X 

25*0.5 

0.45  ±0.1 

1 

0.18 

+91  ±  4 

55  ±  12 

55  *  12 

Blended 

5 

-44.6 

±0.1 

78  ±  9 

2.8  ±0.5 

0,50  ±  0. 1 

1 

0.18 

—  X  to  —IX 

34  ±  10 

105  ±  40  1 

6 

-43,7 

±0.1 

147  ±  18 

1 .9  ±  0.5 

0.35  ±0.1 

1 

0.18 

+83  ±  5 

21  ±  5 

165  ±60 

7 

-43.1 

±0.1 

39  ±  5 

1.4  ±0.5 

0.26  ±0.1 

1 

0.18 

+60  to  +1X 

24  ±20 

90  ±  90 

8 

-41,7 

±01 

69  ±  8 

2  0  ±  0,5 

0,36  ±0.1 

1 

0.18 

+80  to  +1X 

<15 

- 

W28 

1720 

1 

10.8 

±0.1 

119  ±  10 

5.0  ±0.5 

0.9  ±0.1 

1 

0.17 

<8 

<8 

_ 

2 

9.4 

±0.1 

72  ±  7 

5.0±  0.5 

0.9  ±0.1 

1 

0.17 

<16 

<12 

- 

3 

5,2 

±0.4 

17  ±  7 

14  ±5 

2.4  ±0.9 

3 

0.52 

- 

- 

- 

W43 

1612 

1 

2 

40.4 

27.4 

±0.2 

±0,4 

83  ±  9 
X±  10 

6  5  ±0.5 
7.7  ±0,6 

1.2  ±0,1 
1.4  ±0.1 

1 

1 

0.19 

0.19 

<7  t 
<70 

Probably  o  double  feotur* 

W44 

1720 

1 

44.0 

±0.3 

24  ±  2 

9,0  *  1 

1.6  *0.2 

3 

0,52 

<X 

<24 

- 

2 

46.3 

±0.3 

18  ±  2 

6,0  ±  1 

1.0  ±0.2 

3 

0.52 

<X 

<X 

- 

NGC  6334N 

1665 

1 

-12.7 

±0.1 

54  ±  5 

2.2  ±  0.5 

0.40  ±  0.1 

1 

0.18 

+80  to  +1X 

<20 

“  J 

Blended  into  a  tingle 

2 

-12  3 

±0.1 

115  ±  10 

1.8  *  0.5 

0. 33  ±  0. 1 

1 

0.18 

+90  to  +1X 

<10 

-  I 

feotur* 

3 

-9.0 

±0  1 

2X±20 

1  8*0.5 

0.32  ±0.1 

1 

0,18 

-95  to -IX 

16  ±  5 

90  *  10 

4 

-8.1 

±0.1 

29  ±  4 

1.4  ±0.5 

0.25  ±  0.1 

1 

0.18 

-X  to  -  IX 

26  *  15 

175  *  90 

NGC6334N 

1667 

1 

-12.0 

±  0,1 

27  ±  4 

3.0  ±  0.5 

0.54  ±0.1 

1 

0  18 

+75  to  +1X 

40  ±40 

60  ±  90 

2 

-10.6 

±0.1 

29  ±  4 

2  5  ±0.5 

0.45  ±0.1 

l 

0.18 

-50  ±  20 

<50 

-  1 

Blended 

3 

-10.1 

±0.1 

59  ±  6 

1.7  ±0.5 

0.31  ±0.1 

1 

0.18 

-73  to -IX 

18  ±  18 

0  ±  90  I 

NGC6334S 

1665 

1 

-11.8 

±0.1 

29  ±  4 

1.8  ±0.5 

0.32  ±0.1 

1 

0.18 

+70  to  +1X 

50  ±X 

15  *  60  I 

2 

-10.5 

±0.1 

44  ±  4 

3.0  ±0.5 

0.55  ±  0.1 

1 

0.18 

0  to  +44 

60  to  IX 

IX  ±20 

3 

-9.7 

±0.  1 

22  ±  4 

3.3  *  0.5 

0.60  ±0.1 

1 

0.18 

-20  to -IX 

60  to  IX 

X  ±  40 

The  vorioble  feotur**  or* 

4 

-8.1 

±0.1 

49  ±  5 

1.7  ±0.5 

0.X±0.1 

1 

0.18 

+49  ±  15 

26  ±  15 

150  ±  X 

omitted 

5 

-6.0 

±0,1 

61  ±  5 

1.8  ±0.5 

0.32  ±0.1 

1 

0.18 

-25  to  +9 

X  to  IX 

10  ±  10 

6 

-5.6 

±0.1 

34  ±  4 

1.7  ±0.5 

0.30  ±  0.1 

1 

0.18 

+62  to  +1X 

64  ±X 

IX  ±60 

NGC6334S 

1667 

- 

- 

- 

_ 

- 

- 

- 

- 

- 

The  major  feotur**  or* 

time  varying 

W49 

1665 

1 

5.2 

±  0.1 

59  ±  6 

3  9  ±0.5 

0.70  ±0.1 

1 

0.18 

-75  *  25 

<17 

- 

2 

7.7 

±0.1 

54  ±  5 

3.1  ±0.5 

0.55  ±0.1 

1 

0.18 

+75  ±  28 

<19 

- 

3 

11.1 

±0.1 

34  ±  5 

2.5  ±  0.5 

0.45  ±0.1 

1 

0.18 

-50  *44 

<39 

- 

Blended 

4 

12.0 

±0.1 

108  ±  11 

3.3  *0.5 

0.60  ±0.1 

1 

0.18 

-65  ±  14 

<10 

- 

5 

13.6 

±0.1 

49  ±  5 

5.5  ±0.5 

1.00  ±0.1 

1 

0,18 

0  ±31 

<21 

- 

6 

15.0 

±0.1 

44  ±  5 

4.3  ±0.5 

0.77  ±0  1 

1 

018 

+32  to  +1X 

<23 

- 

7 

15.1 

±0.1 

73  ±  7 

1. 7*0.5 

0.X±0.1 

1 

0,18 

-58'to-lX 

<14 

- 

8 

15.6 

±  0.1 

157  ±  16 

2.5  *0.5 

0.45  ±0.1 

1 

0.18 

-20  ±  10 

<7 

- 

9 

16.0 

±0.1 

2X  ±20 

2.2  ±0.5 

0.40  ±0.1 

1 

0.18 

+X±  8 

<6 

- 

10 

16.6 

±0.1 

200  ±20 

3.6  *0.5 

0.65  *0.1 

1 

0.18 

-84  to -IX 

<6 

- 

11 

16,8 

±0.1 

186  ±  19 

3.0  ±0.5 

0.54  ±0.1 

1 

0.18 

-84  to  -IX 

<6 

- 

8lend*d 

12 

17.4 

±0.1 

83  ±  8 

1 . 4  ±  0.5 

0.25  ±0.1 

1 

0.16 

-20  ±  18 

<13 

- 

13 

17.8 

±0.1 

98  ±  10 

4.4  ±0.5 

0.80*0.1 

1 

0.18 

-70  to -IX 

<11 

- 

14 

18.0 

±0.1 

44  ±  5 

2.5  ±0-5 

0.45  *0.1 

1 

0.18 

+36  to  +1X 

<24 

- 

15 

18.5 

±0.1 

64  ±  6 

5.5  ±0.5 

1,00±0.1 

1 

0.18 

-54  to  -IX 

<17 

- 

16 

20.0 

±0.1 

64  ±  6 

1 .9  ±  0.5 

0.35  *0.1 

1 

0.18 

-50  ±23 

<17 

- 

17 

20.9 

±0.1 

226  ±  23 

2.8  ±0.5 

0.50  ±0  1 

1 

0  18 

-90*  7 

<5 

- 

W49 

1667 

1 

2.0 

±0.1 

74  ±  7 

2.5  ±0-5 

0.45  ±0.1 

1 

0.18 

+50  *  16 

<14 

_ 

2 

2.6 

±0.1 

83*  8 

3.9  ±  0.5 

0.70  ±0.1 

1 

0.18 

-60  ±  14 

<12 

- 

3 

2,8 

±0.1 

39  ±  5 

3,6  ±0.5 

0.65  ±  0.1 

1 

0.18 

+70  to  +1X 

<26 

- 

4 

3.3 

±0.1 

24  ±  5 

3.1  ±0.5 

0.55*0.1 

0.18 

-52  to -IX 

<42 

- 

5 

3.9 

±0.1 

59  ±  6 

9.4  ±0.5 

1  70  ±0.1 

1 

0.18 

+X  to  +1X 

<17 

- 

6 

4.9 

±0.1 

88  ±  9 

7.5  *0.5 

1,35  ±0.1 

1 

0.18 

-87  to -IX 

<11 

- 

7 

5.6 

±0.1 

44*  5 

6.8  *0.5 

1.22  *0.1 

1 

0.18 

+X  27 

<23 

- 

8 

6.2 

±0.1 

24  ±  5 

3.1  ±0,5 

1  55  ±0.1 

1 

0.16 

-52  to -IX 

<24 

- 

9 

7.5 

±0.1 

39  ±  5 

11.1  *  0.5 

2  X  *  0.5 

1 

0.18 

+35  ±X 

<26 

- 

10 

19.0 

±  0.1 

98  ±  10 

3  3  ±0.5 

0.59  ±0.1 

1 

0.18 

-X  to  -IX 

<10 

- 

W49 

1720 

1 

15.22  ±0.1 

48  ±  5 

3.9  *0.5 

0.68  *0.1 

1 

0.17 

-84  to -IX 

- 

- 

In  July  1 968 

W51 

1665 

1 

56.8 

±0.1 

29  *  6 

1.8  ±0.5 

0.32  ±0.1 

1 

0.18 

-35  *  38 

<35 

- 

2 

57.7 

±0.1 

127  ±  13 

1. 9*0.5 

0.35  *0.1 

1 

0.18 

+6*9 

<8 

- 

3 

58.1 

±0.1 

49  ±  6 

2.8  ±0.5 

0.50*0.1 

1 

0.18 

+50  ±  23 

<20 

- 

4 

59.1 

±0.1 

78*  8 

4.1  ±0.5 

0.74  ±0.1 

1 

0.18 

+15  ±  14 

<13 

- 

Blended 

5 

59  9 

±0.1 

20*  6 

3.1  ±0.5 

0.55  ±0.1 

1 

0,18 

-42  to -IX 

<50 

- 

6 

61.4 

±0.1 

44  ±  6 

2.2  ±0.5 

0.40  *0.1 

1 

0,18 

-74  to -IX 

<23 

W75S 

1665 

1 

-0.4 

±0.1 

10*  81 

1 

2.5  ±0.5 

0.45*0.1 

1 

0.18 

-57  to -IX 

<27 

- 

2 

+0.5 

±0.1 

48*  e| 

L 

3.3±  0.5 

0.59*0.1 

1 

0.18 

-78  ±  6 

42  *  7 

2  *  10 

3 

+1.0 

±0.1 

34  ±  8 

r 

2.8  ±  0.5 

0.50*0.1 

1 

0.18 

-75  to -IX 

16  *  12 

X  ±  25 

Blended 

4 

+  1  8 

*0.1 

20  ±  8J 

i 

2.2  ±0.5 

0.40*0.1 

1 

0.18 

+64  to  +IX 

42  *  25 

113*  50 

W75N 

1665 

1 

3.0 

±  0.1 

24  ±  4 

1,4  ±  0.5 

0.25  ±0.1 

1 

0.16 

-82  to -IX 

- 

- 

2 

4.9 

±0.1 

14*  4 

2.2  ±0.5 

0.40*0.1 

1 

0.18 

-75  ±  20 

- 

- 

3 

5.3 

±0.1 

15  ±  4 

2.1  ±0.5 

0.38  ±0.1 

1 

0.18 

+60  to  +1X 

- 

- 

Blended 

4 

5.8 

±0.1 

26  *  4 

1.8  ±0.5 

0.32  ±0.1 

1 

0. 18 

+84  to  +1X 

- 

- 

5 

9.1 

±0.1 

15  ±  4 

2.0  ±0.5 

0  36  *0.1 

1 

0.18 

+64  to  +1X 

- 

- 

6 

11.8 

±0.1 

42  ±  5 

1 .9  ±  0.5 

0.34  ±  0.1 

1 

0.18 

+93  to  +1X 

<14 

- 

7 

12.9 

±0.1 

10*  4 

2.0  ±0.5 

0.36*0.1 

1 

0.18 

+73  to  +1X 

*  The  feotur*  summory  for  W3  i*  token  from  Meeks,  *t  ol,  (1966);  for  W28  ond  W44  from  Boll  ond  Stoelin  (1968);  ond  for  NGC  6334,  W49,  ond  W5I  from  Boll  ond  Meek*  (1968),  with 

occosional  correction*. 

t  But  Pc  rise*  to  X  ±  8  percent  ot  40. 7  km/**c  ood  to  -31  *  13  percent  ot  41 .0  km/sec  on 

the  side  of  this  feature  This  is  c 

i  complex  or  multiple  feature. 

I  Th*  velocities  or*  measured  with  respect  to  th*  local  ttondard  of  rest  (tee  Sec.  IV-B).  Th*  avumed  line  r*st  frequencies  or*  those  givon  in  5*C.  I-A. 

§Th*  width*  or*  th*  totol  width*  b«tw**n  half  pow*r  point*  ond  or*  not  corr*ct*d  for  instrumental  resolution. 

ITh*  error  limit*  on  P  (p*rc*nt  clrculor  pol 

arizotion)  ond  Pi  (percent  linear 

polorizotion)  represent  obout  2  ttondard  divisions  estimoted 

on  the  basis  of  m*o*ur*d  receiver  noise.  These 

limit*  apply  only  to  th*  center*  of  the  feature*.  Pc 

i*  positive  for  right  circular  polarization  ond  negative  for  left  circulor  polorizotion. 

#Thl*  source  1*  time  varying;  the  fluxes  her*  or*  from  the  clrculor  doto  token  In  July  1968.  See  Sec.  V-H. 

35 


Figure  12  shows  a  map  of  the  W3  region  made  by  Mezger  [Mezger,  P.  G.  (1965)  unpublished] 
at  1415  MHz  using  the  NRAO  300-foot  antenna.  The  beamwidth  is  about  10  arcminutes.  The  con¬ 
tours  of  this  map  are  shown  superimposed  on  the  National  Geographic- Palomar  Sky  Survey  plate 
of  this  region. 


Fig.  12.  W3  region  mop  mode  of  1.415  GHz  ond  superimposed  on  the  Notionol 
Geographic- Po I omor  Sky  Survey  Print  of  this  region.  The  NRAO  300-foot  on- 
tenno  wos  used  by  P.G.  Mezger  (1965,  unpublished). 


The  SNR  HB-3  is  located  at  2*115m;  62°30'  (1950),  i.e.,  near  the  upper  right  edge  of  Mezger's 
map,  and  is  about  110  arcminutes  in  diameter  [Poveda  and  Woltjer  (1968)].  An  association  be¬ 
tween  HB-3  and  W3  has  been  pointed  out  by  Caswell  (1967).  Moreover,  Mezger's  map  shows 
what  may  possibly  be  another  SNR  arc  extending  from  2*125m;  60°40'  up  around  to  2*118m;  6l°45' 
(1950).  Studies  of  the  spectrum,  polarization,  and  hydrogen  recombination  lines  from  this  arc 
would  establish  whether  it  is  an  SNR  or  not.  However,  if  this  arc  turns  out  to  be  an  SNR,  it 
shares  with  the  SNR  W28  (see  Sec.  V-B)  the  property  of  having  an  OH  emission  source  along  its 
outer  edge  and  the  property  of  having  a  small-diameter  continuum  source  at  its  center.  On 
Mezger's  map  this  center  source  is  at  2nl61?12;  60°55”  (1950).  W41  is  also  an  SNR  with  a  center 

source  and  associated  OH  emission  [Kesteven  (1968a  and  b),  Goss  (1968)]. 

Figure  13  shows  the  1665-MHz  OH  emission  Stokes  parameters  of  W3  measured  at 
Haystack  in  1965  and  1966  This  is  a  re-reduction  of  some  of  the  data  of  Meeks,  e_t  al.  (1966) 
combined  with  some  more  recent  circular  polarization  data  as  indicated  in  Table  VI.  An  error 
in  the  Meeks,  et  al.  spectra,  namely  reversed  labeling  of  S^  and  S^,  has  been  corrected.  A 
summary  of  the  features  seen  on  these  spectra  is  given  in  Table  VIII.  The  properties  of  the  W3 
region  are  summarized  in  Table  IX. 

Figure  14  shows  the  results  of  a  re-measurement  of  the  same  Stokes  parameters  in  1967 
and  1969.  This  re-measurement  was  partly  in  response  to  the  suggestion  by  Coles,  et  al.  (1968) 
that  the  polarization  properties  of  two  of  the  features  in  the  OH  emission  spectrum  of  W3  had 
changed  in  time  without  any  corresponding  change  in  flux.  Coles,  et  al.  found  substantial  linear 
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Fig.  13.  W3  spectrum  at  1665  MHz;  1965  and  1966  data. 
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Fig.  14.  W3  spectrum  at  1665  MHz;  1967  and  1969  data. 


TABLE  IX 

W3  REGION  -  SUMMARY  OF  PROPERTIES 

Distonce 

Excited  Hydrogen 
Line  Observotions 

21 -cm  Hydrogen 

Line  Observations 

OH  Observations  * 

V 

r 

(km/sec) 

Notes 

V 

r 

(km/sec) 

Notes 

V 

r 

(km/sec) 

Notes 

kpc 

Notes 

2.6 

Mezger  ond 
Hoglund 
(1967) 

-43.5 

-50.4 

Mezger  and 
Hoglund 
(1967) 

-48  to  -37 

Absorption 
Rogers  (1967) 

-50  to  -38 

Emission 

Meeks,  et  ol. 
(1966) 

3.  19 

Dieter 

(1967) 

-41.4 

Dieter 

(1967) 

Weaver,  et  al . 
(1968),  etc. 

3.1 

Reifenstein 

(1968) 

-42.3 

Reifenstein 
(1 968) 

-41  to  -36 

Absorption 
Weover,  et  ol . 
(1968) 

*  The  OH  emission  is  disploced  from  the  continuum  sources  to  which  most  of  these  meosurements  opply. 

polarization  in  the  feature  at  —  46 . 4  km/sec,  whereas  neither  Fig.  13  or  14  shows  any  linear  po¬ 
larization  in  this  feature  to  an  upper  limit  of  1  3  percent  (see  also  Table  VIII). 

Coles,  et  al  also  found  no  linear  polarization  in  the  feature  at  —  43. 7  km/sec.  A  compar¬ 
ison  of  Figs.  13  and  14  shows  that  this  feature  does  indeed  seem  to  show  significantly  less  lin¬ 
ear  polarization  on  the  later  date  and  a  small  frequency  shift  as  well  These  data  are  probably 
compatible  with  the  assumption  that  a  weak  but  strongly  linearly  polarized  feature,  which  was 
blended  on  the  high  velocity  side  of  the  —43.7  km/sec  feature,  has  decreased  or  disappeared. 
Otherwise,  the  spectra  of  OH  emission  from  W3  appear  unchanged. 

Some  salient  features  of  these  spectra  are  worth  noting.  The  features  are  highly  polarized  — 
some  are  likely  to  be  nearly  100  percent  Furthermore,  the  complex  features  near  the  center 
may  be  made  up  of  several  Gaussian  features,  each  100  percent  polarized,  which  add  incoher¬ 
ently  in  the  region  of  overlap  [Meeks,  et  al .  (1966)].  The  percent  linear  polarization,  the  per¬ 
cent  circular  polarization,  the  total  percent  polarization,  and  the  linear  polarization  position 
angle  all  vary  as  a  function  of  velocity  across  the  central  region  of  the  spectra. 

The  kinetic  temperatures  implied  by  the  feature  widths  are  15°  to  700° K,  but  the  brightness 

1 2 

temperatures  implied  by  the  measured  fluxes  and  angular  sizes  are  perhaps  10  °K  [Moran, 
e_t  al.  (1968)].  Each  of  these  numbers  must  be  viewed  with  caution,  however.  The  calculated 
kinetic  temperatures  may  be  too  small  if  maser  amplification  has  narrowed  the  feature  widths 
or  too  large  if  microturbulence  or  some  other  line  broadening  mechanism  has  an  appreciable 
effect.  The  non-Gaussian  shape  of  some  of  the  features  probably  means  that  some  influence 
other  than  kinetic  broadening  is  effective.  And  the  calculated  brightness  temperatures  may  be 
too  small  if  interstellar  scintillation  has  broadened  the  apparent  feature  sizes  or  if  absorption 
has  appreciably  reduced  the  flux,  or  too  large  if  refraction  or  reflection  effects  occur  in  a  ma- 
sering  medium. 

The  pointing  position  used  for  the  spectra  and  listed  in  Table  VI  is  from  Rogers,  et  al.  (1967). 
Raimond  and  Eliasson  (1968)  give  a  slightly  more  accurate  position,  namely  2n23ml 6S8±0.  2s; 

61  °38' 54,,±1"  (1950). 
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Fig.  15.  W28-W29  region  map  at  7.795  GHz 
at  Haystack. 


Fig.  16.  W28  region  map  at  15.5  GHz 
at  Haystack. 
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Fig.  17,  W28  region  map  at  15.  5  GHz 
at  Haystack.  Map  contains  same  data 
as  in  Fig.  16;  however,  data  have  been 
smoothed  ta  a  0P072  equivalent  beam- 
width  (see  Table  V  and  Sec.  IV-A-1). 
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The  early  detection  at  8  GHz  of  the  weak  continuum  source  coincident  with  the  W3  OH  emis¬ 
sion  position  (W3y  in  Table  VII)  suggested  that  lower  frequency  continuum  measurements  would 
be  desirable.  Initially,  it  was  not  known  whether  W3y  was  a  thermal  or  nonthermal  source,  or 
whether  it  might  be  a  representative  of  a  new  class  of  source.  Also,  it  was  not  known  whether 
the  OH  emission  and  continuum  positions  were  truly  coincident  or  only  approximately  so.  A 
series  of  measurements  was  undertaken  with  the  Haystack- Millstone  interferometer  to  measure 
the  continuum  emission  from  W3  at  about  18  cm. 

The  initial  measurements  were  made  with  a  4- MHz  window  approximately  centered  on  the 
1665-MHz  emission  from  W3.  Rogers  (1967b)  shows  that  the  vector  average  of  the  fringe  am¬ 
plitude  and  phase  across  the  window  (for  an  add-subtract-autocorrelation  interferometer)  con¬ 
tains  no  contribution  from  any  narrow  line  emission  but  only  a  contribution  due  to  the  continuum. 
This  mode  of  observation  is  ideal  because  the  emission  line  phase  may  be  used  as  a  reference 
with  which  to  compare  the  average  (continuum)  fringe  amplitude  and  phase. 

The  first  series  of  observations,  which  were  made  in  this  mode,  gave  a  clear  indication  of 
a  continuum  source  with  essentially  no  phase  offset  from  the  phase  of  the  OH  emission  line. 
Subsequently,  other  measurements  were  made  with  a  4-MHz  and  a  1.2- MHz  spectral  window  that 
was  offset  to  exclude  the  OH  emission  line.  These  measurements  gave,  by  contrast,  no  indica¬ 
tion  of  any  continuum  emission  near  the  OH  position  in  W3.  For  a  continuum  source  that  is  small 
in  angular  size  compared  with  the  fringe  spacing  and  coincident  in  position  with  the  OH  emission 
source,  these  measurements  imply  an  upper  limit  of  0.4  flux  unit.  This  upper  limit  is  compat¬ 
ible  with  an  extrapolation  of  the  flux  determinations  of  this  source  by  Mezger,  et  al.  (1967)  and 
Aikman  (1968).  There  is  some  indication  in  the  interferometer  data  of  a  rapidly  changing  con¬ 
tinuum  fringe  amplitude  and  phase  which  is  due  to  continuum  emission  well  away  from  the  ref¬ 
erence  position  (the  OH  emission  position). 

It  is  likely  that  some  unknown  instrumental  effect  caused  the  anomalous  results  that  were 
seen  in  the  initial  series  of  measurements. 

B.  W28 

The  M20-W28-W29  (M8)  region  has  been  studied  in  the  radio  continuum  by  Rieu  (1963), 
Mezger  and  Henderson  (1967)  and  Kesteven  (1968a  and  b)  among  many  others.  Table  VII  lists 
the  four  prominent  continuum  peaks  in  this  region,  but  Fig.  15  also  shows  several  subsidiary 
peaks  and  a  ridge  of  continuum  emission  extending  along  the  Galactic  plane.  In  Table  VII  the 
W28  entry  refers  to  the  SNR  which  appears  in  Fig.  15  as  a  semicircular  arc  some  20  arcminutes 
in  diameter.  In  the  center  of  this  arc  is  a  small  diameter  source  listed  as  W28C  in  Table  VII. 

The  marginal  detection  of  W28C  at  15  GHz  shown  in  Figs.  16  and  17  makes  it  likely  but  not  cer¬ 
tain  that  W28C  has  a  thermal  spectrum.  Low  frequency  maps  have  so  far  failed  to  separate  W28 
and  W28C  because  of  lower  resolution  and  probably  also  because  of  the  relative  spectra.  M20 
on  the  north  and  W29  on  the  south  are  known  to  be  H  II  regions.  The  source  G5.9  —  0.4  was,  un¬ 
fortunately,  not  covered  by  the  15-GHz  map;  however,  it  probably  corresponds  to  source  A^  of 
Rieu  (1963),  who  measured  58  fu  at  1430  MHz.  If  this  identification  is  correct,  thenG5.9—  0.4 
is  a  thermal  source. 

Table  X,  which  summarizes  the  properties  of  these  four  sources,  contains  a  number  of  gaps 
and  questionable  entries.  If  the  distance  determinations  are  all  correct,  then  we  conclude  that 
M20  is  separated  in  distance  from  W28  and  W29  which  are  at  about  the  same  distance.  But  the 
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TABLE 

X 

W28-M20-W29  (M8)  -  SUMMARY  OF  PROPERTIES  * 

Excited  Hydrogen 

21- 

cm  Hydrogen 

Line 

Observations 

Line  Observations 

OH  Observations 

Distonce 

V 

v^ 

V 

Region 

kpc 

Notes 

(km/sec) 

Notes 

(km/sec) 

Notes 

(km/sec) 

Notes 

M20  t 

4.15 

Dieter  (1967) 

17.9 

Dieter  (1967) 

- 

- 

- 

- 

4.2 

Goss  (1968) 

16.4 

Reifenstein 

3.8 

Reifenstein 

(1968) 

(1968) 

1.5 

Miller  (1968) 

W28  * 

3.3 

Aizu  ond 

14.2 

Reifenstein 

7.5 

Absorption,  Menan 

6  to  7.5 

Absorption  at  1612,  1665, 

Tobaro  (1967) 

and  Williams  in 

1667  MHz,  Goss  (1968). 

2 

Goss  (1968) 

Goss  (1968) 

3  ta  17 

Emission  ot  1720  MHz, 

1.3 

Povedo  and 
Waltjer  (1968) 

Goss  (1968),  Boll  ond 
Staelin  (1968),  Coles 
et  oJ.  (1968),  Goss  and 

3.6 

Reifenstein 

Robinson  (1968) 

G5.9-0.4  § 

(1968) 

3.0 

Reifenstein 

10.5 

Reifenstein 

- 

- 

-35  to +30 

1612-MHz  emission  at 

(1968) 

(1968) 

* 

—20.5  km/sec, 

1665  MHz  emission  ot 
—20  to  30 km/sec, 

1667  MHz  emission  ot 
—  35  to  27  km/sec, 

1720  MHz  possible  ab¬ 

sorption  at  9  km/sec. 
(Coles,  e_t  ol.  (1968), 
Turner  (1969)) 

W29  ^ 

1.2 

Clark  (1965) 

3.5 

Mezger  ond 

6.3 

Absorption,  Clork 

10.6 

Weok  absorption,  Goss 

1.4 

Gass  (1968) 

Hoglund  (1967) 

0 

(1965) 

(1968) 

1.1 

Mezger  ond 

Hog  1  u  nd 

4.6 

Lilley,  et  al. 
(1966) 

(1967) 

-4 

McGee  ond 

1.0 

Reifenstein 

Gardner 

(1968) 

(1968) 

1.4 

Miller  (1968) 

3.0 

Reifenstein 

(1968) 

*  Refer  to  Toble 

VII  for  additional  information  about  the  corresponding  continuum  sources. 

t  These  meosurements  probably  belong  here 

,  notwithstanding  that  Dieter  (1967)  gives  her  pointing  position  close  ta  W28. 

t  CTB46,  on  SNR  but  there  is  thermol  emission  in  this 

region  olso;  see  text. 

§The  1612-MHz  OH  emission  position  coincides  within  2  orcminutes  with  the  continuum  peak. 

TM8,  NGC  6523,  AMWW3 7,  CTB46,  the 

Logaan  Nebula,  on  H  II  region 

.  Optical 

velocities:  —7.  1  km/sec,  Courtes  (1960),  or— 3  km/sec, 

Campbell  and  Moore  in  ' 

Wilson  (1953). 

picture  is  not  at  all  clear,  because,  for  example,  the  4.2kpc  distance  of  Dieter  (1967)  and  Goss 
(1968)  may  refer  to  the  ridge  of  continuum  radiation  near  W28  rather  than  M20.  Also  the 
1720- MHz  OH  emission,  presumably  associated  with  W28,  extends  almost  up  to  the  velocity 
found  by  Dieter  (1967)  for  M20.  And  in  W29  the  various  hydrogen  recombination  line  velocities 
listed  in  Table  X  suggest  a  complex  or  multiple  source.  The  1612-MHz  OH  emission  velocity 
(about  —  20.  5  km/sec)  for  G5.9-  0.4  is  not  near  any  of  the  other  velocities  listed  and  corresponds 
to  a  kinematic  location  on  the  far  side  of  the  galaxy.  However,  the  1665-  and  1667-MHz  emis¬ 
sion  extends  over  the  range -35  to  +30  km/sec.  The  center  of  this  velocity  range  is  probably  a 
more  reasonable  estimate  for  the  velocity  of  the  OH  source  [see  Turner  (1969)). 


Fig.  18.  W28  spectrum  at  1720  MHz  (S^). 

Figure  18  shows  the  1720-MHz  emission  measured  writh  the  Millstone  antenna.  The  linear 
and  circular  polarization  properties  of  this  OH  emission  were  also  observed  and  upper  limits 
are  listed  in  Table  VIII.  Using  the  Haystack  antenna,  the  position  of  the  two  strongest  1720-MHz 
OH  emission  features  from  W28  was  measured  to  be  17^58m52s±5s,  —  23°18'±1'  (1950)  and  this 
position  has  been  confirmed  by  Goss  and  Robinson  (1968)  who  also  measured  the  circular  and 
linear  polarization  parameters  with  a  better  signal  to  noise  ratio  than  ours  and  who  confirm  the 
low  percentage  of  polarization. 

The  1720-MHz  OH  observations  of  W28  were  made  in  collaboration  with  Professor  D.  H. 
Staelin  of  M.I.T.  (Ball  and  Staelin  (1968)]. 

C.  W43 

The  continuum  emission  around  W43  is  the  most  complex  studied  for  this  report.  Although 
only  two  sources  in  the  W43  region  are  listed  in  Table  VII,  a  complex  extended  region  with  sev¬ 
eral  subsidiary  peaks  is  revealed  in  the  maps  of  Figs.  19  and  20.  Also  the  emission  continues 
on  northward  from  where  the  mapping  stopped  and  perhaps  on  southward  also. 

The  continuum  emission  from  this  region  has  been  studied  by  numerous  other  workers  in¬ 
cluding  Mezger  and  Henderson  (1967),  Kesteven  (1968a),  and  MacLeod  and  Doherty  (1968).  The 
H  I  radiation  in  this  direction  has  been  studied  by  Clark  (1965)  and  others,  and  the  hydrogen 
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Fig.  19.  W43  region  mop  at  7.795  GHz  at  Haystack.  Fig.  20.  W43  region  map  at  15.5  GHz  at  Haystack. 


recombination  line  radiation  from  this  region  has  been  studied  by  Mezger  and  Httglund  (1967), 
Dieter  (1967),  Goss  (1968),  and  McGee  and  Gardner  (1968).  And  the  OH  emission  and  absorption 
in  this  region  have  been  studied  by  Goss  (1968)  and  Coles,  et  al.  (1968).  Figure  21  shows  the 
1612-MHz  OH  emission  observed  for  this  report. 

In  the  direction  of  W43  we  are  probably  seeing  at  least  three  regions  separated  in  distance 
but  approximately  coincident  on  the  sky.  Table  XI  summarizes  the  properties  of  these  regions. 
The  recent  observations  of  Coles,  et  al.  (1968)  have  elucidated  the  nature  of  the  OH  source  in 
region  II  by  showing  that  the  1665-  and  1667-MHz  OH  spectrum  is  made  up  of  narrow  polarized 
emission  lines  masked  by  strong  absorption  at  about  the  same  velocity. 

The  continuum  radiation  from  W43  is  not  known  to  have  a  nonthermal  component  and  there 
are  no  known  nonthermal  sources  nearby.  The  region  around  the  south  source  (G30.0  —  0.0)  has 
not  been  well  studied.  The  W43  region  clearly  deserves  much  more  work. 

On  11  July  1968,  using  the  Millstone  antenna,  the  40.5km/sec  OH  emission  feature  at 
1 6 1 2  MHz  in  W43  was  mapped  to  determine  its  position  and  angular  size.  These  data  show  that 
this  emission  is  coming  from  18^45I^10i’0r^5;  —  1  ^45!  5±6 *  (1950).  This  position  is  about  15  arc- 
minutes  north  of  the  continuum  maximum  and  does  not  agree  with  the  position  given  by  Turner 
(1969)  for  the  same  feature.  Turner’s  position  for  this  OH  feature  is  essentially  coincident  with 
the  continuum  maximum.  Although  the  Millstone  antenna  occasionally  gives  trouble  with  pointing, 
no  errors  as  large  as  15  arcminutes  have  been  noted  on  other  occasions.  Only  additional  obser¬ 
vations  will  resolve  this  discrepancy. 

D.  W44 

The  W44  (3C392)  region  has  been  studied  in  the  radio  continuum  by  Scheuer  (1963),  Hollinger 
and  Hobbs  (1966),  Kesteven  (1968a  and  b)  and  Milne  [in  Goss  and  Robinson  (1968)]  among  many 
others  The  neutral  hydrogen  absorption  from  this  region  has  been  studied  by  Menon  and 
Williams  [in  Goss  (1968)]  and  the  OH  emission  by  Goss  (1968),  Ball  and  Staelin  (1968),  Goss 
and  Robinson  (1968),  and  Coles,  et  al.  (1968).  Excited  hydrogen  line  emission  from  this  region 
was  not  detected  by  Dieter  (1967),  but  was  found  by  Reifenstein  (1968)  in  the  nearby  source  W44W. 
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Figure  22  shows  a  7.795-GHz  map  and  Fig.  23  a  15.5-GHz  map  of  this  region.  The  central 
extended  source,  W44  proper,  is  an  SNR  and  shows  the  arc  structure  typical  of  most  such 
sources.  However,  the  small  diameter  source  designated  W44W  in  Table  VII  has  a  thermal 
spectrum  and  the  ridge  of  continuum  emission  connecting  the  two  may  also  be  thermal. 
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Fig.  22.  W44  region  map  at  7.  795  GHz  at  Haystack.  Fig.  23.  W44  region  map  at  15.  5  GHz  at  Haystack. 


Table  XII  summarizes  the  properties  of  W44  and  divides  it  into  two  regions  separated  by 
20  or  30  km/sec  in  radial  velocity  and  probably  also  separated  in  angle  as  indicated  in  the  notes 
in  Table  XII. 

The  17  20- MHz  OH  emission  from  W44  is  shown  in  Fig.  24  and  upper  limits  for  the  polariza¬ 
tion  of  these  features  are  listed  in  Table  VIII.  The  spectrum  in  Fig.  24  is  not  identical  with  the 
spectrum  of  Goss  and  Robinson  (1968)  who  list  a  feature  at  43.4km/sec  which  does  not  agree 


Fig.  24.  W44  spectrum  at  1720  MHz  (S^). 
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with  the  44  0  km/sec  feature  in  Fig.  24;  however,  the  46.3km/sec  feature  agrees  well.  Both  the 
antenna  beamwidths  and  the  pointing  positions  are  different,  so  the  most  likely  explanation  of  the 
discrepancy  involves  the  existence  of  two  or  more  features  overlapping  in  velocity  but  separated 
in  position.  Another  possibility  involves  the  existence  of  a  feature  extended  in  angle.  Such  a 
feature  would  appear  relatively  more  pronounced  in  the  Millstone  spectrum  as  compared  with  the 
spectrum  made  with  the  Parkes  210-foot  antenna.  The  recent  observations  of  Turner  (1969)  con¬ 
firm  the  existence  of  spatially  separated  features. 

The  distances  for  W44  listed  in  Table  XII  do  not  agree.  The  distance  determination  of  Poveda 
and  Woltjer  (1968)  is  based  on  their  empirical  relation  between  the  radio  surface  brightness  and 
the  linear  diameter  of  an  SNR.  There  will,  of  course,  be  some  scatter  in  this  relation.  The 
kinematic  distance  may  also  be  in  error  because  of  peculiar  motions.  The  OH  absorbing  cloud 
is  certainly  in  front  of  the  SNR,  so  a  hypothetical  spatial  separation  between  them  does  not  rec¬ 
oncile  the  discrepancy.  For  subsequent  calculations,  Goss'  distance  is  used,  but  definitive  con¬ 
clusions  about  this  region  must  await  further  work.  An  excited  hydrogen  line  measurement  would 
be  particularly  helpful. 

E.  NGC  6334 

NGC  6334,  a  complex  H  II  region  at  a  declination  of  about— 36°,  is  difficult  to  observe  from 
the  latitude  of  Massachusetts.  However,  it  has  proved  to  be  an  interesting  source,  both  in  the 
continuum  and  in  OH  emission,  and  well  worth  the  extra  effort  necessary  to  make  accurate  meas¬ 
urements.  Radio  continuum  observations  of  this  region  have  been  made  by  Haddock,  et  al.  (1954), 
Wilson  and  Bolton  (I960),  Wilson  (1963)  (CTB39),  Kesteven  (1968a),  and  others.  Excited  hydro¬ 
gen  lines  from  this  region  have  been  observed  by  Dieter  (1967),  McGee  and  Gardner  (1968),  and 
Reifenstein  (1968).  The  OH  emission  from  this  region  has  been  observed  with  a  single  antenna 
by  Weaver,  et  al.  (1965),  Rogers  (1967a),  Gardner,  et  al.  (1967),  McGee  and  Gardner  (1967), 
Weaver,  et  al.  (1968),  and  Ball  and  Meeks  (1968),  and  with  an  OH  interferometer  by  Rogers, 
et  al.  (1967),  and  Raimond  and  Eliasson  (1969). 

This  region  is  not  too  heavily  obscured  at  optical  wavelengths  and  optical  photographs  have 
been  shown  by  Rogers,  et  al .  (1967),  Gardner,  et  al.  (1967),  Weaver,  et  aT  (1968)  and  others. 

A  continuum  map  made  at  Haystack  at  7  795  GHz  is  shown  in  Fig.  25  and  an  enlarged  portion 
of  this  map  is  shown  in  Fig.  26.  A  15-GHz  map  made  at  Haystack  is  shown  in  Fig.  27.  The  main 
ridge  of  continuum  emission  was  divided  into  three  principle  sources  on  the  basis  of  the  15-GIIz 
map  (Fig.  27).  These,  together  with  three  other  continuum  sources  from  Fig.  25,  are  listed  in 
Table  VII.  The  Howard  and  Maran  (1965)  catalog  entries  for  this  source  (number  497)  are  con¬ 
fusing  but  seem  to  indicate  that  the  source  designated  NGC  6334£  is  the  component  seen  at  low 
frequencies  and  is  therefore  possibly  nonthermal.  However,  Reifenstein  (1968)  found  an  excited 
hydrogen  line  from  this  source  at  a  velocity  very  different  from  the  velocity  of  NGC  6334  (see 
Table  XIII)  and  this  velocity  implies  that  NGC  6334£  is  much  farther  away.  No  prominent  optical 
object  can  be  seen  at  the  position  of  NGC  6334£  on  a  photograph  taken  with  the  ADH  Baker- 
Schmidt  telescope  in  South  Africa.  It  seems  unlikely  that  NGC  6334£  is  associated  with  NGC  6334 
itself. 

The  OH  emission  in  NGC  6334  comes  from  two  small-diameter  regions  separated  by  about 
15.5  areminutes.  The  northern  region  emits  most  strongly  at  1665  MHz  and  the  southern  region 
most  strongly  at  1667  MHz.  In  Table  VI  and  in  the  Berkeley  papers  these  regions  are  referred 
to  as  NGC  6334N  and  NGC  6334S,  but  in  Ball  and  Meeks  (1968)  they  are  referred  to  as  NGC  63341 
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Fig.  25.  NGC  6334  region  map  at  7.  795  GHz 
ot  Haystack. 
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Fig.  26.  NGC  6334  region  map  ot  7.  795  GHz 
at  Haystack  (enlorged  portion  af  Fig.  25). 
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Fig.  27.  NGC  6334  region  map  at  15.  5  GHz 
at  Haystack. 
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ANTENNA  TEMPERATURE  (*K)  ANTENNA  temperature  CK) 


Fig.  28.  NGC  6334N  restored  spectrum  ot  1665  MHz;  Haystack  dota  of  September  ond  October  1966. 
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Fig.  29.  NGC  6334S  restored  spectrum  ot  1665  MHz;  Haystock  doto  of  September  ond  October  1966. 
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and  NGC  6334II,  and  Gardner,  et  ad.  (1967)  refer  to  them  as  NGC  6334A  and  NGC  6334B.  Ap¬ 
parently  it  has  not  occurred  to  anyone  to  call  them  a  and  0. 

Figures  28  through  31  show  restored  spectra  (see  Sec.  IV-B-4)  of  NGC  6334  taken  with  Hay¬ 
stack  in  1966.  These  data  have  been  presented  in  another  form  and  discussed  in  some  detail  by 
Ball  and  Meeks  (1968).  The  circular  and  linear  data  near —6.6  and  —7.8  km/sec  at  the  south  po¬ 
sition  at  1665  MHz  are  mutually  incompatible  because  a  time  variation  occurred  between  the  dates 
on  which  these  data  were  taken. 

PMgures  32  through  35  show  restored  spectra  of  NGC  6334  taken  with  Haystack  in  1968.  A 
comparison  of  these  two  sets  of  figures  shows  no  evidence  for  time  changes  in  any  of  the  features 
from  the  north  position,  but  the  —  8.9km/sec  feature  at  1667  MHz  from  the  south  position  has  de¬ 
creased  over  this  period.  Figure  36  shows  more  recent  data  taken  at  Millstone  that  also  indicate 
that  the  —9.8  km/sec  feature  has  increased.  The  Millstone  data  were  not  restored  because  the 
accuracy  of  the  pointing  of  the  Millstone  antenna  was  inadequate. 

A  very  considerable  amount  of  additional  data  was  taken  at  Millstone  and  Haystack  and  these 
data  were  examined  carefully  for  the  presence  of  any  additional  time  variations.  For  the  most 
part,  the  features  at  both  positions  are  very  stable  not  only  in  flux  but  also  in  polarization.  These 
unchanged  features  lend  confidence  in  the  ability  to  make  accurate  and  reproducible  flux  meas¬ 
urements.  The  time  variation  situation  is  summarized  below. 

Time  variations  in  the  OH  emission  from  NGC  6334  have  been  reported  by  Dieter,  et  al. 

(1966a  and  b),  Weaver,  et  al.  (1968),  Ball  and  Meeks  (1968),  and  Goss,  et  al .  (1968).  At  Hay¬ 
stack  four  features  have  been  seen  to  vary.  At  1665  MHz,  two  features  at  about  —6.6  and 
—  7.9km/sec  showed  a  dramatic  increase  in  flux  on  one  and  probably  two  occasions.  The  time 
scale  for  the  increase  was  9  days  or  less  on  one  occasion.  Within  a  month  or  so,  the  features 
were  back  to  "normal."  At  1667  MHz,  two  features  at  about  —8.9  and  —9. 8  km/sec  varied  on  a 
longer  time  scale,  the  former  decreasing  and  the  latter  increasing  slowly  over  a  period  of 
months  All  four  of  these  variable  features  are  from  NGC  6334S;  none  of  the  features  from 
NGC  6334N  have  been  seen  to  vary  at  Haystack. 

An  important  question  is  whether  the  polarization  properties  vary  in  time  as  suggested  by 
Goss,  et  a_l.  (1968)  or  whether  the  variation  may  be  interpreted  as  a  change  in  the  flux  of  a  fea¬ 
ture  whose  polarization  properties  remain  constant  in  time  as  suggested  by  Ball  and  Meeks  (1968). 
The  spectra  show  that  the  variable  features  are  blended  with  other  features  which  remain  con¬ 
stant.  Since  these  nearby  constant  features  have  different  polarization  properties  from  the  var¬ 
iable  features,  the  polarization  of  the  resulting  blend  will  change  in  time  even  if  the  variable 
feature  itself  maintains  constant  polarization.  This  question  will  be  answered  with  certainty 
only  with  more  complete  data  preferably  made  with  higher  frequency  resolution. 

The  original  Haystack-Millstone  interferometer  measurements  of  Rogers,  et  al.  (1967)  lo¬ 
cated  NGC  6334N  at  17h17m33S5±2S;  -  35 °45’  35" ±1 0"  (1950).  More  recently  Raimond  and 
Eliasson  (1969)  determined  this  position  to  be  17^1 7m32Sl±0s4;  — 35°44’ 1 5"±1 0"  (1950).  To  re¬ 
solve  this  disagreement  and  to  verify  the  position  given  by  Raimond  and  Eliasson  for  NGC  6334S 
[  17h16m36s3±0s4,  35 ° 54' 57" ±4"  (1950)],  the  Haystack- Millstone  interferometer  was  used  to  ob- 

serve  these  sources  at  1665  MHz  on  3,  4,  and  5  January  1968.  A  120-kHz  spectral  window  was 
used  with  3-kHz  resolution  and  both  antennas  were  set  for  vertical  linear  polarization.  The  ref¬ 
erence  position  for  NGC  6334N  was  the  Rogers,  et  al.  position  above,  and  the  reference  position 
for  NGC  6334S  was  the  Raimond  and  Eliasson  position  above.  The  fringe  phase  for  the  two  strong¬ 
est  features  from  NGC  6334N  and  the  strongest  feature  from  NGC  6334S  is  plotted  in  Fig.  37. 


53 


ANTENNA  TEMPERATURE  (*K)  ANTENNA  TEMPERATURE  CK) 


Fig.  30.  NGC  6334N  restored  spectrum  ot  1667  MHz;  Hoystock  doto  of  September  ond  October  1966. 


Fig.  31.  NGC  6334S  restored  spectrum  ot  1667  MHz;  Hoystock  doto  of  September  ond  October  1966. 
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Fig.  32.  NGC  6334N  restored  spectrum  at  1665  MHz;  Haystack  data  of  May  1968. 
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Fig.  34.  NGC  6334N  restored  spectrum  at  1667  MHz;  Haystack  data  of  May  1968. 


57 


Fig.  35.  NGC  6334S  restored  spectrum  at  1667  MHz;  Haystack  data  of  May  1968. 
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Fig.  37.  NGC  6334,  1665-MHz  Haystack- 
Millstone  interferometer  data. 


The  offset  in  fringe  phase  between  the  two  features  from  the  north  position  is  probably  an 
instrumental  effect  related  to  their  opposite  circular  polarizations.  For  example,  a  20°  offset 
between  the  polarization  position  angles  of  the  two  linear  feeds  could  cause  this  effect;  however, 
the  actual  error  in  setting  the  feeds  to  be  vertical  was  almost  certainly  much  smaller  than  20°. 

Figure  37  also  includes  a  line  representing  the  calculated  fringe  phase  for  a  feature 
coming  from  the  position  given  by  Raimond  and  Eliasson  for  the  north  source.  This  line  illus¬ 
trates  that  it  is  very  difficult  to  measure  declination  offsets  for  southern  declination  sources  with 
the  Haystack-Millstone  interferometer  or  with  any  two-element  fixed  baseline  interferometer 
with  an  approximately  north-south  baseline.  The  measurements  of  Raimond  and  Eliasson  made 
with  a  variable  interferometer  do  not  suffer  from  this  problem.  The  present  measurements  do 
not  allow  a  distinction  between  these  two  positions. 

In  conclusion,  it  is  likely  that  this  discrepancy  results  from  an  underestimation  of  the  dec¬ 
lination  errors  for  NGC  6334  (and  W24  also)  by  Rogers,  et  al.  and  the  positions  given  by  Raimond 
and  Eliasson  are  probably  correct. 

F.  W49 

The  source  W49,  a  combination  of  a  very  large  H  II  region  and  a  supernova  remnant,  has 
been  studied  in  the  radio  continuum  by  Mezger  and  Henderson  (1967),  Mezger,  et  ah  (1967)  and 
others.  The  excited  hydrogen  emission  from  this  region  has  been  studied  by  Mezger  and  Hoglund 
(1967),  Dieter  (1967),  McGee  and  Gardner  (1967),  and  Reifenstein  (1968).  Since  its  discovery, 
referred  to  in  Sec.  I-B,  this  OH  emission  from  W49  has  been  studied  with  a  single  antenna  by 
many  workers  including  Palmer  and  Zuckerman  (1967),  Weaver,  et  al.  (1968),  and  Ball  and 
Meeks  (1968),  and  with  an  OH  interferometer  by  Rogers,  et  al.  (1966),  Cudaback,  et  al. .  (1966), 
Davies,  et  al.  (1967),  Moran,  et  al.  (1968),  and  by  Raimond  and  Eliasson  (1969).  The 
J  l/2  line  of  excited  OH  was  detected  from  this  source  by  Zuckerman,  et  al.  (1968).  The  region 
is  heavily  obscured  in  the  optical,  and  photographs  have  failed  to  reveal  this  source. 


59 


t 


o 

an 


> 

an 

< 

£ 

3 

LO 

I 

z 

o 

o 


o 

I 


c  c 
&-2 
2  o 

■o  £ 

-p  <u 

«s 
1  0> 
—  c 

CN  — , 


C  C 

&  .a 

2  o 

x  £ 

x  Si 

■go 


l213 


u 

<D 

E 


E 


<u 

>“t 


& 


■g~ 

O  00 
u.  <> 
4)  O'* 

E  O 
o  c 

Q-  O 

g  I 
.2  -* 
u 

.12  => 

E  N 


CO  CO 
'O  £ 
CN  CN 


— 1_* 
ol  <u 
<u 

<ul  2 

Si 

h- 

<u  — 

5:  <2 


co 

CM 


CM 

I 


<S. 

u 


■-  -I 

ol 


E  «l 


MO 

CN 


i-i 

O  <u 

II 

3^ 


l\ 

<> 

CN 

I  -o 

8>i 

M  :o 

II 


l\ 

£ 

CN 


f\ 

CN 


cS 


(5 


00 

CN 


c 

J2 

£ 


xt- 

r\ 


CM 

u-j 


M5 

00 


r\ 

<> 

CN 

I  -g 

II 

|X 


<> 

CN 


4> 

2 

Q 


c 

* v  ^ 

3  00 
c  "O 
^  O' 


E 

<u 

X 

O 


60 


A  thorough  discussion  of  most  of  the  known  properties  of  W49  has  been  given  by  Mezger, 
et  al.  (1967),  and  so  these  properties  will  only  be  summarized  here  (see  Table  XIV).  Compo¬ 
nent  A  of  W49  has  a  complex  thermal  spectrum  which  Mezger,  et  al.  (1967)  have  analyzed  into 
two  components,  a  large  low-density  H  II  region  (Al)  and  a  small  high-density  H  II  region  (A2) 
almost  coincident  in  position.  One  of  the  two  OH  emission  positions  coincides  with  A2  and  the 
other  lies  within  Al  along  what  may  be  a  ridge  extending  outward  from  A2.  Figure  38  shows  a 
7.795-GHz  map  and  Fig.  39  a  15-GHz  map  made  at  Haystack.  The  OH  emission  positions  are 
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Fig.  38.  W49  region  at  8  GHz  at  Haystack. 


Fig.  39.  W49  region  at  15.25  GHz  at  Haystack. 


indicated  in  Fig.  39  by  squares.  Figure  40  shows  the  1665-MHz,  Fig.  41  the  1667-MHz,  and 
Fig.  42  the  1720-MHz  radiation  from  W49  as  measured  at  Haystack  and  Millstone.  The  OH  emis¬ 
sion  from  W49  is  notable  for  containing  almost  no  linear  polarization.  The  flux  of  the  1720-MHz 
OH  emission  feature  in  Fig.  29  is  time  varying  as  pointed  out  by  Zuckerman,  et  al.  (1969).  The 
excited  hydrogen  velocity  for  W49  (  +  5  to  +7km/sec)  agrees  within  about  lOkm/sec  with  the  cen¬ 
troid  of  the  OH  emission. 

The  other  continuum  source,  W49B,  is  a  supernova  remnant  at  about  the  same  distance  as 
W4  9A  (14kpc)  as  pointed  out  by  Mezger,  et  al.  (1967).  This  SNR  is  apparently  too  far  from  the 
OH  emission  position  for  any  causal  relationship. 

G.  W51 

The  W51  region  has  been  studied  in  the  radio  continuum  by  Mezger  and  Henderson  (1967)  and 
McLeod  and  Doherty  (1968)  and  others,  and  in  excited  hydrogen  by  Lilley,  et  al.  (1966),  Mezger 
and  Hoglund  (1967),  Dieter  (1968),  and  McGee  and  Gardner  (1968).  The  OH  emission  from  this 
region  has  been  studied  with  a  single  antenna  by  Weaver,  et  al.  (1968),  Ball  and  Meeks  (1968), 
and  others,  and  with  an  OH  interferometer  by  Raimond  and  Eliasson  (1969).  A  summary  of  W5l 
properties  is  given  in  Table  XV. 

Figure  43  shows  an  8-GHz  and  Fig.  44  a  15-GHz  map  of  this  region  made  at  Haystack.  How¬ 
ever,  the  region  over  which  appreciable  emission  exists  is  much  larger  as  shown  in  the  maps  of 
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Fig.  42.  W49  spectrum  at  1720  MHz. 
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Fig.  43.  W51  region  at  8  GHz  at  Haystack. 


Fig.  44.  W51  region  at  15.25  GHz  at  Haystack. 
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Mezger  and  Henderson  (1967)  and  McLeod  and  Doherty  (1968).  The  OH  emission  position  deter¬ 
mined  by  Raimond  and  Eliasson  (1969)  is  shown  by  a  square  on  Fig.  44.  Figure  45  shows  the 
1 665-MHz  emission  from  W5l.  This  emission  is  notable  for  containing  two  nearly  unpolarized 
features  as  noted  in  Table  VIII.  W51  also  shows  OH  emission  at  1612  and  1720  MHz  to  be  about 
half  as  strong  as  the  1665-MHz  emission. 

The  two  continuum  sources  shown  in  Figs.  43  and  44  are  thermal  and  it  is  clear  that  W51  is 
a  large  and  complex  H  II  region.  What  is  not  clear  is  whether  it  also  contains  an  SNR  as  sug¬ 
gested  by  Mezger  and  Henderson  (1967)  and  others.  The  spectrum  does  turn  up  at  low  frequen¬ 
cies;  however,  McLeod  and  Doherty  (1968)  suggest  that  the  spectrum  of  W51  can  be  understood 
in  terms  of  several  thermal  components  becoming  optically  thick  at  different  frequencies.  Shaver 
(1969)  interpreted  his  recent  observations  at  408  MHz  as  indicating  the  presence  of  a  nonthermal 
loop  well  south  of  the  region  mapped  in  Figs.  43  and  44  It  would  be  very  helpful  to  have  a  meas¬ 
urement  of  linear  polarization  on  that  part  of  the  W5l  region  which  is  suspected  of  being 
nonthermal. 

H.  W75 

The  W75  region  has  been  studied  in  the  radio  continuum  by  Pike  and  Drake  (1964),  Downes 
and  Rinehardt  (1966),  Ryle  and  Downes  (1967),  and  others.  The  OH  emission  from  W75  has  been 
studied  by  Palmer  and  Zuckerman  (1967),  Ball  and  Meeks  (1968),  Weaver,  et  al.  (1968),  and 
Zuckerman,  et  al.  (1969).  Raimond  and  Eliasson  (1969)  studied  the  OH  emission  from  this  re¬ 
gion  interferometrically,  but  only  one  of  the  two  OH  emission  sources  near  W7  5  was  then  known. 
A  summary  of  properties  is  given  in  Table  XVI. 


Fig.  45.  W51  spectrum  at  1665  MHz. 


65 


Cl. 

o 

CtC 

Cl. 


C 

o 

■*- 

o 

£ 

0) 

j 

X 

o 


c  c 
fl)  O 
D)  •“ 

S  O 

0) 


E 
u 

-i  s 

CN  3 


c  c 
8.  .2 
2  o 

"5.  ^ 

X  S 

■go 


O 


O 

Z 


E 


o 

Z 


o 

Z 


u 

fl) 

E 


<D 

O  w 
D  _* 

N  8 
c  2 
O  -D 

L-  C 

<D  O 

£  — 


C 

o  r\ 

•«  s 

v>  O' 


CO 

<> 

O' 

^  cP 
o‘l* 
tl_ 

j-  ol 

|ti 

*  c 

*  i 

CO  _y 

« s 

ON 


co 

T 

o 

04 


'  CN 
^  <> 
-  Os 


0 - 1 

8  Z\ 

-0  fl)[ 

< 


CN  N" 
CO  ^ 


<> 

O' 


0) 

% 

b 


N 

0)  o 

5 


c 

4)  ^ 

2  CO 
c  'O 
^fl>  (> 


£' 


00 

N 


CO 

CN 


c 

_  "O  ^ 

N  E  IN 

-  O  sQ 


ol 

ait 

N 

4>  O  “T 

2  £ 


0)  O' 


c 

*  fl)  ^ 

"5  co 

I  £ 
ic 


0“ 

o 


c 

* 

o 

c 

_* 

c 

D 


o 

fl) 


fl) 

o 

c 

a 


"o 

fl) 


"O 

c 

o 


1 


o 

JO 

(5 


c 

o 

“O 

fl) 


2L 


66 


Figure  46  shows  a  7.795-GHz  map  and  Fig.  47  a  15.5-GHz  map  made  at  Haystack.  Because 
this  source  was  used  as  a  calibrator  {see  Secs.  II-B  and  IV-A),  neither  the  flux  nor  the  position 
of  W75-DR2 1  given  in  Sec.  V  can  be  considered  as  independent  measurements.  In  Fig.  46  the 
two  OH  emission  positions  are  shown  by  squares.  The  position  of  the  southern  OH  emission 
source  (W75S)  was  taken  from  Raimond  and  Eliasson  (1969)  and  the  position  of  the  northern  source 
(W75N)  was  measured  at  Haystack  in  May  1968  to  be  20^36m51  ±9  ;  42°25!5±1'  (1950)  (see 
Zuckerman,  et  al_.  (1969)). 
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Fig.  46.  W75  map  at  7.  795  GHz  at  Haystack. 


Fig.  47.  W 75  map  at  15.  5  GHz  at  Haystack. 


In  Fig.  46,  the  contour  that  touches  W75N  is  marginally  significant,  but  even  if  it  is  real,  it 
may  correspond  to  a  ridge  or  extended  emission  north  of  DR21  rather  than  a  point  source  at 
W75N.  Any  8-GHz  continuum  source  that  is  coincident  with  W75N  is  much  weaker  than  the  con¬ 
tinuum  source  coincident  with  the  W3  OH  position  (G133.9  +  1.1)  which  gives  about  0.25°K  of  an¬ 
tenna  temperature  with  Haystack  at  8  GHz  (see  Table  VII).  The  data  in  Fig.  47  show  that  there 
is  no  continuum  source  at  15  GHz  near  W75N  to  a  limit  of  about  0.05°K  of  antenna  temperature. 

Figure  48  shows  the  spectrum  of  W75N  at  1665  MHz  as  measured  with  the  Millstone  antenna. 
This  spectrum  was  not  restored  and  the  features  between -1  and  +1.5km/sec  are  from  W75S 
seen  on  the  side  of  the  antenna  beam.  An  attempted  restoration  of  these  spectra  failed  because 
the  pointing  of  the  Millstone  antenna  was  not  accurate  enough.  This  leaves  open  the  possibility 
that  some  W75N  features  exist  in  the  range  - 1  to  +  1.5km/sec.  W75N  has  a  1°K  left  circular 
feature  at  1667  MHz  at  a  velocity  of  5.4km/sec.  No  features  were  observed  at  either  1612  or 
1720  MHz  to  within  about  the  same  noise  level  and  over  the  same  velocity  range  as  shown  in 
Fig.  48.  There  is  no  linear  polarization  in  the  strongest  feature,  which  is  at  11.9km/sec,  to  an 
upper  limit  of  14  percent. 

Figure  49  shows  the  1665-MHz  spectrum  of  W75S.  The  strongest  feature  in  this  spectrum, 
which  is  at  0.5  km/sec,  is  time  varying  [see  Zuckerman,  et  ad.  (1969)).  Figure  50  shows  the 
details  of  the  spectrum  of  this  source  on  two  particular  occasions.  Despite  the  flux  variations, 
the  linear  position  angle  has  probably  remained  constant  in  time,  although  it  does  vary  as  a 
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Fig.  49.  W75S  spectra  at  1665  MHz  (right  and  left  circular  palarizatian). 


tn 


Fig.  50.  W75S  spectral  details  at  1665  MHz 
for  the  twa  dates  indicated. 
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flux  units 


function  of  velocity  across  the  feature.  The  percent  linear  and  the  percent  circular  polarization 
also  vary  across  the  feature  and  this  variation  makes  it  difficult  to  compare  results  made  with 
different  receivers;  however,  it  is  likely  that  these  polarization  parameters  have  also  remained 
constant  in  time. 

The  Cygnus  X  region  in  which  these  OH  emission  sources  are  located  is  very  complex.  At 
least  some  components  of  the  continuum  emission  are  probably  nonthermal.  Downes  and 
Rinehart  (1966)  concluded  that  DR- 4  has  a  nonthermal  spectral  index,  but  both  Dieter  (1967)  and 
Reifenstein  (1968)  found  excited  hydrogen  line  emission  in  this  source.  However,  Reifenstein 
found  no  109a  line  in  DR  3,  11,  and  12.  However,  Mrs.  Dieter  found  the  158c*  line  in  DR  3. 

These  confusing  and  partly  contradictory  results  may  arise  from  varying  amounts  of  smearing 
by  the  various  antenna  patterns  in  a  region  in  which  both  thermal  and  nonthermal  sources  prob¬ 
ably  exist  along  some  lines  of  sight.  Definitive  results  will  come  from  higher  resolution  studies 
and  polarization  measurements.  However,  it  seems  likely  that  at  least  some  components  of  this 
region  will  prove  to  be  nonthermal. 

Some  of  these  observations  of  W75  were  undertaken  in  collaboration  with  Dr.  Ben  Zuckerman 
[Zuckerman,  et  al.  (1969)]. 


J.  Cas  A 


Cas  A  is  a  relatively  young  supernova  remnant  and  the  brightest  continuum  radio  source  in 
the  sky  at  most  wavelengths.  The  angular  diameter  of  the  radio  source  is  about  4  arcminutes 
and  there  is  considerable  structure.  A  21.3-cm  continuum  map  made  by  Ryle,  et  al.  (1965) 
shows  a  ring  structure  such  as  would  be  seen  from  a  shell  source  with  a  ratio  of  thickness  to 
radius  of  about  0.25.  However,  there  are  numerous  bright  spots,  some  of  which  are  probably 
not  resolved  by  Ryle's  interferometer. 

The  OH  absorption  against  Cas  A  as  observed,  for  example,  by  Rogers  and  Barrett  (1967) 
and  Goss  (1968),  is  in  two  parts,  one  near  zero  and  one  near  -45km/sec  radial  velocity.  Each 
of  these  parts  breaks  up  into  further  components;  there  are  at  least  four  absorption  features  be¬ 
tween  —  50  and  —  30km/sec.  These  features  are  due  to  OH  near  Cas  A  but  the  zero  velocity  ab¬ 
sorption  is  due  to  OH  much  closer  to  the  sun. 

The  problem  of  the  distribution  of  OH  absorption  across  an  extended  source  such  as  Cas  A 
can  be  approached  by  at  least  two  methods.  The  first  method  involves  a  sum  rule  for  absorption 
and  can  be  used  to  get  evidence  for  mottling  in  the  distribution  of  OH  absorption  without  specif¬ 
ically  measuring  that  OH  distribution.  The  second  method  involves  the  use  of  a  spectral-line 
interferometer.  These  two  methods  are  discussed  below. 

Rogers  and  Barrett  (1967)  derived  a  sum  rule  for  OH  absorption  which  says  that  at  a  given 
frequency  offset  (from  the  line  rest  frequencies),  the  opacities  in  the  four  transitions  must  add 
up  as  follows: 


dTl667  dT1665 
9  5 


dTl6l2  +  dT1720 


(42) 


For  local  thermodynamic  equilibrium  (i.e.,  all  four  excitation  temperatures  equal)  we  have  also 
the  stronger  rule 


dTl667  dTl 665 

9  5  dTl 61 2  dTl 720 


(43) 
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Equation  (42)  will  also  hold  for  the  total  opacity  along  a  given  line  of  sight  to  the  source,  and  if 
the  OH  along  all  lines  of  sight  to  the  source  that  are  accepted  by  the  antenna  beam  is  the  same 
in  excitation  temperature  and  density,  then  Eq.  (42)  will  hold  for  the  measured  opacities. 

Rogers  and  Barrett  (1967)  checked  this  sum  rule  for  their  observations  of  the  OH  absorption 
features  near  zero  velocity  in  the  direction  of  Cas  A.  They  showed  that  Eq.  (42)  holds  for  the  ap¬ 
parent  opacities  to  within  limits  set  by  noise  in  spite  of  the  existence  of  emission  at  1720  MHz. 
The  emission  at  1720  MHz  is  associated  with  unusually  strong  absorption  at  1 6 1 2  MHz  at  the  same 
radial  velocity  as  required  by  Eq.  (42).  However,  for  some  other  sources  the  sum  rule  does  not 
hold  for  the  apparent  opacities  [see  Robinson  and  McGee  (1967)].  This  may  be  caused  by  a  mot¬ 
tled  distribution  of  OH  across  the  face  of  the  source.  For  example,  if  a  small  but  dense  cloud 
of  OH  covers  only  a  part  of  the  source,  the  apparent  opacity  (over  the  whole  source)  will  be  less 
than  the  true  opacity  through  the  cloud,  as  discussed  below. 

Let  us  consider  the  simplified  case  in  which  OH  clouds  cover  a  fraction  x  of  the  face  of  the 
so.  'ce.  Then  we  have 

“Tt 

T^  a  T(,  [1-x  +  xe  ]  +  terms  proportional  to  (44) 

where  T^  is  the  measured  antenna  temperature,  T^,  is  the  continuum  source  temperature,  is 
the  true  opacity  of  the  OH  in  a  cloud,  and  is  the  excitation  temperature  of  the  OH.  We  may 

define  an  apparent  opacity  as  that  value  of  which  makes  Eq.  (44)  hold  for  an  assumed  x  =  1. 
Then 


~Tt 

x  ( 1  —  e  )  -  1  -  e 


(45) 


Since  0^  x<  1,  we  underestimate  T ^  if  we  take  for  t^,  and,  for  a  given  x  we  underestimate 
t ^  by  an  amount  which  increases  as  increases.  For  example,  for  x  =  \  and  we  get 

=  0.22,  but  for  =  10,  we  get  -  0.69.  Equation  (42)  will  not  hold  for  the  Ta’s  unless  either 
x  =  1  or  t  «  1.  This  model  is  essentially  a  two-piece  model  —  uniform  OH  over  a  fraction  of 
the  source  and  no  OH  over  the  remainder.  The  real  situation  is  likely  to  be  more  complicated. 

We  now  consider  the  question  of  detecting  the  presence  of  mottling  in  OH  absorption  using 
the  Haystack- Millstone  interferometer.  For  a  continuum  source  overlaid  with  a  uniform  slab 
of  absorbing  OH,  the  fringe  amplitude  across  the  absorption  will  be  proportional  to  the  single¬ 
antenna  spectrum  and  the  fringe  phase  will  be  constant  as  a  function  of  velocity  although  it  may, 
of  course,  vary  as  a  function  of  local  hour  angle.  For  a  mottled  distribution  of  OH,  the  fringe 
amplitude  will  depart  from  the  single  antenna  spectrum  and  the  fringe  phase  will  no  longer  be 
constant.  The  departure  from  constancy  in  the  fringe  phase  will  be  used  as  the  most  sensitive 
test  for  mottling. 

The  problem  of  obtaining  data  on  the  OH  absorption  against  Cas  A  is  discussed  in  Sec.  Ill— C . 
Initial  measurements  were  made  on  31  November  and  1  and  2  December  1967,  and  again  on  3  to 
5  January  1968,  using  a  400-kHz  window  centered  at— 25km/sec  to  cover  both  the  zero  and 
—  45  km/sec  absorption  features.  The  400-kHz  window  was  wide  enough  to  include  several  fre¬ 
quency  fringes  and  the  instrumental  effect  mentioned  in  Sec.  111-C  was  small.  Most  of  the  data 
from  this  series  of  ob  nervations  showed  no  variation  of  phase  with  radial  velocity  except  for  that 
attributed  to  the  instrumental  effects,  but  a  few  runs  at  certain  local  hour  angles  showed  hints 
of  the  sought-for  effect  as  shown  in  Fig.  51. 
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Fig.  51.  Interferometer  fringe  phase  on  Cas  A  for  the  indicated  local  hour  angles. 


A  further  series  of  Haystack- Millstone  interferometer  observations  of  the  OH  absorption 
against  Cas  A  was  undertaken  27  February  to  2  March  1968  using  a  120-kHz  and  a  40-kHz  win¬ 
dow  centered  on  those  regions  of  the  spectrum  which  had  shown  some  hint  of  fringe  phase  var¬ 
iation  at  the  lower  resolution.  For  this  last  series  of  observations,  the  continuum  fringe  com¬ 
pensation  scheme  described  in  Sec.  III-C  was  used.  Although  this  scheme  was  only  partially  suc¬ 
cessful  and  so  the  data  were  difficult  to  interpret,  no  certain  cases  of  fringe  phase  variation  were 
observed.  The  upper  limit  that  these  measurements  place  on  the  mottling  of  the  OH  distribution 
is  discussed  below. 

In  order  to  calculate  an  upper  limit  to  the  mottling  of  the  OH  absorption  against  Cas  A,  it 
is  necessary  to  have  some  model  for  the  OH  distribution.  Aside  from  signal- to-noise  consider¬ 
ations,  there  are  some  spatial  distributions  to  which  the  interferometer  does  not  respond.  These 
are  distributions  whose  transform  is  zero  along  the  ellipses  of  Fig.  5.  Of  the  other  spatial  dis¬ 
tributions,  some  will  be  below  the  limit  of  detectability  because  of  noise.  If  all  the  absorbing 
OH  at  some  frequency  were  concentrated  in  a  spot  small  compared  with  the  interferometer 
fringes,  and  if  this  spot  were  placed  so  that  its  fringe  phase  were  in  quadrature  with  the  contin¬ 
uum  fringe  phase  (at  some  local  hour  angle),  then  the  observed  fringe  phase  offset  would  be  just 
equal  to  the  percentage  absorption  expressed  in  radians,  that  is,  to  about  1.5  degrees  at  most, 
when  the  continuum  fringe  amplitude  is  unity.  However,  the  observed  continuum  fringe  amplitude 
was  as  little  as  10  percent,  at  certain  local  hour  angles,  so  if  the  condition  above  had  existed  at 
such  a  local  hour  angle  the  fringe  phase  offset  would  be  15  degrees.  No  fringe  phase  variation 
was  seen  on  the  narrow  band  runs  to  an  upper  limit  of  about  5  degrees. 

VI.  INTERPRETATION 

A.  OH  Source  Classification  Schemes 

Several  dozen  OH  emission  sources  are  now  known.  These  sources  were  found  by  a  number 
of  different  observers  engaged  in  a  variety  of  programs;  only  two  limited  but  systematic  searches 
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for  OH  emission  sources  over  selected  regions  of  the  sky  have  been  made  [Zuckerman,  et  al. 
(1969),  Elld€r,  et  al.  (1969)],  and  these  were  not,  of  course,  regions  chosen  at  random.  One 
must  be  wary  of  the  contribution  of  observational  selection  to  phenomenological  associations  and 
also  be  wary  of  statistical  fluctuations  in  so  small  a  sample.  Table  XVII  lists  all  the  OH  emis¬ 
sion  sources  known  to  the  author  with  a  classification  scheme  based  on  suggestions  by  Turner 
(1967)  and  Goss  (1968),  and  used  by  Ball  and  Staelin  (1968). 

Class  I  OH  emission  sources  emit  most  strongly  at  1665  or  1667  MHz.  The  features  are 
narrow,  typically  0.6  to  4kHz,  and  are  usually  strongly  circularly,  or  less  commonly,  linearly 
polarized. 

W3  and  W49  are  typical  Class  I  sources  and  the  descriptions  in  Secs.  V-A  and  V-F  can  be 

taken  as  typical  of  this  class.  For  example,  these  sources  have  low  apparent  kinetic  tempera- 

12 

tures  (15  to  700°K)  but  very  high  apparent  brightness  temperatures  (perhaps  10  °K  for  W3). 

The  polarization  of  the  features  is  often  complex  with  the  percent  circular  polarization,  the  per¬ 
cent  linear  polarization,  and  the  linear  polarization  position  angle  varying  as  a  function  of  ve¬ 
locity  across  a  feature.  The  shape  of  most  features  is  complex  (asymmetrical)  suggesting  either 
two  or  more  overlapping  features,  or  a  line  broadening  mechanism  which  creates  asymmetric 
line  shapes.  These  characteristics  are  discussed  further  in  Sec.  VI-C. 

Weaker  emission  in  the  wing  transitions  (1612  and  1720  MHz)  is  sometimes  present  inClass  I 
sources  and  when  present  it  is  also  narrow  and  polarized.  All  known  time  variable  sources  are 
in  Class  I,  but  this  may  be  due  to  the  fact  that  the  Class  I  sources  were  more  carefully  studied. 
Almost  all  Class  I  sources  seem  to  be  associated  with,  but  not  usually  coincident  with  H  II  re¬ 
gions.  One  Class  I  OH  emission  source  (Orion  A)  also  has  an  IH  star  [Raimond  and  Eliasson 
(1967)]  and  recent  reports  indicate  that  others  may  also  have  one  [Neugebauer,  et  al.  (1968)]. 

Many  Class  I  sources  also  have  SNR's  or  NTS's  somewhere  around  them,  but  these  are  far  enough 
away  so  that  one  would  not  suspect  any  connection  if  this  association  did  not  occur  so  frequently. 
Even  if  these  phenomenological  associations  stand  the  test  of  future  observations,  one  cannot 
necessarily  infer  an  interaction  mechanism. 

Class  IIa  OH  emission  sources  emit  most  strongly  at  1720  MHz.  The  features  are  wider, 
typically  4  to  20  kHz,  and  not  strongly  polarized.  There  is  usually  no  emission  at  the  other  tran¬ 
sitions.  There  are  only  a  few  sources  in  this  class,  but  an  association  with  SNR's  seems  likely. 
This  association  has  also  been  pointed  out  by  Goss  and  Robinson  (1968).  In  Sec.  VI-B,  one  pos¬ 
sible  interaction  mechanism  between  SNR's  and  OH  emission  sources  is  considered. 

Class  11^  OH  emission  sources  emit  most  strongly  at  1612  MHz.  The  1612-MHz  features 
from  Class  11^  sources  are  wider  than  the  features  from  Class  I  sources,  perhaps  even  wider 
than  those  from  Class  11^  sources,  and  are  not  strongly  polarized.  Several  of  these  sources 
also  show  emission  at  the  main  transitions  (1665  and  1667  MHz)  and  in  at  least  one  case  (NML 
Cygnus)  the  1665-MHz  emission  is  similar  to  that  from  Class  I  sources  in  being  circularly  po¬ 
larized  and  quite  narrow  in  velocity.  Some  of  the  Class  11^  sources  are  associated  with  IR  stars 
and  some  with  H  II  regions.  But  since  the  IR  star  surveys  are  incomplete,  it  may  turn  out  that 

all  the  Class  II,  sources  have  IR  stars.  The  similarities  between  Class  II,  and  Class  I  lead  one 
b  b 

to  speculate  that  the  differences  may  be  of  degree  and  not  of  kind.  No  Class  11^  sources  are 
known  to  have  1720-MHz  emission. 

This  classification  scheme  does  not  contain  the  "Heiles  Clouds"  [Heiles  (1968),  Cudaback 
and  Heiles  (1969)]  —  regions  emitting  apparently  normal  (i.e.,  LTE)  OH  emission.  These  clouds 
are  associated  with  dusty  regions  and  Heiles  interprets  his  results  to  mean  that  the  OH  in  dusty 
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TA8LE  XVII 

CLASSIFICATION  OF  OH  EMISSION  SOURCES 


Tronution  (MHz) 

Ncarwt 

18-cm 

l" 

b" 

Sourc*  Nom#(i) 

1612 

1665 

1667 

1720 

Clou 

Continuum 

(d+fl) 

(dm) 

Ditto  nc« 

R#f*r*ne«t  t 

CIT3 

IR0103  +  12C 

U 

O 

U 

O 

"b 

128.6 

-50.1 

20,21 

G 1 27.  7  +  14.0 

X 

HC 

127  7 

+  14.0 

5 

W3 

1C  1795 

C,  A 

E.A 

C.A 

C.A 

1 

a 

H  II 

133. 95 

+  1.06 

2  6  to 

3.2 

1,15 

G159.7-  19.6 

X 

HC 

159  7 

—  1?.  6 

5 

NML  Tauru* 

U 

O 

U 

O 

"b 

<20  fu 

178.0 

-31.4 

20,21 

G168.7-  15.5 

X 

HC 

168.7 

-15.5 

5 

G173.3  -  16.3 

X 

HC 

173.3 

-16.3 

5 

GI73.9-  15.9 

X 

HC 

173.9 

-15.9 

5 

Cloud  2 

G174.6-  13.7 

X 

HC 

174.6 

-13.7 

5,12 

G171.4-  10.6 

X 

HC 

171.4 

-10.6 

5 

IR050  7  +  52 

U 

O 

O 

O 

"b 

<20  fu 

156.4 

+  7.8 

21 

G178.9  -6.7 

X 

HC 

178.9 

-6.7 

5 

Orion  A 

c 

C,T 

X 

A 

1 

a 

H  II 

209.0 

-19.4 

17,19 

G214.5-  19.9 

X 

HC 

214.5 

-19  9 

5 

G203.3-  11.2 

X 

HC 

203.3 

-11.2 

5 

G201.5  +0.6 

X 

HC 

201.5 

+0.6 

5 

VY  Coni* 

Major!* 

u 

X 

X 

o 

"b 

<20  fu 

240 

-6 

6,21 

RCW38 

A,U? 

A,U? 

A 

1  ? 

0 

H  II 

267.9 

-1.1 

0.5 

13 

RCW36 

X,A 

1  ? 
a 

H  II? 

265.  1 

♦  1.5 

1.5, 

3.5 

13 

IR0942  -  21 

u 

o? 

O? 

O 

"b 

<20  fu 

255.8 

+23.35 

21 

RCW46 

A,X 

1  ? 
o 

H  II? 

282.0 

-1.2 

4? 

13 

G284.2  -  0  8 
(ICW48-RCW49) 

u 

u,c 

U,A 

O 

1  ? 
o 

H  11 

284.2 

-0.8 

3.  8, 

4.9? 

14 

1308 -62 S 
(RCW74) 

E,C 

A 

1  ? 
o 

H  II 

305.3 

+0.2 

3,  8.5 

13 

1308-62N 

E 

A 

1  ? 
o 

H  II 

305.3 

♦0.2 

3,  8.5 

13 

1404-61 

X 

1  ? 
a 

311.9 

♦0. 1 

12.5 

13 

CIT7 

1 R 1 525  +  19G 

u 

u 

U 

O 

"b 

<20  fu 

29.5 

+53.5 

20,21 

1541-53 

A,E 

i  ? 
a 

326.7 

+0.6 

15.5 

(or  1) 

13 

G4.2  +35.8 

L134 

X 

HC 

4.2 

+  35-8 

5.12 

1608-51 

(MHR47) 

? 

A,C 

A 

1  ? 
a 

331.5 

-0.1 

11  (or  6.5) 

13 

1617-50N 

C ,  A 

A 

1 

o 

332.2 

-0.4 

3  to  4 

13 

1617-50S 

A 

C.A 

'b? 

332.2 

-0.4 

3  to  4 

13 

1618-49 

A,X 

i  ? 
a 

333.6 

-0.2 

2  to 

2.5 

13 

G352.9  +  16.  7 

X 

HC 

352.9 

♦  16.7 

5 

G353.4  +  16.9 

Cloud  4 

X 

HC 

353.4 

♦  16.9 

5,12 

NGC  6334S 

? 

E,  T,A 

E.T.A 

A(X?) 

'b 

H  II 

351.2 

+0.7 

0.7  to 

1.3  (?) 

1,2,8,11,13,19 

*  $•«  note*  at  *nd  of  tobl*. 
t  S««  r*f«r*snce»  at  «nd  of  toblv. 
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TABLE  XVII  (Continued) 

Transition  (MHz) 

Nearest 

18-cm 

b" 

Source  Nomt(i) 

1612 

1665 

1667 

1720 

Class 

Continuum 

(deg) 

(d«S) 

Distonce 

References  1 

NGC  6334  N 

? 

E,A 

C ,  A 

X,  A 

1 

a 

H  11 

351. 

+0.5 

0.7  to 
1.2(7) 

1,2,8,11,13,19 

NGC  6357 

W22 

A 

A,X 

A 

A 

1  ? 
a 

H  II 

353.2 

+0,67 

0.6  or 

1.3 

9,13 

IR 1 73-4  -32 
(NGC  6383?) 

U 

O? 

O? 

O 

"b? 

H  II? 

356.  4 

-0.3 

21 

Sog  A-| 

o? 

X,  A 

A? 

o? 

i 

a 

NTS? 

359.95 

-0.  15 

10? 

19 

Sog  A-ll 

A 

A 

A 

A,y 

II 

NTS? 

0 

0 

10? 

10 

Sgr.  62 
(W24-B2) 

U  or 

C,A 

C ,  A 

C,A 

A 

'b 

H  11? 

0.7 

-0.0 

10? 

17,19 

W28S 

G5. 9-0.4 

W2SA2 

u 

c,u 

C,U 

A? 

"b 

H  II 

5. 94 

-0.40 

3.0 

4,18 

W28N 

A 

A 

A 

u 

11 

a 

SNR  +  H  II 

6.  76 

-0.27 

1.3  to 

3.6 

1,3,4,9,10 

W33 

IC4701 

O 

E,T,A 

C,A 

A 

1 

a 

H  II  +  SNR 

12.7 

-0.0 

4.  7  or 

14.  8 

4,9,18 

W41 

X 

A 

A 

? 

"b 

SNR  +  H  II 
(?) 

23.2 

-0.3 

5  8  or 

12.6 

9,18 

G23.5  +  8.2 

X 

HC 

23.5 

+  8.2 

5 

W42 

X? 

A 

X 

0 

‘b 

NTS? 

24.6 

-0.1 

9,18 

G27.3  +4.3 

X 

HC 

27.3 

+4.3 

5 

G27.0  +3.5 

X 

HC 

27.0 

+3.5 

5 

G27. 1+3.5 

X 

HC 

27.1 

+  3.5 

5 

W43I 

A 

A 

A 

X,A 

II 

30.8 

-0.  1 

0.6  to 

1.0 

9,Ns 

W43II 

X,A 

E,A 

E,A 

X,  A 

‘b 

H  II 

30.8 

-0. 1 

6. 5  to 

7.5 

4,9,18 

W43I1I 

u 

X 

X 

O 

"b 

H  II 

30.8 

-0.1 

14? 

4,9,18 

W44I 

A 

A 

A 

U 

IIq 

SNR  +  H  II 

34.6 

-0.4 

1.6  to 

3.4 

1,3,4,10,18 

W44II 

X 

A 

A 

A 

"b 

? 

34.6 

-0.4 

1 

9,10,18 

IR  1903  +08 

u 

? 

? 

O 

"b 

<20  fu 

41.95 

+0.45 

21 

W49S 

c? 

C 

C? 

C,T?,A 

1 

a 

H  II 

43.2 

-0.0 

14 

1,2,13,17,19,22 

W49N 

C? 

c 

c? 

C,T?,A 

1 

a 

H  II 

43.2 

-0.0 

14 

1,2,13,17,19,22 

W51-I 

C,A 

u,c,a 

A(X?) 

C(E?) 

1 

a 

H  II  (NTS?) 

49.5 

-0.4 

6.5 

1,2,9,19 

W51-II 

A 

A 

A 

X 

11 

a 

H  II? 

49.5 

-0.4 

0.43 

1,2,9, 19 

ON3 

W58,  CTB82 

NGC  6857 

O 

c 

C 

c 

II 

a 

H  II 

70.37 

♦  1.60 

8.2 

7 

ONI 

o 

c 

C(?) 

o 

1 

a 

None 

69. 58 

♦  0.96 

1.8  or 

5.2 

7 

ON2 

W64,  CTB90 

o 

c 

c 

o 

1q 

H  II? 

75. 75 

+  0.35 

5.7 

7 

ON4 

u 

c? 

o 

o 

"b 

H  11 

77.  93 

+0.  17 

8. 1 

7 

W75N 

o 

c 

c? 

o 

i 

H  11  (or  none) 

81.86 

♦  0.76 

1.5? 

1,19,22 

W75S 

A 

E,T 

A 

A? 

i 

H  11 

81.75 

+0.54 

1.5? 

(or  6?) 

1,2,17,19. 22 

NML  Cygnus 

u 

C 

c 

0 

"b 

<  5  fu 

80.8 

-1.9 

?<300  pc? 

21,22 

Cat  A 

A 

A 

A 

A,l 

ii 

111.7 

-2.  1 

9,18 

G1H.5  +  14.6 

Cloud  IB 

X 

HC 

114.5 

+  14.6 

5,12 

‘Notes 

A  Abtorptian 

C  Circularly  polarized  emission 

L  Linearly  polarized  emission 

E  Bath  circular  and  linear  emission  present 

U  Subitantiolly  unpolarized  (<20  percent)  emission 

X  Emission,  but  polarization  unknown 
?  Questionable 

T  Time  variable 

0  No  emission  or  absorption  detected 

Underline  indicates  transition  of  strongest  em 

iuion. 
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Fig.  52.  Location  of  OH  emission  sources  In  the  galactic  plane. 
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regions  has  a  higher  excitation  temperature  compared  with  OH  in  the  general  interstellar  me¬ 
dium,  and  perhaps  a  higher  relative  density  also.  The  polarization  of  these  regions  has  not  been 
measured  and  the  line  ratios  are  known  only  very  crudely.  If  the  OH  in  Heiles'  clouds  does  turn 
out  to  be  in  LTE,  it  will  be  almost  the  only  extraterrestrial  OH  known  to  be  in  such  a  state. 

The  cases  of  anomalous  absorption  [Goss  (1968)]  which  maybe  combinations  of  absorption 
and  anomalous  emission  are  also  not  included  and  not  listed  in  Table  XVII.  Goss  states  that 
there  is  no  case  in  his  survey  for  which  the  OH  absorption  can  be  characterized  by  LTE  (i.e.,  a 
single  excitation  temperature  for  the  four  18-cm  OH  transitions).  It  is  clear  that  some  mech¬ 
anism  must  be  effective  in  the  general  interstellar  medium  which  produces  these  unequal  exci¬ 
tation  temperatures  in  OH.  It  is  not  clear  what  this*  mechanism  is  or  whether  the  same  mech¬ 
anism  produces  the  population  inversion  in  some  or  all  of  the  anomalous  OH  emission  sources. 

Figure  52  shows  a  projection  onto  the  galactic  plane  of  all  the  OH  emission  sources  (from 
Table  XVII)  that  have  some  estimate  for  their  distance.  No  obvious  patterns  can  be  seen  in  this 
plot,  in  particular: 

(1)  The  more  distant  (and  intrinsically  stronger)  sources  do  not  all  belong 
to  a  given  class. 

(2)  The  closer  sources  are  not  more  strongly  polarized,  although  there  may 
perhaps  be  some  tendency  for  distant  sources  to  have  a  lower  percentage 
of  linear  polarization. 

(3)  There  are  at  least  a  few  sources  which  do  not  lie  in  spiral  arms  as  traced 
by  neutral  hydrogen  emission. 

These  OH  emission  sources  are  quite  strongly  concentrated  toward  the  galactic  plane;  the 
largest  well-established  Z  distance  is  2 30 pcs  for  W58  (ON- 3).  And  no  sources  are  known  to  be 
farther  than  about  1  3  kpc  from  the  galactic  center. 

B.  Action  of  Ions  on  OH  Molecules 

1.  Introduction 

Because  of  the  apparent  phenomenological  association  between  OH  emission  sources  and 
SNRs,  it  is  desirable  to  reexamine  the  interaction  between  ions  and  OH  molecules  to  determine 
whether  this  interaction  can  provide  either  a  pumping  mechanism,  or  a  spoiler  mechanism  to 
explain  the  peculiarities  of  those  OH  sources  (Class  11^)  most  closely  associated  with  SNRs  The 
problem  of  the  interaction  of  ions  and  OH  molecules  has  been  considered  by  Johnston  (1967)  and 
Turner  (1967)  using  an  approach  that  made  use  of  the  Born  Approximation,  and  by  Rogers  (1964), 
Rogers  and  Barrett  (1968),  and  Goss  and  Field  (1968)  using  time-dependent  perturbation  theory 
based  on  an  approach  by  Purcell  (1952).  Using  the  latter  approach  this  study  extends  the  previ¬ 
ous  work  by  (a)  employing  the  relativistically  correct  form  of  the  equations  and  taking  relativ¬ 
istic  ions  into  consideration,  (b)  dropping  some  of  the  approximations  made  in  previous  work, 
and  (c)  using  instead  of  the  Maxwell  equation,  a  form  for  the  ion  energy  density  which  is  thought 
to  be  appropriate  for  the  vicinity  of  a  supernova  remnant. 

2.  Order-of- Magnitude  of  the  Effect 

It  is  perhaps  best  to  begin  by  doing  an  order-of-magnitude  calculation  to  indicate  that  it  is 
at  least  conceivable  that  the  charged  particle  —  OH  molecule  interaction  in  the  vicinity  of  an  SNR 
can  significantly  affect  the  observed  OH  radio  emission. 
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We  choose  the  1720-MHz  OH  emission  from  W28  because  it  is  well  known  (but  probably 

stronger  than  the  average  type  II  OH  emission  source)  and  we  calculate  the  spectral  flux  for 

a  _2i  2 

the  two  strong  features  as  9.55  x  10  watts/m  [Ball  and  Staelin  (1968)].  Using  the  distance 

of  4200 pcs  suggested  by  Goss  (1968),  which  is  much  larger  than  the  1300 pcs  suggested  by  Poveda 

2 1 

and  Woltjer  (1968),  we  calculate  a  spectral  luminosity  of  1.9  x  10  watts  and  this  corresponds 
45  / 

to  about  1.7  x  10  photons/sec.  This  calculation  assumes  that  the  OH  emission  is  isotropic  or 
at  least  that  the  spectral  intensity  in  the  direction  of  the  earth  is  typical  of  the  source.  If  we  had 
used  1300  pcs  for  the  distance,  the  inferred  spectral  luminosity  would  be  smaller  by  about  an 
order  of  magnitude. 

We  have  only  a  rather  large  upper  limit  on  the  angular  size  of  the  OH  emission  source  in 
W28  so  we  assume  a  diameter  of  4000 au  which  is  approximately  the  value  obtained  by  Moran, 

et  al.  (1968)  for  the  total  size  of  the  emitting  region  in  W3.  From  this  we  calculate  a  volume 

50  3 

of  9  x  10  cm  . 

~2  3  3 

Inside  this  volume  we  assume  there  are  10  OH  molecules/cm  and  500  electrons/cm  with 
speeds  of  10  to  100  km/sec.  This  OH  density  can  only  be  a  guess,  but  it  is  clear  that  something 
is  special  about  the  OH  emitting  region  to  distinguish  it  from  other  nearby  regions  that  do  not 
have  OH  emission.  This  high  an  OH  density  was  used  by  other  workers  who  calculated  maser 
models.  The  electron  density  of  500  electrons/cm  is  lower  than  the  maximum  estimate  of 
Minkowski  (1968)  and  others,  but  the  speed  is  also  slower  than  the  speed  of  mass  motions  ob- 
served  in  SNRs.  We  now  read  off  a  transition  rate  of  10  cm  /sec  from  Fig.  4  of  Rogers  and 
Barrett  (1968)  and  with  the  numbers  above  we  calculate  4.5  x  104^  transitions/sec  in  the  region. 
This  is  about  an  order  of  magnitude  larger  than  the  number  of  photons/sec  calculated  above. 

We  conclude  that  if  the  conditions  assumed  above  are  reasonable,  then  the  interaction  be¬ 
tween  OH  molecules  and  passing  ions  will  certainly  be  important  and  may  in  fact  be  sufficient 
to  explain  the  source  of  energy  for  the  OH  emission. 

For  the  case  in  which  the  ions  are  the  sole  energy  source,  another  point  is  worth  noting. 
Because  the  mechanism  under  consideration  is  not  a  maser,  i.e.,  no  inversion  of  the  population 
levels  is  expected,  the  brightness  temperature  of  the  18-cm  radiation  cannot  exceed  the  equiv¬ 
alent  kinetic  temperature  of  the  ions.  For  the  case  in  which  the  ions  have  a  non- Maxwellian  ve¬ 
locity  distribution,  we  must  at  least  have  a  substantial  number  of  ions  with  energies  comparable 
to  k  times  the  brightness  temperature.  To  see  whether  this  is  the  case,  we  first  calculate  the 

7 

peak  brightness  temperature  for  Feature  1  of  W28  at  1720  MHz  as  2.5  x  10  °K  using  the  size  as¬ 
sumed  above.  We  have  no  way  of  knowing  whether  this  is  reasonable  for  a  type  II  source,  but 
the  type  I  source  in  W3  has  much  higher  brightness  temperatures  [Moran,  et  al.  (1968)].  To 
correspond  to  this  kinetic  temperature,  an  electron  must  move  at  about  3  x  10^  cm/sec  and  a 
proton  about  a  fortieth  of  this.  This  electron  speed  is  nearly  relativistic,  which  suggests  the 
need  for  a  relativistic  analysis,  and  is  much  higher  than  we  assumed  previously,  which  illus¬ 
trates  that  we  do  not  have  the  whole  story. 

3.  Equivalent  Electric  Field  Energy  Density  and  the  Transition  Probability 

In  this  section  we  rederive  the  appropriate  equations  of  Rogers  (1964),  Rogers  and  Barrett 
(1968)  and  Goss  and  Field  (1968)  in  their  relativistically  correct  form. 

Consider  an  ion  with  charge  q  and  velocity  v  passing  an  OH  molecule  with  impact  param¬ 
eter  b.  Figure  53  shows  the  situation.  We  assume  that  the  ion  is  not  appreciably  deflected  by 
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the  OH  molecule.  According  to  Jackson  (1962),  the  nonzero  components  of  E  and  B  are 


E ,  = 


(b2  + 


yqb 

2  2.2. 3/2 
y  v  t  ) 


(46) 


E3-- 


_ iaxt _ 

,.2  7  2  2, 2, 3/2 

(b  +  y  v  t  ) 


(47) 


B2  =  /3E1  (48) 

2  _  i  /2 

where  /?  =  v/c  and  y  =  (1  —  (3  )  /  .  Note  that  in  the  relativistic  limit  (v  -*  c)  these  equations  are 

the  same  as  the  equations  describing  a  pulse  of  plane  polarized  electromagnetic  radiation  prop¬ 
agating  in  the  direction. 


Fig.  53.  Ion  passing  an  OH  molecule. 


To  use  the  approach  suggested  by  Purcell  (1952)  we  need  the  Fourier  spectrum  of  the  E  field. 
If  we  define 


oo 

F(u>)  E  \  F(t)  elaJt  dt 
°LoO 


(49) 


then  we  get 


K.(2^) 

1  yv 


E3(w) 


^4  K  (^) 
2  2  o  yv 
y  v 


(50) 

(51) 


(see  Jackson  (1962)].  The  K  functions  are  the  modified  Bessel  functions,  the  properties  of 
which  are  summarized  in  Abramowitz  and  Stegun  (1964).  Subsequent  formulae  involving  these 
functions  are  either  taken  from  Abramowitz  and  Stegun  or  are  derivable  from  equations  given 
by  them. 

The  energy  density  and  also  the  interaction  energy  with  a  harmonic  oscillator  are  propor¬ 
tional  to  |  E  (a? )  |  ^  which  is 


E(u>) 


.  2  2 
2  4q  a) 

=  2  4 

y  v 


K,2(— )  +  A  K  2 (rr~) 
1  'yv  2  o  yv 

y 


(52) 
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Several  properties  of  this  equation  are  worth  noting 

-*■7  1  /2 

(a)  |E(u>)|  falls  exponentially  at  large  b  [because  Kq(x)  —  K^(x)  -*  (tt/2x) 

exp  {— x)  for  x  —  °°]. 

(b)  Although  the  E^  spectrum  has  a  peak  at  **  1.35  yv/b  [because  xKq(x)  = 

maximum  at  x  ^  1.35]  there  is  no  peak  for  the  total  |E(a>)|2  [because 

xK.  (x)  =  maximum  at  x  =  0]. 
l 

(c)  Although  this  equation  has  a  singularity  for  b  —  0,  it  is  well  defined  for 
Cl?  -►  0  and  also  for  y  —  °°. 

Using  Eq.  (52)  we  may  write  the  relativistically  correct  equivalent  of  Eq.  (18)  in  Rogers  and 
Barrett  (1968)  as 


P. 

r 


♦  (b )  = 


.  2  2 
4q  cj 


U  . 

3  fi  y  v 


L2(^)  +  4  k2(^)1 
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(53) 


where  p. .  is  the  electric  dipole  moment  connecting  the  two  states  in  question,  and  P.  .  is  the 
ij  J 

probability  of  producing  a  transition. 

In  the  low  frequency  or  high  energy  limit  u>  «  yv/b  we  get  the  approximation 


|  E  (oj  )  | 


v  b 


1  ru>b  ,  ,u>bn2 
1  +  — ~  ln( 

2  lyv  yv 
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vb4 


(54) 


[because  xk^(x)  —  1,  xKq(x)  —  0  as  x  —  0).  Note  that: 


(a)  The  spectrum  becomes  flat  (the  u>  dependence  cancels  out)  correspond¬ 
ing  to  the  pulse  character  of  E(t). 

^  2  o 

(b)  |  E  (co )  |  falls  as  l/b  rather  than  exponentially  as  above. 

In  the  high  frequency  or  low  speed  limit  a)  »  yv/b  we  get  the  approximation 


(55) 


This  spectrum  is  illustrated  schematically  in  Fig.  54. 

For  the  conditions  under  which  Eq.  (54)  holds  (o>  «  yv/b)  we  have  the  approximation 


Pj~i(b)  “  ^ 


2  .  2 
4q 


„  ITT 

3  b  vb 


(56) 


Fig.  54.  Equivalent  energy  spectrum 
of  a  passing  ion. 


80 


provided  Pj_^.(b)  <  |  Following  Rogers  and  Barrett  (1968)  we  define  r^  by  P.^.Jr^)  =  i  in  Eq.  (56) 
and  calculate 


r 


2 

1 


8q 


ni 

ITT 


3ft  v‘ 


(57) 


(The  equivalent  of  Eq.  (57)  in  Rogers  and  Barrett  contains  a  misprint;  it  should  have  an  exponent 
of  }.)  The  use  of  Eq.  (56)  to  calculate  would  be  incorrect  under  circumstances  in  which 
P.  ~  }  but  co  >  yv/r^  We  shall  see  below  that  this  is  unlikely. 

In  summary,  provided  co  «  yv/r^,  we  use  P._^.(b)  =  \  for  b  <  r^  and  Eq.  (56)  for  b  >  r^  un¬ 
less  co  >  yv/b  in  which  case  we  must  use  Eq.  (53).  If  co  >  yv/r  we  must  use  Eq.  (53)  both  for 
calculating  r ^  and  for  Pj_^(b)  for  b  >  r^. 

We  now  assume  that  the  OH  molecule  is  in  a  stream  of  mono-energetic  ions  characterized 
by  Nj,  their  number  density  per  unit  volume.  Then  the  number  of  ions  passing  the  OH  molecule 
per  unit  time  with  impact  parameters  in  the  range  b  to  b  +  db  will  be 

dN4  =  N 1  v2 Trb  db  (58) 


We  wish  to  evaluate 

H^v) 


\  P.  .(b)  N  .  v2?rb  db  =  2ttINLv 

J  J-i  1  1 


P.  .(b)  b  db 
J-i 


(59) 


which  is  the  rate  of  j  —  i  transitions  caused  by  the  ions.  Using  Eq.  (53)  [as  we  must  because 
the  approximation  used  to  get  Eq.  (56)  breaks  down  somewhere  in  the  range  of  integration],  the 
integral  in  Eq.  (59)  is 
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These  integrals  can  be  done; 


(  xK2(x)  dx  =  [Kq2(x)-  K2(x)] 


(  xK  2(x)  dx  -  [K2(x)  -  K2(x)]  -  xKq(x)  K^x) 

so  Eq.  (60)  becomes,  after  some  manipulation, 
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However,  we  have  already  assumed  that  cj  «  yv/r^,  so  Eq.  (63)  reduces  to 


P.  .(b)  b  db 


2 »  2 

rm. 4q 


(64) 


[because  Kq(x)  ■*  — ln(x)  for  x  -»  0].  Now  we  may  write  Eq.  (59)  as 
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This  equation  is  comparable  with  Eq.  (24)  of  Rogers  and  Barrett  (1968)  and  Eq.  (12)  of  Purcell 
(1952),  but  Eq.  (65)  is  for  mono-energetic  ions  and  so  differs  in  detail. 

We  now  assume  that  the  OH  molecule  is  in  a  stream  of  ions  characterized  by  an  energy  spec¬ 
trum  such  as  that  shown  in  Fig.  55(a),  i.e.,  a  power  law  for  energies  above  a  limiting  energy  Eq 

and  a  constant  below  E  .  We  assume  that  E  «  m  C2,  i.e.,  that  the  cutoff  energy  E  is  well 
o  o  o  o 

below  the  relativistic  range,  and  we  assume  1  <  ft  “  1.5  or  2. 

The  form  of  this  energy  spectrum  and  the  assumption  about  Eq  perhaps  deserve  some  com¬ 
ment.  The  power  law  energy  spectrum,  and  values  for  a  and  K  but  not  E  ,  are  determined  for 
electrons  in  supernova  remnants  from  radio  observations  of  the  synchrotron  continuum  radiation 
and  some  assumption  about  the  magnetic  field.  The  energy  spectrum  is  therefore  known  only  for 
electrons  in  a  limited  range  of  relativistic  energies.  It  is  also  not  clear  just  what  part  of  a  su¬ 
pernova  remnant  is  being  observed  in  the  radio  continuum.  It  is  likely  that  conditions  near  an 
OH  emission  source  on  the  edge  of  a  supernova  remnant  are  very  different.  We  shall  see  be¬ 
low  that  ions  with  low  energies  are  very  important,  and  the  energy  spectrum  in  Fig.  55(a) 
should  be  considered  only  as  plausible  but  certainly  not  uniquely  specified  by  any  observations. 


(a) 


(b) 


Fig.  55.  Schemotic  illustration  of  the  cross-section  integral. 
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Furthermore,  we  need  to  determine  the  third  parameter  Eq.  We  do  this  by  specifying  that  the 
total  ion  density  —  the  integral  under  the  curve  in  Fig.  55(a)  —  be  some  reasonable  but  assumed 
value. 

The  total  ion  density  jfrom  Fig.  55(a)]  is 


H-r 


N  (E )  dE  =  KE 


1  -  a 


— — — i 
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(66) 


We  now  wish  to  perform  the  total  rate  integral  which  is 

r4[v(e)] 

n  =  \  — ^ -  N(E)  dE  (67) 


This  will  give  us  the  number  of  transitions  per  OH  molecule  per  second  for  an  OH  molecule  in 
the  ion  field  corresponding  to  N(E).  It  will  be  necessary  to  break  up  the  range  of  the  integral 
in  Eq  (67)  into  three  pieces: 


(1)  E  <*  E  for  which  N(E)  =  KE  a  =  constant  and  for  which  E  =  l/2m  v^; 

o  o  o’ 

2  _  Q, 

(2)  E  <  E  <  l/2m^c  for  which  N(E)  =  KE  and  for  which  we  assume 

E  »  i/2m  v^,  and 
o 

2  -a 

(3)  E  >  l/2mQC  for  which  N(E)  =  KE  and  for  which  we  assume  v  -  c 
(but  use  the  correct  value  for  y,  of  course). 

The  results  of  these  three  integrals  are 
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2 

For  any  reasonable  value  of  N  (several  hundred/cm  )  and  with  values  of  K  and  ot  usually 
2 

used,  we  get  Eq  «  mQc  and  this  immediately  shows  that  and  the  first  term  of  are  negli¬ 
gible  with  respect  to  and  the  second  term  of  f 2^,  which  are  comparable.  This  result  means 
that  the  low  energy  electrons  contribute  almost  all  of  the  interaction.  And  the  form  for  N(E)  is 
probably  not  correct  at  low  energies.  The  sought-for  effect  (coupling  of  the  OH  microwave  tran¬ 
sitions  to  passing  high  energy  ions)  exists  but  it  is  completely  dominated  by  the  interaction  with 
low  energy  ions  for  the  energy  spectrum  assumed.  The  interaction  between  OH  molecules  and 
passing  ions  can  explain  neither  the  1720- MHz  emission  in  type  Ha  sources  nor  the  nonequilib¬ 
rium  excitation  temperature  in  the  OH  clouds  observed  by  Goss  (1968). 

In  conclusion  it  seems  likely  that  the  apparent  association  between  OH  emission  sources  and 
SNRs  is  an  association  in  the  sense  that  certain  regions  of  the  galaxy,  perhaps  because  of  the 
chemical  composition  of  the  interstellar  medium,  the  presence  of  large  amounts  of  dust,  or  other 
causes,  tend  to  form  stars  with  the  kind  of  circumstellar  clouds  that  show  OH  emission  and  also 
tend  to  form  the  kind  of  stars  that  explode  as  supernovae.  Thus  the  two  kinds  of  objects  are 
proximate  in  space  but  probably  not  interacting.  Very  likely  the  association  between  OH  emis¬ 
sion  sources  and  H  II  regions  is  of  the  same  character.  The  association  between  OH  emission 
sources  and  IR  stars  will  probably  turn  out  to  be  a  more  direct  causal  relationship. 

C.  Physical  Pictures  of  OH  Emission  Sources 

It  was  pointed  out  above  that  most  OH  emission  features  that  contain  linear  polarization  and 
also  some  that  contain  only  circular  polarization  show  a  variation  of  the  polarization  properties 
(the  percent  circular  and  the  percent  linear  polarization  and  the  linear  polarization  position  angle) 
as  a  function  of  velocity  across  the  feature.  There  are  several  possible  causes  for  this  behavior- 

(1)  The  feature  itself  may  have,  intrinsically,  polarization  properties  that 
vary  as  a  function  of  frequency  or  velocity.  Until  the  mechanism  which 
produces  the  polarization  is  understood,  this  remains  a  possibility  at 
least  for  some  features  in  some  sources. 

(2)  If  we  say,  rather  arbitrarily,  that  an  individual  feature  is  to  be  Gaussian 
or  at  least  symmetrical  in  shape,  then  most  spectra  are  made  up  of  com¬ 
posite  features  and  the  variation  of  the  polarization  properties  may  be 
due  to  the  incoherent  overlap  of  separate  features,  each  of  which  has  its 
own  polarization. 

(3)  Faraday  rotation  in  the  intervening  medium  rotates  the  plane  of  linear 
polarization  as  a  function  of  frequency.  However,  this  is  unlikely  to  be 
the  cause  of  the  rotation  of  the  linear  polarization  position  angle  across 
the  spectrum  if  the  Faraday  rotation  is  to  take  place  either  in  the  inter¬ 
stellar  medium  or  in  the  earth's  ionosphere,  because  the  rotation  meas¬ 
ures  are  far  too  small.  For  example,  to  Taurus  A  the  total  rotation 
from  X  =  0  to  X  =  20  cm  is  only  about  60°  [see  Kraus  (1966),  p.  146].  It 
is  conceivable  that  the  OH  emission  originates  inside  a  region  of  much 
higher  rotation  measure. 

Regardless  of  which  of  these  possibilities  turns  out  to  be  correct,  it  is  likely  that  some  of 
the  features  in  some  of  the  sources  (e.g.,  W49  and  W51)  will  turn  out  to  be  more  highly  polarized 
when  examined  with  higher  frequency  and  spatial  resolution. 

Moran,  et  al.  (1968)  show  a  map  of  the  1665-MHz  OH  emission  from  the  W3  region.  They 
find  that  each  feature  comes  from  a  spot  or  small  region  0.02  arcsecond  or  less  in  size,  but  the 
separation  among  these  spots  is  typically  0.2  to  2  arcseconds.  At  the  distance  of  W3  (which  we 
have  taken  to  be  3kpc)  this  corresponds  to  a  projected  linear  size  of  50 au  or  less  for  the  size  of 
the  spots  and  a  few  hundred  to  several  thousand  au  for  their  separations. 
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It  is  not  known  whether  the  length  of  the  OH  emitting  region  along  the  line  of  sight  is  about 
the  same  as  the  spot  size  or  the  separation,  or  has  some  other  value.  For  some  models  of  the 
OH  emission  it  is  a  weak  general  rule  to  say  that  the  time  scale  of  any  time  variations  is  equal 
to  the  light  travel  time  through  the  emitting  region  along  the  line  of  sight.  This  is  true,  for  ex¬ 
ample,  of  a  saturated  maser  that  has  its  pump  approximately  at  right  angles  to  the  line  of  sight. 
The  output  of  such  a  maser  is  insensitive  to  changes  in  the  microwave  signal  input,  but  a  change 
in  the  pump  power  will  change  the  output.  The  change  in  output  of  various  regions  along  the  line 
of  sight  will  not  reach  the  observer  at  the  same  time,  however,  and  the  result  will  be  a  rise  or 
decay  no  faster  than  the  time  scale  mentioned  above. 

Although  the  OH  emission  from  the  W3  region  has  been  reported  to  be  time  varying  [Coles, 
et  al .  (1968),  Menon  (1967b)  and  Sec.  V-A  of  this  report],  this  has  been  only  partially  confirmed. 
The  OH  emission  from  NGC  6334S,  a  source  in  the  same  class  as  W3,  has  varied  dramatically 
in  a  time  of  9  days  or  less  (see  Sec.  V-E).  Using  the  rule  above,  we  get  a  length  along  the  line 
of  sight  of  about  1600  au.  Other  time  varying  OH  emission  sources  seem  to  vary  more  slowly. 
There  is  no  conflict  with  the  assumption  that  the  lengths  of  the  OH  emitting  regions  are  much 
larger  than  their  apparent  projected  diameters.  In  any  case,  some  models  for  the  OH  emission 
sources  have  no  restriction  on  the  time  scale  of  the  observed  time  variations.  For  example, 
almost  arbitrarily  rapid  variations  may  occur  with  a  model  that  has  the  OH  emission  confined  to 
a  narrow  angular  beam  which  rotates  onto  and  off  of  the  line  of  sight  to  the  observer. 

This  illustrates  the  value  of  searching  for  time  variations  on  a  time  scale  of  hours.  Such 
variations  have  been  reported  by  Menon  (1967);  however,  they  have  not  been  confirmed  by  other 
observers  and  in  any  case  were  attributed  to  interstellar  scintillations  rather  than  to  changes  in 
the  sources.  Robinson  and  Goss  (1968)  have  reported  variations  with  a  time  scale  of  millisec¬ 
onds  due  to  interplanetary  scintillations. 

It  now  seems  likely  that  some  or  all  of  the  OH  emission  sources  are  circumstellar  clouds 
around  a  single  star  or  at  least  objects  whose  characteristic  scale  size  is  appropriate  for  solar 
system  dimensions  rather  than  typical  interstellar  distances.  Still  open  is  the  question  of  whether 
these  are  young  objects  in  which  the  circumstellar  cloud  is  in  the  process  of  contracting  to  form 
a  star,  or  old  objects  in  which  the  cloud  has  been  ejected  from  the  star.  If  it  can  be  shown  that 
the  OH  velocity  field  corre sponds  to  contraction  or  rotation,  we  would  conclude  that  the  object 
is  young;  expansion  would  correspond  to  an  old  object.  Furthermore,  if  the  space  velocities  of 
the  OH  turn  out  to  be  larger  than  the  escape  velocity  for  some  measured  or  assumed  mass,  then 
we  would  also  conclude  that  the  object  is  old.  The  mass  of  an  object  must  be  28  solar  masses 
for  the  escape  velocity  from  a  distance  of  2000  au  to  be  5  km/sec.  These  numbers  are  appro¬ 
priate  for  W3  [Moran,  et  al.  (1968)].  Arguments  presented  below  lead  us  to  suspect  that 
5  km/sec  is  less  than  the  true  space  velocity  for  the  OH  in  objects  such  as  W3  and  so  this  mass 
is  a  lower  limit  if  the  OH  is  gravitationally  bound  corresponding  to  a  young  object.  We  have  no 
way  of  knowing  whether  28  solar  masses  is  a  reasonable  number. 

It  is  interesting  in  this  connection  to  examine  the  spread  in  Doppler  velocities  among  features 
from  each  source  and  compare  the  median  or  centroid  of  these  velocities,  which  probably  should 
be  interpreted  as  the  Doppler  velocity  of  the  cloud  as  a  whole,  with  the  Doppler  velocities  of  the 
objects  with  which  the  OH  emission  is  presumably  associated  (usually  H  II  regions).  In  all  cases 
the  OH  emission  velocity  spread  includes  the  velocity  of  the  associated  object.  Furthermore, 
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the  offset  of  the  velocity  of  the  associated  object  from  the  median  of  OH  emission  velocities  is 
often  much  less  than  the  spread  of  the  OH  emission  velocities.  This  suggests  that  the  cloud  of 
OH  is  moving  rather  rapidly  (up  to  several  tens  of  km/sec)  with  respect  to  the  star  it  surrounds 
(expanding,  contracting,  or  possibly  rotating)  but  that  the  star  and  its  envelope  as  a  whole  are 
moving  rather  slowly  (typically  less  than  10  km/sec)  with  respect  to  nearby  objects  [see  also 
Menon  (1967a)].  These  latter  velocities  are  typically  less  than  the  random  velocities  of  individ¬ 
ual  stars  and  this  may  be  interpreted  to  mean  that  the  OH  emission  sources  are  more  massive 
than  most  stars  and  so  their  velocities  on  an  equipartition  basis  are  smaller.  Menon  (1967)  has 
pointed  out  that  the  Doppler  velocity  of  the  OH  emission  from  Orion  is  closer  to  the  Doppler  ve¬ 
locity  of  the  cluster  stars  than  to  the  Doppler  velpcity  of  the  H  II  region,  and  he  interprets  this 
to  mean  that  the  OH  emission  source  is  a  member  of  the  same  class  of  objects  as  the  cluster  stars. 

With  the  discovery  of  more  OH  emission  sources  associated  with  visible  objects  such  as 
VY  Canis  Majoris,  it  should  be  possible  to  explore  these  correlations  in  detail. 

With  the  trepidation  appropriate  to  a  considerable  extrapolation  of  our  present  knowledge, 
the  physical  picture  in  Fig.  56  is  proposed  to  explain  some  of  the  peculiarities  of  the  OH  emission 

sources.  Figure  56  shows  an  expanding  enve¬ 
lope  of  OH  surrounding  a  star  which  provides 
the  pump  energy  to  invert  the  OH  and  also  pro¬ 
vides  one  source  of  microwave  signal  input  to 
the  resulting  maser.  The  OH  in  the  inner  thick 
shell  either  is  not  saturated  or  is  saturated  by 
the  microwave  signal  from  the  star.  Thus  only 
the  cross-hatched  region  is  visible  to  the  ob¬ 
server.  If  the  shell  is  expanding,  the  Doppler 
velocity  of  this  visible  OH  will  be  negative  with 
respect  to  the  Doppler  velocity  of  the  star  itself. 
This  thick  shell  corresponds,  for  example,  to 
the  single  feature  at  1612  MHz  seen  in  W28S. 

If  it  is  more  strongly  inverted,  the  outer 
thin  shell  may  be  saturated  by  OH  emission 
traveling  approximately  transverse  to  the  direc¬ 
tion  to  the  star.  This  is  because  of  the  geometry  which  allows  amplification  of  spontaneous 
emission  along  the  longer  transverse  path  until  it  dominates  over  microwave  radiation  from  the 
star.  Thus  only  the  cross-hatched  regions  of  the  thin  shell  will  be  observed  in  OH  emission  and 
this  emission  is  likely  to  be  saturated  and  so  polarized.  The  observed  Doppler  velocity  spread 
will  be  small  because  the  expansion  (or  contraction)  velocity  is  largely  transverse  to  the  line  of 
sight  to  the  observer.  Because  of  saturation  and  irregular  motions,  the  appearance  of  the  re¬ 
sulting  OH  emission  will  be  patchy  and  irregular  over  the  projected  ring.  The  thin  shell  model 
corresponds  to  type  I  sources  and  in  particular  to  the  interferometer  picture  of  the  OH  emission 
from  W3  by  Moran,  et  al .  (1968). 

Intermediate  cases  also  occur,  as,  for  example,  a  relatively  thin  shell  which  is  not  saturated 
and  so  can  be  seen  both  in  amplified  spontaneous  emission  and  amplified  stellar  microwave  con¬ 
tinuum.  The  Doppler  velocity  spread  in  this  case  includes  twice  the  expansion  velocity  (because 
even  the  far  side  of  the  shell  can  be  seen)  and  this  case  corresponds,  for  example,  to  the  very 
broad  and  relatively  unpolarized  emission  from  W28S  at  1665  and  1667  MHz. 


Fig.  56.  Schematic  model  of  an  OH  emission  source. 
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There  are  a  number  of  obvious  features  of  some  of  the  OH  emission  sources  which  cannot 
be  understood  in  terms  of  this  simple  picture.  One  of  the  more  disturbing  is  the  frequent  occur¬ 
rence  of  two  groups  of  features  at  well  separated  Doppler  velocities  (for  example,  NML  Cygnus 
and  VY  Canis  Majoris  at  1612  MHz)  which  may  correspond  to  the  front  and  rear  of  an  expanding 
shell  as  suggested  by  Wilson  and  Barrett  (1968b);  what  is  not  clear  is  why  the  rear  of  the  shell 
should  be  prominent.  It  may  be  that  the  velocity  field  is  such  that  the  path  length  for  amplifica¬ 
tion  of  a  given  frequency  is  longer  in  the  radial  direction  either  toward  or  away  from  the  star, 
than  in  any  transverse  direction.  This  will  be  the  case  if  there  is  little  dispersion  in  velocity 
radial  to  the  star,  as  might  occur  if  the  envelope  were  ejected  gradually  over  an  extended  time 
rather  than  by  a  single  explosion. 

From  this  picture  several  predictions  may  be  made  to  be  tested  by  future  observations. 
Among  these  predictions  are:  (1)  The  1612  MHz  feature  in  W28S  should  appear  to  come  from 
the  center  of  a  larger  region  defined  by  the  1665  and  1667  MHz  emission.  (2)  The  features  from 
type  I  sources  (with  narrow  velocity  spreads  and  highly  polarized  features)  should  define  a  rough 
ring  when  examined  with  high  angular  resolution.  (3)  For  sources  with  a  large  Doppler  velocity 
spread  and  relatively  unpolarized  features  (for  example,  W28S  at  1665  and  1667  MHz)  the  central 
features  in  velocity  should  define  a  region  on  the  sky  inside  which  should  lie  the  higher  velocity 
features  (if  the  cloud  is  indeed  expanding). 

These  predictions  must  be  tested  using  the  techniques  of  very  long  baseline  interferometry. 

VII.  CONCLUSIONS  AND  SUGGESTIONS  FOR  FUTURE  WORK 
A.  Conclusions 

1.  New  Observational  Constraints  on  OH  Source  Models 

A  considerable  amount  of  specific  new  information  about  OH  emission  sources  has  been  ob¬ 
tained  in  this  study.  Some  of  this  information  is  probably  related  only  to  an  individual  source  in 
the  sense  that  each  source  is  unique.  However,  a  number  of  definite  patterns  are  emerging  and 
it  is  these  patterns  that  must  first  be  explained  by  OH  source  models.  A  completely  satisfactory 
OH  source  model  must  now  explain  not  only  the  high  brightness  temperature,  the  small  angular 
sizes,  and  the  strong  linear  and  circular  polarization  of  the  OH  emission  sources,  but  also  the 
variation  of  the  percent  circular  and  linear  polarization  and  the  linear  polarization  position  angle 
with  velocity  across  spectral  features.  In  every  feature  that  contains  strong  linear  polarization, 
there  is  evidence  for  such  a  variation  of  the  polarization  position  angle. 

The  dramatic  and  unambiguous  time  variations  in  the  OH  emission  from  NGC  6334S  in  a  time 
of  9  days  or  less  must  also  be  explained  by  source  models.  Some  source  models  (in  general 
those  with  isotropic  OH  emission  output)  must  be  restricted  in  size  by  this  time  scale  to  9  light 
davs  (1600  au)  in  thickness  along  the  line  of  sight.  Alternatively,  the  time  variations  might  be 
explained  by  a  model  which  emits  OH  photons  only  into  a  small  solid  angle  which  may  then  rotate 
on  and  off  the  line  of  sight  to  the  observer.  A  lack  of  time  variation  in  the  polarization  properties 
of  the  variable  features  seems  to  be  indicated  by  the  data  and  must  be  interpreted  in  terms  of  a 
lack  of  angular  variation  of  the  polarization  properties  of  the  OH  emission  for  this  model.  This 
point,  however,  deserves  further  work 

Bv  contrast,  some  of  the  features  in  some  of  the  sources  have  remained  remarkably  constant 
in  time.  Most  of  the  features  in  W3,  for  example,  have  not  changed  in  any  way  after  a  time  of 
more  than  three  years. 
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Some  source  models  imply  the  presence  of  detectable  amounts  of  thermal  microwave  emis¬ 
sion  to  be  coincident  with  the  OH  emission  position.  This  idea  is  supported  by  the  continuum 
observations  of  W3  and  some  other  sources.  However,  W7  5N  and  other  sources  have  much  less 
continuum  emission  or  perhaps  none  at  all.  An  OH  source  model  which  predicts  a  large  amount 
of  microwave  continuum  must  be  incorrect  for  these  sources. 

2.  Associations  Between  OH  Emission  Sources  and  Nearby  Objects 

Observational  selection  effects  may  contribute  to  apparent  associations  between  OH  emission 
sources  and  nearby  objects.  Furthermore,  phenomenological  associations  do  not  necessarily 
imply  causal  relationships.  However,  some  associations  seem  to  be  well  established  and  per¬ 
haps  some  causal  relationships  implied  thereby. 

Two  OH  emission  sources  (W28S  and  W49I)  coincide  with  local  maxima  in  the  thermal  con¬ 
tinuum  within  the  error  limits  of  the  positional  measurements.  These  two  OH  sources  may  be 
inside  dense  H  II  regions.  However,  most  Class  I  sources  are  displaced  at  least  a  few  arcmin- 
utes  from  nearby  strong  thermal  continuum  maxima.  There  is  probably  some  association  be¬ 
tween  Class  I  OH  emission  sources  and  H  II  regions,  but  such  an  association  cannot  involve  ex¬ 
act  positional  coincidence. 

A  similar  statement  may  be  made  concerning  OH  emission  sources  and  supernova  remnants. 
In  the  Class  11^  sources  W28  and  W44,  the  OH  emission  lies  near  the  edge  of  an  SNR.  From  the 
analysis  of  Sec.  IV-B  we  conclude,  however,  that  there  is  probably  no  causal  association  between 
SNR's  and  OH  emission.  Also  many  Class  I  OH  emission  sources  have  an  SNR  nearby  but  prob¬ 
ably  not  close  enough  for  any  interaction. 

We  noted  in  several  regions  of  the  sky  (W43,  W44,  and  others)  that  we  are  probably  observ¬ 
ing  two  or  more  separate  sources  close  together  in  angle  but  probably  separated  in  distance.  If 
this  separation  is  correct,  then  interferometers  or  larger  antennas  should  resolve  these  sources 
in  angle  as  well.  In  some  regions  (W49,  NGC  6334,  and  W75)  we  are  clearly  seeing  two  inde¬ 
pendent  sources  close  together  on  the  sky  and  also  spatially  proximate  to  each  other. 

These  seemingly  unrelated  facts  may  be  understood  in  terms  of  a  breeding  ground  theory. 
Certain  regions  of  the  galaxy  may  contain  the  right  combination  of  density,  chemistry,  and  mag¬ 
netic  field  to  promote  the  formation  of  OH  emission  objects  and  also  supernovae.  Certain  dusty 
H  II  regions  are  probably  such  breeding  grounds. 

In  Sec.  VI-C  an  OH  source  picture  was  presented  which  indicated  that  at  least  the  Class  11^ 

OH  sources  are  probably  old  objects  because  the  simplest  explanation  of  their  spectrum  involves 
an  expanding  shell.  An  expanding  circumstellar  cloud  probably  corresponds  to  an  old  object  al¬ 
though  a  young  star  just  turned  on  might  be  seen  in  the  process  of  blowing  away  its  shell.  It  is 
interesting  to  note  that  the  same  class  of  OH  sources  (11^)  that  seem  most  likely  to  be  old  are 
also  not  closely  associated  with  H  II  regions  or  with  the  galactic  plane.  It  is  possible  that  these 
objects  originally  formed  in  the  same  breeding  grounds  as  Class  I  sources,  but  have  since  es¬ 
caped.  There  is  little  evidence  that  any  of  the  OH  sources  trace  spiral  arms,  but  this  needs 
more  research. 

This  study  has  not  solved  the  problem  of  the  OH  emission  sources;  this  report  can  only  de¬ 
scribe  the  progress  to  date.  The  area  of  interest  has  widened  to  include  some  optically  visible 
objects  and  the  recently  discovered  water  vapor  emission.  Rather  than  elucidating  the  situation, 
this  study  seems  to  have  uncovered  new  anomalies  and  new  complexities.  Perhaps  our  present 
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situation  with  respect  to  the  OH  emission  sources  is  similar  to  the  early  days  in  the  study  of 
stellar  spectra.  If  so,  we  have  a  great  deal  of  work  yet  to  do.  The  following  section  describes 
some  of  the  work  which  seems  most  likely  to  yield  new  insight  into  these  problems. 

B.  Suggestions  for  Future  Work 


1  Extragalactic  OH 


Recently  Roberts  (1967)  and  Radhakrishnan  (1967)  reported  negative  results  in  their  searches 
for  OH  emission  from  external  galaxies.  This  section  summarizes  calculations  on  the  detect¬ 
ability  of  extragalactic  OH  in  emission  and  absorption. 

Roberts  (1967)  based  his  calculations  on  so-called  "normal"  OH  emission,  that  is,  emission 
from  extended  regions  of  OH  with  excitation  temperatures  of  the  same  order  as  the  kinetic  tem¬ 
peratures  (Roberts  actually  calculates  upper  limits  for  the  OH  to  HI  density  ratios  on  the  as¬ 
sumption  of  equal  excitation  temperatures).  But  such  "normal"  emission  in  our  galaxy  has  only 
recently  been  detected  (by  Heiles,  1968)  and  is  very  weak.  Let  us  assume  that  the  external  gal¬ 
axy  M31  is  basically  similar  to  our  own  galaxy,  and  calculate  the  detectability  of  OH  absorption 
against  the  nucleus  of  M31  using  the  figures  for  the  strong  absorption  in  Sag  A.  We  pick  M31 
because  it  is  the  nearest  spiral  galaxy. 

5  4 

We  take  the  distance  to  M31  to  be  6.3  x  10  pc  and  the  distance  to  Sag  A  to  be  10  pc.  The 
OH  absorption  in  Sag  A  is  about  80° K  in  the  Parkes  210- foot  antenna  [Bolton,  et  aL  (1964a  and  b)], 
and  is  several  hundred  kHz  wide.  If  we  simply  scale  this  to  the  distance  of  M31  and  to  the  area 
of  the  NRAO  140-foot  antenna  we  get 


\6.3  x  10  / 

As  an  aside,  note  that  we  get  the  same  number  if  we  use  the  results  of  Goldstein,  et  aL  (1964). 
They  measured  about  6.5°K  of  absorption  with  the  Agassiz  60-foot  telescope.  We  have  avoided 
calculations  based  on  the  measured  continuum  temperature  of  M31  and  assumed  opacities,  be¬ 
cause  it  is  difficult  to  estimate  the  fraction  of  the  measured  continuum  temperature  that  arises 
in  the  nucleus. 

To  have  a  figure  to  compare  with,  we  calculate  the  theoretical  rms  temperature  fluctuation 
for  a  200° K  system  temperature,  a  100-kIlz  bandwidth,  and  an  8-hour  integration  time  as 


=  0.  008  °  K 


The  upper  limit  placed  by  Roberts  (1967)  on  the  emission  (and  presumably  also  the  absorption) 
in  the  direction  of  the  nucleus  of  M31  was  0.04°K,  based  on  a  45  minute  integration  with  approx¬ 
imately  the  above  system  parameters. 

Several  comments  need  to  be  made.  It  is  clearly  possible  to  increase  the  bandwidth,  per¬ 
haps  to  500  kHz,  and  also  increase  the  integration  time,  so  the  signal  is  detectable  if  the  above 
figures  are  accurate.  Also  both  transitions  can  be  measured  simultaneously  and  this  increases 
the  detectability.  However,  the  integration  time  requirements  depend  on  the  fourth  power  of  the 
distance  to  M31  and  this  distance  is  not  too  accurately  known.  And  finally,  we  are  looking  at 
M31  in  a  direction  well  out  of  the  plane  of  that  galaxy.  If  the  absorbing  OH  is  confined  to  the 
plane  of  the  galaxy,  as  seems  likely,  then  the  absorption  in  our  direction  may  be  much  less  than 
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that  between  the  sun  and  Sag  A.  It  is  possible  to  find  other  spiral  galaxies  that  are  seen  edge-on, 
but  they  are  considerably  farther  away  and  the  distance  dependence  mentioned  above  puts  them 
out  of  the  question. 

Now  let  us  assume  the  existence  of  a  W49-type  object  in  M31.  W49  produces  about  a  100°K 

peak  in  the  NRAO  140-foot  telescope  [Palmer  and  Zuckerman  (1967)].  If  we  take  the  distance  to 
W49  to  be  15kpc  and  scale  the  temperature  to  the  distance  of  M31  we  get  0.057°K.  However, 
this  peak  temperature  will  only  be  gotten  if  we  use  a  1-  or  2-kHz  resolution  and  this  increases 
the  rms  temperature  fluctuations  by  a  factor  of  7  to  10  for  the  same  integration  time.  Further¬ 
more,  we  know  neither  the  position  nor  the  velocity  of  this  hypothetical  W49-type  object  in  M31, 
so  a  search  is  necessary.  And  finally  we  should  notQ  that  W49  is  very  probably  unique  in  our 
galaxy;  its  spectral  luminosity  is  considerably  higher  than  any  other  known  Oil  emission  source; 
there  may  be  no  comparable  object  in  M31.  It  seems,  then,  that  a  search  for  W49-type  OH 
emission  objects  in  M31  is  less  likely  to  succeed  than  a  search  for  OH  absorption  against  the 
galactic  nucleus  of  M31.  Ball  and  Zuckerman  have  proposed  to  use  the  NRAO  140-foot  telescope 
to  look  for  extragalactic  OH. 

1 8 

2.  O  H  in  Emission 

16  18 

In  the  terrestrial  environment  the  isotopes  O  and  O  are  in  the  ratio  of  about  490  to  1. 

Barrett  and  Rogers  (1964),  Rogers  (1967a),  and  Rogers  and  Barrett  (1966)  calculated  the  fre- 

1 8 

quencies  of  the  four  18-cm  transitions  in  O  H  and  obtained  the  values  shown  in  Table  XVIII. 

1 8 

They  then  observed  the  1639-MHz  absorption  of  O  H  against  the  galactic  center  and  confirmed 
the  line  rest  frequency  and  the  relative  abundances.  Similar  results  have  been  reported  from 

Australia  [Robinson  (1967)],  and  by  Wilson  and  Barrett  (1968). 

18  1 6 
The  strength  of  radio  emission  from  O  II  is  more  difficult  to  calculate  because  the  O  H 

emission  sources  are  very  far  from  thermodynamic  equilibrium.  We  obtain  the  1/490  ratio  for 

the  018II  to  016II  radio  emission  only  if  the  OH  radio  emission  comes  from  a  maser  for  which 

(a)  the  terrestrial  oxygen  abundance  ratio  (1/490)  holds  for  the  observed  OH,  (b)  the  pumping 

mechanisms  are  independent  of  isotopic  species,  and  (c)  the  maser  at  both  the  0*hII  and  the 


TABLE  XVIII 

o18hand  o16h  frequencies 

O^H  Frequency* 

1 8 

O  H  Frequencyt 

F  -F 

(kHz) 

(kHz) 

1  -2 

1612231 

1584330 

1  - 1 

1665401 

1637460 

2-2 

1667358 

1639460 

2  -  1 

1720533 

1692630 

*As  measured  by  Radfard  (1964). 

|As  calc 

ulated  by  Barrett  and  Rogers  (1964)  with 

expected  errors  af  about  100  kHz. 
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1 8 

O  H  transitions  is  fully  saturated.  Assumption  (a)  is  likely  to  hold,  at  least  approximately; 
however  (b)  and  (c)  are  less  likely.  Some  of  the  proposed  maser  pumping  models  are  sensitive 
to  the  isotopic  species  or  at  least  the  isotopic  density.  For  example,  models  involving  selective 
absorption  or  resonance  trapping  of  pump  photons  [Litvak,  et  al.  (1966),  Litvak,  et  al.  (1969)] 
will  be  less  efficient  for  the  less  abundant  species  unless  the  absorption  lines  at  the  pump  fre¬ 
quencies  overlap  due  to  Doppler  broadening.  Pumping  models  which  involve  the  formation  of 
OH  in  excited  states  are  probably  independent  of  isotopic  species.  It  is  difficult  to  evaluate  the 
likelihood  of  assumption  (c);  however  the  strong  polarization  of  the  1665-  and  1667-MHz  emitters 
(e.g.,  W3  and  W49)  is  probably  related  to  a  maser  saturation  effect.  The  1612-MHz  emission 
from  NML  Cygnus  [Wilson  and  Barrett  (1968b)]  and  VY  Canis  Majoris  [Eliasson  and  Bartlett 
(1969)]  is  less  strongly  polarized. 

1 8 

Wilson  and  Barrett  (1968a)  reported  negative  results  in  their  search  for  O  H  emission  at 
1637  MHz  from  W3  and  W49,  and  gave  an  upper  limit  significantly  better  than  the  l/490  terres¬ 
trial  oxygen  abundance  ratio  times  the  1665-MHz  emission.  Because  the  1612-MHz  emitters 

1 8 

probably  belong  to  a  different  class  of  objects,  O  H  emission  at  1584  MHz  from  NML  Cygnus 
and  YY  Canis  Majoris  at  velocities  corresponding  to  their  1612-MHz  emission  should  be  sought. 
Such  a  project  is  under  way  (Ball  and  Penfield,  to  be  published). 

3.  High  Resolution  Stokes  Parameters 

Although  most  Oil  emission  features  seem  to  be  resolved  with  1-kHz  instrumental  resolution, 
it  was  pointed  out  in  Sec.  Vl-C-1  that  many  features  show  a  variation  of  polarization  properties 
as  a  function  of  velocity  across  the  feature.  This  suggests  that  some  of  the  emission  which  ap¬ 
pears  to  be  only  partially  polarized  may  contain  a  higher  degree  of  polarization  and  additional 
complexities  when  examined  with  higher  spectral  resolution.  For  example,  it  may  be  possible 
to  understand  the  unpolarized  features  in  W51  either  in  terms  of  overlapping  polarized  features, 
or  as  cases  in  which  the  variation  of  the  linear  polarization  position  angle  is  so  rapid  that  it  is 
smeared  by  the  1-kHz  instrumental  resolution. 

4.  Further  Studies  of  Time  Variations 

It  was  pointed  out  in  Sec.  VI-C  that  studies  of  the  time  variations  may  help  us  understand 
both  the  kinematics  and  the  physics  of  the  OH  emission  sources.  It  is  particularly  important  to 
determine  whether  the  variable  features  change  their  flux  without  changing  their  degree  and  sense 
of  polarization  as  suggested  by  Ball  and  Meeks  (1968)  and  Zuckerman,  et  al.  (1969)  or  whether 
the  polarization  percentages  also  change  as  suggested  by  Coles,  et  al.  (1968)  and  Goss  and 
Robinson  (1969). 

Another  important  area  of  study  is  now  possible  for  those  OH  emission  sources  associated 
with  visible  objects.  Any  time  variations  in  the  OH  emission  should  be  compared  with  variations 
at  optical  and  infrared  wavelengths.  Such  a  comparison  along  with  long-baseline  interferometer 
data  may  eventually  allow  us  to  specify  more  precisely  the  physical  conditions  near  the  OH  emis¬ 
sion  source  in  a  complex  object  such  as  VY  Canis  Majoris  which  apparently  consists  of  several 
stars  and  extended  nebulosity  [Wallerstein  (1958  and  1969)]. 

Ball  and  Zuckerman  have  proposed  to  use  the  NRAO  140-foot  telescope  and  the  413-channcl 
digital  correlator  to  measure  the  Stokes  parameters  of  well  known  sources  with  high  resolution 
in  order  to  look  for  the  effects  mentioned  in  Sec.  VII- B- 3  and  to  make  comparisons  with  previous 
work  as  a  check  on  the  character  of  the  time  variations. 
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5.  Search  for  New  OH  Emission  Sources 


Evidence  was  presented  in  this  report  for  the  existence  of  some  association  between  OH 
emission  sources  and  (1)  supernova  remnants  (or  galactic  nonthermal  sources),  (2)  H  II  regions 
(or  galactic  thermal  sources),  (3)  IR  stars,  and  (4)  the  galactic  plane.  One  should  therefore  look 
for  new  OH  emission  sources  near  objects  of  this  kind  which  have  not  already  been  observed. 
Several  such  programs  are  already  in  progress  or  proposed:  (1)  near  SNR's  or  galactic  NTS's 
by  Ball  (to  be  published)  and  by  Downes  (private  communication),  (2)  near  H  II  regions  or  galac¬ 
tic  thermal  sources  by  Downes  (private  communication),  (3)  near  IR  stars  by  Eliasson  and  Groth 
(1968)  and  by  Wilson  and  Barrett  (private  communication),  and  near  selected  optical  stars  by 
Ball  (to  be  published),  and  (4)  selected  sections  of  the  galactic  plane  by  Ellder,  et  al.  (1969). 

It  was  pointed  out  in  Sec.  VI-A  that  some  of  the  apparent  phenomenological  associations  be¬ 
tween  OH  emission  sources  and  other  objects  may  be  the  result  of  observational  selection  effects. 
To  check  this,  a  search  for  OH  emission  sources  in  randomly  selected  regions  of  the  sky  should 
be  conducted.  Only  if  such  a  search  turns  up  enough  OH  emission  sources  to  give  statistically 
significant  correlations  with  other  objects  can  we  be  certain  that  such  correlations  are  meaningful. 

6.  Interferometry  and  Correlations  with  H^O  Emission 

Very  likely  the  next  major  step  in  understanding  the  OH  emission  sources  will  come  from 
VLB  interferometry  (see  Sec.  VI-C)  and  from  a  study  of  the  H^O  emission  sources  [Cheung, 
et  al.  (1969),  Knowles,  et  al.  (1 969 )]  which  are  probably  the  same  objects  as  those  seen  in  OH 
emission.  VLB  interferometry  on  OH  emission  sources  is  being  pursued  by  Moran  (private 
communications),  and  Lilley,  Penfield,  and  Ball  have  proposed  to  use  the  Haystack- Westford 
interferometer  to  study  the  spatial  distribution  of  H^O  emission  sources. 
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